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Abstract: The IR dynamics of effective holographic theories capturing the interplay be- 
tween charge density and the leading relevant scalar operator at strong coupling are an- 
alyzed. Such theories are parameterized by two real exponents (7, 5) that control the IR 
dynamics. By studying the thermodynamics, spectra and conductivities of several classes 
of charged dilatonic black hole solutions that include the charge density back reaction fully, 
the landscape of such theories in view of condensed matter applications is characterized. 
Several regions of the (7, 5) plane can be excluded as the extremal solutions have unaccept- 
able singularities. The classical solutions have generically zero entropy at zero temperature, 
except when 7 = 5 where the entropy at extremality is finite. The general scaling of DC 
resistivity with temperature at low temperature, and AC conductivity at low frequency 
and temperature across the whole (7, 8) plane, is found. There is a codimension-one region 
where the DC resistivity is linear in the temperature. For massive carriers, it is shown that 
when the scalar operator is not the dilaton, the DC resistivity scales as the heat capacity 
(and entropy) for planar (3d) systems. Regions are identified where the theory at finite 
density is a Mott-like insulator at T=0. We also find that at low enough temperatures the 
entropy due to the charge carriers is generically larger than at zero charge density. 
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1. Introduction and summary of results 

The AdS/CFT correspondence has provided a new look at the connection of gauge theories 
and string theories. It also implements a new way to model strong-coupling physics using 
string theory and its low energy approximations, namely gravitational theories. 

Despite several attempts, the relevant (super)string theories, involving Ramond-Ramond 
(RR) backgrounds, have so far defied solution, although important progress in this direction 
has been achieved recently. Most of this progress is related to the AdS/CFT correspon- 
dence and its avatars, and relies on effective low-energy string theories, emerging when 
the strings can be approximated as point-like objects. In the most studied example, that 
of N=4 super Yang-Mills (sYM), this approximation becomes reliable when the 't Hooft 
coupling is large and the dual quantum field theory is strongly-coupled. The intuition 
developed in this limit suggests that the "gap" between the gravitational masses and the 
string excitation masses becomes large, making the gravity approximation of the string 
theory reliable. This translates in the dual theory as the existence of a large gap in the 
spectrum of anomalous dimensions. 

This new paradigm of strongly-coupled theories has been applied to understand the 
low-temperature dynamics of strongly-correlated electron systems (for recent progress, see 
pi |2L EjL H, |5| ) . The idea behind this approach is that of universality and non-trivial emergent 
behavior. We already know that the physics of condensed matter (CM) systems is very 
varied, and involves many examples when emerging behavior (effective quasi-particles and 
associated interactions) may be described by very different-looking theoretical models. A 
class of such systems exhibits strongly-coupled degrees of freedom, in the border with 



magnetism, high-temperature (T c ) superconductors, heavy fermion systems, graphene and 
others, often with a behavior that is at odds with the traditional Fermi liquids. 

The application of holographic ideas to such CM systems relies on the development 
and quantitative description of universality classes at low temperatures, which is often con- 
trolled by zero-temperature critical behavior (also known as quantum criticality, [|j). The 
selection criteria for the appropriate theories rest on the identification and characteriza- 
tion of appropriate phases via their thermodynamic characteristics, their phase transitions 
and critical behavior, and most importantly, their generalized transport coefficients which 
connect directly to what is measured by experiments. In particular, in this paper we will 
search for strange metal characteristics. 

A basic and important ingredient of realistic CM systems is the presence of a finite 
density of charge carriers. In holographic CFTs, this implies the existence of a conserved 
global charge, and the duality between the associated conserved global current and a (mass- 
less) bulk gauge field. Finite density ground states therefore correspond to asymptotically 
AdS gravitational solutions "charged" under the appropriate bulk gauge field. This iden- 
tification relies on an important assumption that real-life photons are weakly-coupled in 
strongly-correlated systems, a fact supported in most cases by data. 

The initial study of a holographic finite density system involved charged Reissner- 
Nordstrom (RN) black holes of Einstein-Maxwell gravity with a cosmological constant. 
Such solutions proved to be the simplest laboratory for the study of phase transitions at 
finite density, indicating an interesting phase structure, || 0. Moreover, a study of the 
fermionic two-point function in this background suggested the presence of fermionic quasi- 
particles with non- Fermi liquid behavior, B ^], as well as an emerging scaling symmetry 
at zero temperature, [flC|| , associated with the AdS2 symmetry of the extremal RN black 
hole. The addition of a charged test scalar led to superfluid phase transitions, [Jllj, [T^, 13], 
that upon weak gauging of the boundary U(l) symmetry led to superconductivity. 

The RN laboratory has however an important disadvantage from a CM point of view: 
it has a finite (and large) entropy at extremality. Such a fact is considered an anathema at 
zero temperature and seems to contradict the third law of thermodynamics. In particular 



this is a signal of instability, since the Breitenlohner-Freedman bound, [14, 15 1, for the IR 
AdS2 is more stringent than that for the UV AdS4. Furthermore, there are scalar operators 
that violate the bound and are expected to condense, destabilizing the RN saddle point, 

The next step in modeling strongly-coupled holographic systems at finite charge density 
is to include the leading relevant (scalar) operator in the dynamics. This is generically 
uncharged, and it drives the renormalization group flow of the system from UV to IR. 
The corresponding holographic theory is an Einstein-Maxwell-Dilaton system with a scalar 
potential. In the case of zero charge density such theories have been analyzed in the 
recent past as they capture some essential features of pure YM theory in four dimensions, 



1 16, 17, 18, |J|, £CJ] • I n the finite charge density case, solutions to such theories have been 



reconsidered in view of holographic applications, [21, 22, 23, |J|, 25, 26, [J7|, 28, 29, 3C, 31 



3|, |3|. 



The goal of the present paper is to describe a general framework for the discussion 



of the holographic dynamics of Einstein-Maxwell-Dilaton systems with a scalar potential. 
This is a phenomenological approach based on the concept of Effective Holographic Theory 
(EHT), in analogy with Effective Field Theories low energy approximations to QFT. This 
strategy has advantages and disadvantages. It allows a parametrization of large classes of 
IR dynamics and a survey of important observables. On the other hand, it is not obvious 
if concrete EHTs can be embedded in well-defined string theories, or if concrete classical 
solutions that are (mildly) singular are acceptable as saddle points of the dual CFTs, 



|3j, ||,g3, |0 



We will parametrize EHTs of the Einstein-Maxwell-Dilaton kind in terms of the IR 
asymptotics of the two scalar functions that enter in their two-derivative action: the scalar 
potential and the Maxwell (non-minimal) coupling. Firstly, we will investigate classical 
solutions |36| that can be eventually completed into bona-fide asymptotically AdS solutions. 
This will be done in two different directions: by analyzing exact solutions in special cases 
and by studying the IR (extremal) asymptotics of the solutions, in order to classify the 
low-temperature asymptotics of the associated dual quantum field theories. 

Secondly, using this holographic setup, we will compute the conductivity both at zero 
(DC) and non-zero (AC) frequency, which along with the finite temperature entropy and 
specific heat are important experimental observables of finite charge density CM systems. 
Indeed, some of the landmark experimental signals of strange metal behavior include re- 
sistivity linear in temperature and AC conductivity scaling as an inverse power of the 
frequency in some regimes. 

The setup of the paper is as follows: For the rest of this section we will introduce 
and motivate the effective holographic theories we will study as well as some important 
definitions, giving at the end a consistent summary of the results and outlook of our work. 
In sections 2 and 3, we will study EHTs' dynamics at zero and non zero charge density. 
Conductivity will be treated in section 4 including the case of DBI dynamics. Then in 
sections 6-9 we will proceed to study in detail the equilibrium and non-equilibrium dynamics 
of EHT solutions. Several complementary definitions and calculations can be found in the 
Appendix. 

1.1 Effective Holographic Theories 

The AdS/CFT correspondence is a correspondence between strongly-coupled adjoint the- 
ories and weakly-coupled string theories. It is expected to be valid more generally even 
when the dual field theory does not have an adjoint description. An important distinction 
is between "color" (involving adjoint degrees of freedom that are gauged) with prototype 
examples SU(N) YM or N=4 sYM, and "flavor" with prototype examples the Gross-Neveu 
model or flavor in QCD. In the first case all observables are singlets, and are described by 
closed strings in the dual theory. In the case of flavor, there are singlet and non-singlet 
observables. The singlet observables can be described by closed strings, the non-singlet 
ones by open strings. 

The relevant string theories live in diverse dimensions, and have generically non-trivial 
RR backgrounds. They are currently beyond our computational control, although this may 
change in the future. In the best studied example, N=4 sYM, the string theory at strong 't 



Hooft coupling A can be approximated by the effective supergravity theory as the curvature 
of the background is weak at strong coupling. The signal in the dual QFT is that at large 
A, there is a gap in the spectrum of anomalous dimensions. Indeed, BPS operators have 
dimensions of 0(1), while non-BPS operators (corresponding to string excitations) have 
dimensions of 0{\^). 

Truncating a string theory to a finite spectrum of low-lying states is the central point 
of the EHT approximation. In a large-N CFT it can be a controllable approximation when 
the spectrum of dimensions has a gap that can be made parametrically large. This is 
analogous to the existence of mass gaps that make integrating out massive states a good 
approximation in usual QFTs. In non-conformal cases the situation is more complicated. 
At the intuitive level one can argue that a truncation of the string theory spectrum can 
be a reasonable approximation if none of the states that have been omitted can become 
relevant in the UV or the IR. 

We can think about integrating-out massive string modes, in order to obtain an effective 
description of a small number of low energy modes. Whether this is a good approximation 
also depends on the particular classical solution that describes the dual theory ground state. 
If fluctuations have suppressed higher derivative interactions then this has good chances of 
providing reliable information on the dual dynamics. 

In a holographic string theory, string states that correspond to dual theory operators 
having non-zero couplings in the Lagrangian, will certainly have non-zero profiles in the 
string vacuum solution as their sources will be non-trivial. This may change however for 
string states that are non-sourced in the AdS boundary. This is decided by the dynamics 
of the theory, which determines their vevs, and therefore their non-trivial profiles in the 
string vacuum solution. Even if such vevs are non-trivial, neglecting them may not affect 
much the dynamics for other operators. 

There may be cases where spectrum (and interaction) truncations give reliable results 
although there is no small parameter to justify their neglect. There are two complementary 
examples of such a state of affairs. One is the technique of level truncation in string field 
theory, used in order to calculate the effective tachyon potential in tachyon condensation 



studies, [37]. In these, the level truncation gives convergent results although it is not really 
understood why. A complementary example from QFT is YM. This is a theory with a 
single parameter, Aqcd which is expected to provide a universal order of magnitude for 
vevs up to dimensionless numbers of order one. However, the phenomenological success of 
SVZ sum rules has provided unexpectedly small values for the relevant vevs, suggesting 
that (at least for some observables) truncating the spectrum of string theory states may 
be a good approximation. 

Our strategy will be to select a set of operators (or dual string fields) that are expected 
to dominate the dynamics, and parametrize their EHT in terms of a general two-derivative 
action. Higher-derivative corrections are expected if the interactions of other states cannot 
be neglected in the structure of the semi-classical vacuum solution. We will however assume 
that this is not the case for simplicity. 

The next question is to decide the set of string fields to keep in the discussion of the 
classical solution that should describe the ground state or the finite-temperature saddle 



points of the dual strongly-coupled field theory. The minimal case includes only the metric. 
This is always important to consider as: it always encodes the energy distribution in the 
ground state; it is always sourced by the boundary metric of the dual QFT. In this minimal 
case, the two-derivative action is AdS gravity, with a single acceptable solution at zero 
temperature, namely AdS p +i. 

In the case of a system that contains a global U(l) symmetry (particle number con- 
servation), a relevant operator that may be important is the associated conserved U(l) 
current J M dual to a bulk U(l) massless gauge field 1 A^. At the two-derivative level, the 
action is that of Einstein-Maxwell with a cosmological constant. If the ground state of the 
theory contains a finite density of the global charge, then the gauge field must be non-trivial 
and the relevant solution is the AdS-RN black hole. The Einstein-Maxwell description is 
relevant for U(l) symmetries that are acting on adjoint fields, like the R-symmetries of 
N=4 sYM. They are closed string fields and they originate in the bulk. 

There can be also U(l) symmetries that are generated by fundamental fields and are 
more like baryon number in large-N QCD. The origin of the associated U(l) gauge field is 
from flavor branes, and its effective action is given by the DBI action. The non-linearities of 
the DBI action implement interactions between charge carriers. Only in the dilute charge 
limit the action linearizes to a Maxwell action. In this context, the relevant charged solution 



is the AdS-DBI-RN solutions first discussed in [[3g, 39[. 

In many cases the important dynamics of strongly coupled theory is driven by a scalar 
leading-relevant operator. This is the case in YM, where the operator is Tr[F 2 ] when 
the number of flavors Nf ~ 0{1). On the contrary, in the Veneziano limit, Nf ~ 0(N C ) 
another scalar operator namely qq becomes important and it is the interplay of these two 
operators that determines the structure of the vacuum. 

In such cases the minimal set of important bulk fields to consider for the vacuum 
solution are g^ controlling the energy, A^ controlling the charge density and <f> (dual to 
the relevant operator) controlling the coupling constant and its running in the IR. 

For a CM system, <p maybe representing the strong interactions of the ion lattice, or or 
the effect of spins on the charge carriers. In a sense it could describe the glue responsible 
for inducing interactions on the charge carriers. If the charge carriers are dilute, we expect 
that they do not backreact on the surrounding system and the dynamics can be treated in 
the probe approximation: solve first the equations for the Einstein-dilaton system, and in 
the background of this solution we can then solve the gauge field equations. 

On the other hand, at sufficiently large doping, the back reaction cannot be neglected, 
and the full set of equations must be solved. 

As mentioned earlier, we have also two options for the dynamics of the carriers. If 
they correspond to adjoints, then they arise from the closed string sector and a Maxwell 
action may be sufficient. If they are fundamentals then generically they originate on flavor 
branes and we should consider the DBI action. In the latter case, at weak charge densities 
we expect also the charge dynamics to linearize and be expressed by the linear Maxwell 
action. 



1 The gauge field may acquire a bulk mass if the global U(l) symmetry is spontaneously broken or 
anomalous, like U(1)a in QCD. 



At the two-derivative level, the EHT in the Maxwell case is described in the Einstein 
frame by the action 



S = Sg + S MaX , Sg = MP' 1 J <P +l 



xv-g 



R-\{d4>) 2 + V{<j>) 



(i.i) 



S Max = -MV- 1 j dP +1 x^-g^F, v F^ 



Indeed, the most general two-derivative action containing a metric, a vector and a single 
scalar can be brought to this form via field redefinitions. 

On the other hand, in the fundamental charge case the linear Maxwell action, SMax 
must be replaced by the DBI action as follows 

S DBI = -MP' 1 f <P +1 xe 2k ^Z{<P) (j- det(g + e~ k *F) - J-detg ) (1.2 



1 I <P +1 xe 2k ^Z{<j)) (yj- det(s + e- k <t>F) - y 7 - det g 



Here g = e~ ®g a is the Einstein frame metric, and Z has been normalised such that for 
F 2 <C 1 the DBI action reduces to the Maxwell one. The parameter k encodes the dilatonic 
nature of the scalar field (/>. If (/> is the dilaton, k = — \ ~z\- 

The dimension p + 1 of the holographic bulk space-time has been left arbitrary here. 
Indeed, even for CM applications, p can in principle take several values. For layered systems 
that are effectively 1 + 2-dimensional, without any extra adjoint fields, we need p = 3 to 
describe them. For real 3+1 dimensional systems without extra adjoint fields we would 
need p=4. However in the presence of extra adjoint fields more space-time dimensions are 
holographically generated and p can take higher values (as in the case of N=4 sYM). 

1.2 The UV region 

It is standard in holographic theories to assume that there is a non-trivial (strongly coupled) 
fixed point in the UV. This guarantees a globally controllable behavior, and an appropriate 
UV boundary where the theory is defined. 

This assumption is however not necessary if our only interest is to understand the IR 
dynamics. There are two concurring reasons for this. The first is the Wilsonian definition 
of quantum field theory that requires a UV cutoff that can be anything, provided a full 
definition of dynamics is in effect at the cutoff energy. The second is that the non-trivial 
dynamical condition that defines the "vacuum" and effectively defines all expectation values, 
is imposed in the IR. Typically this is a regularity condition that that ties together the 
two linearly-independent solutions in the far IR. Even in mildly singular (and therefore 
acceptable) IR asymptotics such a regularity condition is in order. 

To compute therefore correlators in the cutoff holographic theory, the only thing that 
remains is to evolve the solution containing the single integration constant from the IR to 
the cutoff using the equations of motion and then at the cutoff to distinguish the source 
from the vev contribution. This is a priori non-trivial because although the source and 
the vev linearly-independent solutions start hierarchically different near the boundary, the 
eventually become similar and divergent in the IR. It is the linear combination that cancels 
the divergence and the proper and acceptable solution. 



However, one can start from the IR with the two linearly independent solutions and 
extend them toward the UV. If there is an intermediate scaling region then a cutoff there, 
can be used to separate the solution into a source and vev piece and eventually compute the 
correlator. Otherwise only the correlator at low frequencies and momenta can be computed. 

Note that the real parts of correlators typically suffer from scheme dependence because 
of renormalization choices whereas imaginary parts do not. 

In the rest of the paper we will pretend that there is an AdS completion to the ge- 
ometries we are discussing. As argued above this is not necessary for computing IR data. 
We will assume though that the potential decreases toward the UV signaling a good IR 
approximation. 

1.3 The IR asymptotics 

Another important question concerns the form of the two scalar-dependent functions in the 
actions above, the scalar potential V(cj)) and the gauge coupling constant function Z(4>). 

The potential must be non-zero in all realistic examples as generic scalar operators 
have non-trivial n-point functions. Because of this the attr actor mechanism that is playing 
an important role in spherically symmetric solutions in supergravity in the absence of the 
potential is not at work here 2 . 

What will be our main focus is the IR behavior of the functions, V, Z. Motivated by 
controlled examples from string theory we will assume a pure exponential behavior at the 
IR, 

lim V(<p) ~ e- s * . (1.3) 

(j>—¥<X> 

This is a generic behavior in the absence of IR fixed points with AdS p +i asymptotics. It 
turns out that for generic values of the exponent 5, the asymptotic behavior in ( |1.3[ ) is 



enough to capture the important IR features of the theory, as shown in [17, |2Q|| . In the 
same references it was also shown that there are "transition values" for 5, that we will 
review in the next section, where important IR properties change. In the neighborhood of 
such transition values, a finer parametrization of the asymptotics is necessary in order to 
capture all the ramifications of the IR behavior. 

In this work, we will take V(<p) = Vfoe *, remembering that this is an approximation 
valid in the IR. To enforce that our solutions and analysis is well grounded in the UV, we 
will demand that e ^ — > in the UV for our solutions. This will guarantee that simple 
modifications of the potential like V — > ^ + V$e~ ^ will lead to asymptotically AdS 
solutions without modification of the IR behavior. This will be the working assumption in 
this paper. As argued however above, our results are independent of the UV completion., 

We now focus on the gauge-coupling function Z{4>). Again experience from both gauge 
supergravity actions stemming from string theory and tachyon condensation indicates that 
we may take it to have exponential behavior in the IR: 

lim Z{4>) ~ e 7 *. (1.4) 



2 S. Kachru has suggested that even in this case one can place some constraints on the flows however. 



In particular, in U(l)s arising in the adjoint (closed string) sector Z diverges in the IR, 
|22[| , indicating that the charge sector is driven to zero coupling in the IR. On the other 
hand if the U(l) arises from the fundamental sector (open strings, D — D systems) then 
4> is the analogue of the tachyon and the general arguments of Sen, [37], indicate that 



Z — > in the IR and the gauge field is driven to strong coupling. This is triggered by 
tachyon condensation. It has been suggested, |I(], El]], that this is the mechanism for chiral 
symmetry breaking in holographic QCD. 

Once we have the appropriate action as above we can address the question of the 
ground state solution of the system both at zero and finite temperature, by solving the 
appropriate equations of motion. Our main preoccupation in this paper is to unravel the 
dynamics in the cases where the presence of charge degrees of freedom backreact on the 
neutral sector, g^ u ,4>. 

From the solution one can study thermodynamics and in the presence of multiple 
candidates for the ground state of the system, the presence and nature of phase transitions. 
Moreover, the correlation function of currents can be computed at low energy and from 
this the AC and DC conductivities can be computed. They are important observables of 
the system, and can be used to characterize different models. 

1.4 On naked singularities 

As is well known by now, solutions to EHTs at zero temperature have typically naked sin- 
gularities. From the GR point of view naked singularities are anathema and are therefore 
exorcized. In holography however it is known that not all naked singularities are unac- 
ceptable. Well-defined and controllable perturbations of N=4 sYM lead to (mild) naked 
singularities in the bulk. Such singularities are generically resolved by lifting them to higher 
dimensions or eventually by the inclusion of the stringy states. 

It was Gubser first who asked the question: when is a naked singularity in a solution 
of an EHT acceptable, [gU? One criterion that he introduced in scalar-tensor backgrounds 
can be phrased simply as follows: 

• A naked singularity is unphysical if the scalar potential is not bounded below when 
evaluated in the solution 3 . 

However, as postulated in ]T7| , Eo[ there are further requirements that must be im- 
posed in order to accept naked singularities in the context of an effective holographic 
setup. 

• A naked singularity describes holographic physics that is unreliable if the second- 
order equations describing the spectrum of small fluctuations around the solution are 
not well-defined Sturm- Liouville problems 4 . The spectrum around a given vacuum 
solution is calculated by studying small fluctuations. Such fluctuations, once taken to 
be eigenstates of the p-dimensional Laplacian, satisfy a second-order equation in the 
holographic coordinate. Such an equation has two linearly independent solutions near 

3 Note that we are using the opposite sign for the potential compared to pM. 

4 This has also been observed in tensorial perturbations of higher codimension braneworlds, |42|. 
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the AdS boundary. Typically only one is normalizable. The same should also hold 
true in the IR. In the presence of naked singularities it may happen that both solutions 
are normalizable. Then an extra boundary condition is needed at the IR singularity 
in order to have a well-defined spectrum. Such cases are deemed unreliable. 

• A naked singularity is unreliable if the minimal world-sheet of a string contributing 
to the calculation of Wilson loop expectation values descends arbitrarily close to the 
singularity. Wilson loops expectation values at large- N are saturated by fundamental 
string world-sheets that start at the loop at the AdS boundary and end somewhere 
in the bulk. The ending point is determined dynamically by minimizing the Nambu- 
Goto action, which equivalently minimizes the world-sheet area given the boundary 
Wilson loop. In common cases discussed in the literature for large Wilson loops, the 
endpoint of the string world-sheet is at a minimum of the string frame scale factor of 



the metric, [43]. We demand that this remains a finite distance from the singularity. 



Both the previous criteria, that we define to determine a repulsive singularity, resemble 
a state of affairs that has been studied in string theory in a different context, namely the 
linear dilaton vacuum of string theory in two dimensions. In that case, the effective string 
coupling constant, e , diverges on one side of space. The theory however can be rendered 
well-defined and computable if a suitable potential (cosmological constant) is added which 
raises a barrier in front of the strong coupling singularity. Any finite energy process can 
approach only at a finite distance from the strong coupling singularity and physics is 
therefore well-defined. This analogy is very close to what happens in Einstein-dilaton 



holographic models for YM, g§ [E7|, [1J, gg. 

In EHTs the last criterion is ambiguous unless an embedding to string theory is known, 
or other information is available. The reason is that the relevant metric scale factor is that 
of the string frame metric and this is related to the Einstein metric via a dimension- 
dependent exponential for the dilaton, <£. Unless we know the dilaton dependence of our 
bulk scalar <p we cannot calculate the string frame metric from the Einstein metric. 

There are two extreme situations: 

• The scalar <f> has nothing to do with the dilaton. In that case g® = g a . 

• The scalar <f> is the dilaton. In that case, taking account of the non-standard normal- 
ization of the kinetic term in (|1.3[) we obtain 



2 



C = eV^'j; (1.5) 

In the general case when <f> contains an admixture of the dilaton we have the relation 



9% = e- k4, 9l v , 0<fc< A /— . (1.6) 
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1.5 The conductivity 

An important observable in a theory with a conserved charge at finite density is the con- 
ductivity. It is defined as the response of the system to an external electric field, measured 
by the induced current proportional to electric field. This is a linear response function 
and in general is a tensor. Since the backgrounds that will be studied in this paper are 
homogeneous and isotropic, the conductivity a is parametrised by a scalar. 

a can be calculated by linear response theory using a Kubo formula from the two-point 
function of the fluctuations of the spatial components of the gauge field in the appropriate 
background. Here we will consider fluctuations with very long wavelengths compared to 
their frequency. This is a good approximation for optical CM experimental measurements of 
the conductivity. Therefore we will turn on a electric field with harmonic time dependence 
proportional to e lut . We will define the frequency dependent (AC) conductivity in terms 
of the transverse retarded current-current correlator as is standard, 

&{w) s g-W-o) (1 . 7) 

The DC conductivity on the other hand is the zero frequency limit of the AC conductivity. 



This limit may be ill-defined in several cases, [44]. A direct definition can be obtained by 
turning on a constant electric field, calculating the induced current and thus obtaining the 
ratio, in the linearized regime. This was discussed in detail in [21 1 for probe solutions using 



the DBI action. For massive charge carriers, a string drag calculation is also possible as 
was pointed out in pi] ]. We will come back to this in detail at sections 5 and 6. 

1.6 Setup and Results 

We consider a class of Effective Holographic Theories containing the metric, a U(l) gauge 



field, and a scalar (1.1). The action contains only second derivatives of the metric, and both 



the scalar potential and the gauge coupling function are parametrized by exponentials: 

Z(<j>) = e^ , V(4>) = -2Ae - ** . (1.8) 

These parametrizations are string theory-motivated and describe the IR asymptotics of the 
dynamics, which depend on two real parameters 7, 8. Requiring that the IR potential V((f>) 
vanishes in the UV will allow to promote IR solutions to asymptotically AdS solutions. 

We first study the bulk physics of uncharged solutions (see |17], ^0[ for p = 4) adding 
the dynamics of a test vector field. In particular, we provide power solutions for any 8, at 
zero and finite temperature. We develop criteria for the generic naked singularities that 



appear in the extremal solutions based on previous works, |34|, 17, 2C] and show that at 
zero temperature the following classification holds: 



< \8\ < * , 2 -,s . In this class of theories the extremal (T=0) solution has a good 
singularity and the spectra are continuous without a mass gap. At finite tempera- 
ture there is a continuous phase transition at T = + , whose order is given by the 
parameter 8 2 . At any positive temperature the system prefers to be in the black hole 
phase. Moreover there is a single black hole solution at any temperature. 
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\$\ = \ 7 _i) - This is a marginal case: the singularity is good and the spectrum is 
continuous with a mass gap. This is a case where the IR asymptotics of the potential 



can be further refined, 17]. At finite temperature there is a minimum temperature 
T m i n below which there is no black hole solution. At T m \ a the system undergoes 
a phase transition to the black hole phase. The order of the phase transition is 



continuous and depends on subleading terms in the scalar potential, |4q| ,|46l. 



2(p+2) 



t _ 1 s < 6 < \ 3( -i) • The IR singularity is stronger here but is still acceptable 
holographically. The spectrum is discrete with mass gap. The asymptotic behavior 
of the masses of the n-th excited state is m n ~ n. At finite temperature there is 
a minimum temperature T m i n below which there is no black hole solution. Above 
T m i n there are generically two black hole solutions, a small and a large one. The 
small is unstable while the large one is thermodynamically stable. At a temperature 
T c > T m ; n the system undergoes a first-order phase transition to the large black hole 
case. In special cases, there can be more than two black holes above T m ; n and more 
than one phase transition is possible. 



8 > \l 3 ) P _i) ■ Such cases are unreliable in the effective holographic description. 



* $ > J (n-i) • Such cases violate the Gubser criterion. They are not acceptable 
holographically. 

Always at zero charge density, we analyzed the spectra of probe vector fluctuations 
with the following conclusions: 

1. ^ + / _\\$i > I • When p > 3 or p = 3 and and the UV dimension of the scalar A < 1 
then the potential diverges both in the UV and the IR and the spectrum is discrete 
and gapped. This resembles to an insulator. 

On the other hand if p = 3 and and the UV dimension of the scalar A > 1 then 
the potential vanishes in the UV and the spectrum is continuous. This resembles a 
conductor. 

2. — | < ^ + / _jJ2 < f • The spectral problem is unacceptable and therefore the 
holographic effective field theory unreliable. 

3- j + (_d$2 < — \ i This is holographically acceptable and the same remarks apply 
here as in 1. 

An important observable of such theories in view of applications to condensed matter 
physics is the conductivity, both AC and DC. There are several techniques for calculating 
them depending on the theory. At low charge density they can be computed in the probe 
approximation, where the charge density does not backreact on the system. At finite 
charge densities which is the main preoccupation of this paper the backreaction cannot be 
neglected and therefore we must study the conductivity in the full setup. We will calculate 
AC conductivities at extremality (T=0). This is the regime uj/T 3> 1. The relevant 
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background solution is the extremal solution that will provide the leading power behavior 
of the AC conductivity as a function of frequency. 

The calculation of the DC conductivity is much trickier in the fully backreacted 
Maxwell case we study here. We will calculate it assuming massive carriers and we will 
perform a drag- force calculation as in pi] ]. As we shall see the "drag" DC conductivity 
calculation makes an interesting prediction in the case where the scalar operator is not the 
dilaton: the DC conductivity is tightly correlated to the entropy and therefore also to the 
heat capacity. In particular, in the planar case p = 3, the DC resistivity p verifies the very 
interesting relation 

p ~ S ~ C v , (1.9) 

that seems to be valid in the strange metal region. 

At small charge densities we can use the probe approximation around the solutions 
mentioned above to study the conductivities: 

• At the lowest densities the DC conductivity a = 1/p behaves as 

2(p-l) 7 a+2(p-3) 

p~T (p-i)^-2 . (1.10) 

In particular, in this dilute regime it is independent of the dilatonic nature of the 
scalar field, which can be parametrized by a parameter k in the passage between the 
Einstein and the string metric as 

9% = e- k fg%. (1.11) 

The conductivity can be a linear function for 

8 (p-2) , 

7 = 7linear = 2 " W^T) ■ (L12) 

As the zero charge density system never exhibits an entropy linear in temperature, 
we conclude that when linear resistivity appears, it is due to the presence of charge 
carriers. 



For higher charge density, 



i 2fc(p-i)a+4 

p= T 2-(p-D^ (J*) . (1.13) 

a 



In this case, the resistivity can never be linear. 

We have then analyzed several classes of exact solutions with Maxwell finite (charge) 
density [ |36| . For two families we studied the full class of solutions. The first is a two- 
parameter family with 7<5 = 1 which contains two branches in the 7, 5 plane. The second 
is a two-parameter family with 7 = 5. In both cases the solutions have arbitrary boundary 
data (temperature, scalar charge and charge density). Finally we have also analyzed a 
family of scaling (near-extremal solutions) for arbitrary 7, 5 where the scalar charge is 
related to the charge density. 
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We start with jS = 1 where the general solution is known |36|]. Black holes exist for 
5 2 < 3 and have zero entropy at extremality. We have analyzed the thermodynamics and 
phase structure of these solutions, which depend on the specific values of 5. We focus here 
on the canonical ensemble 5 : 

• < 5 2 < 1. There is a single branch of black holes for any positive temperature. This 
is like the analogous zero charge case. The temperature rises without bound (except 
when 5 2 = 1 where it has a maximal value) with the size of the horizon. At T = + 
there is a continuous phase transition whose order depends on 5 (either second or third 
order in this range). At any positive temperature the system thermodynamically 
prefers to be in the black hole phase. 

• 1 < S 2 < 1 -\ — t=. In this case, apart from the extremal solution, at any finite 
temperature below a maximal temperature T max , there are two black hole solutions, 
a small and a large one. At any temperature < T < T max , it is the small black 
hole that is thermodynamically dominant. It is also the stable one, while the large 
black hole is unstable. The transition at T = + between the extremal and the small 
black hole solution is continuous, with its order depending on 5. In particular, all 
orders of critical behaviour can appear. Above T max it is the extremal solution that 
dominates as it is the only one present. The transition is of zeroth-order and suggests 
that in a UV-complete setup there will be also another black hole solution (of the 
RN type) that will dominate this region, changing the transition to first-order, and 
if appropriately tuned to higher-order. This is however not visible with our simple 
IR potentials. 

• 1-| — y= < 5 2 < 3. In this case the extremal solution (T=0) is behaving differently. 

For 1 + -t= < S 2 < -±f — , the spectrum of the charged carriers is gapped and discrete 
when the UV dimension if the scalar operator is A < 1 and continuous otherwise. 
These systems appear as Mott insulators. In the rest of the range -^j — < 5 2 < 3 
the spectrum is gapless. 

At T > there is a single black hole solution at every non-zero temperature in the 
IR effective action. It is always unstable and its temperature diverges at extremality, 
while it vanishes when the black hole becomes large. It is always thermodynamically 
subdominant to the extremal solution. This situation is similar to the chargeless 
case with discrete spectrum. Indeed the black holes here correspond to the small 
black holes there. This correspondence suggests also the completion of this picture 
when one adds back the AdS asymptotics. In that case a large black hole branch 
will appear and will match to the small black hole case at some T m i n . Therefore for 
T > T m ; n there will be generically two black holes (plus the extremal solution). At 
some higher temperature a phase transition to the large black hole will take place. 
So the system will be a Mott-like insulator up to T c and then there will be a first 
order phase transition to a conducting phase. 



D The grand canonical ensemble is also analyzed in the text. 
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In the first two regimes the AC conductivity at extremality scales as 



- 1 *£3£iH- 



at 



The exponent is always larger than 5/3 in the region, < <5 < 1 + -4= and diverges 
5 2 = I+-75 (see Figj2l| for details). The system behaves as a conductor. The situation in the 
case, 1 + -js < S 2 < 3, is different: for S 2 < |(5 + \/33), the Schrodinger potential diverges 
near the horizon at extremality, and hence the system is insulating, while it becomes 
conducting away from extremality (see the left plot in Fig. 22); for 4(5 + \/33) < 5 2 < 3, 



the system is conducting with a continuum of states at extremality. 

The DC conductivity for massive charge carriers is depicted in Fig. 23. In the range 
< 5 2 < 1, the DC resistivity vanishes at zero temperature and then rises in the black hole 
phase. For a special value of 5 in this range the rise is linear at low temperatures. The 
validity of linear regime is determined by the IR dynamical scale, £. In the intermediate 
range 1 < 5 2 < 1 + —t= , the resistivity drops steeply with temperature in the unstable large 
black hole branch and rises slowly with temperature in the thermodynamically dominant 
small black hole branch. In this latter branch, the low-temperature resistivity is linear for 
a special value of 6. Again the validity of linear regime is determined by the IR dynamical 
scale. Finally, in the upper branch 1 + -y= < 5 2 < 3, the resistivity is decreasing with the 
temperature. 

In the special case S 2 = 1, the resistivity is finite at zero temperature, then increases, 
and finally diverges at a higher temperature signaling potentially a critical behavior. 

The second class of solutions analyzed in full detail are the charged solutions of the 
theory with 7 = 5. This class is special as it seems to be the only possible case on the 
7,(5 plane that have finite entropy at extremality (beyond AdS-RN black holes). These 
are solutions with the full set of parameters (charge density and temperature). Unlike 
the previous case, they may not be the most general solutions. We have studied their 
thermodynamics and transport properties and the results are as follows. 

• < 5 2 < 1. There is a single branch of black holes for any positive temperature. This 
is like the analogous zero charge and jS = 1 case. The temperature rises without 
bound (except when 5 2 = 1 where it has a maximal value) with the size of the 
horizon. There is no phase transition as T — > + . At any positive temperature the 
system thermodynamically prefers to be in the black hole phase. There is however 
a continuous transition from the charged black holes to the neutral ones at finite 
temperature, as the charge is sent to zero. 

• 1 < S 2 < 3. Apart from the extremal solution, at any finite temperature below a 
maximal temperature T max , there are two black hole solutions, a small and a large 
one. For all temperatures < T < T max , the small black hole is thermodynamically 
dominant. It is also the stable one, while the large black hole is unstable. Again, 
there is no phase transition as T — > + is approached for finite charge, but there is 
a first order transition to the extremal solution as Q — > + for finite temperature. 
Above T max , the extremal solution dominates as it is the only one present. 
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The zero temperature AC conductivity for all such black holes behaves as to 2 indepen- 
dent of 5. 



The DC conductivity for massive charge carriers is depicted in Fig,3£. In both ranges 
the resistivity is non-zero at zero temperature and has a regular low-temperature expansion 
in integer powers of the temperature starting with a linear term. Again the span of the 
linear regime is determined by the IR dynamical scale, I. 

We have also investigated near extremal solutions for the whole 7, 5 plane in the fully 
backreacted case. Such solutions agree for 7<5 = 1 and 7 = 5 with the near extremal limits 
of the solutions mentioned above. They exist for a special value of the scalar charge and 
the charge density. They have a single parameter, the temperature. 

Such solutions exist when 

2(p - 1) + pj 2 - 2j5 - (p - 2)5 2 > , 7 2 -75 + 2 >0 , 2 - 5 2 + 7<5 > (1.15) 
The condition above define the Gubser bound for the whole (7, S) plane, as shown in figure 



42. The spin-2 IR spectra of such solutions can be reliably calculated from the effective 



holographic action when 



2 + 6(6-y) + t_JL >1 (1.16) 



or 

(7 - 5) 2 



2 + 5(5 - 7) + ^ J - < 1 and w p u > 4(p - 1) (1.17) 

2{jp - 1) 

with w p u given in ( |8.12j ). The spin-1 spectra are always reliable except when 

2 + 5(5-7) + ^^ <1 and (7 - 5) 2 > ^^ (1.18) 

2(p - 1) p 



when the spectrum is reliable above a critical value of the charge density, (p. 18] ) . 

For 5 = 0,7 = — y ^_ 1 , they provide Lifshitz solutions at extremality with Lifshitz 
exponent z > 1. This suggests that generically Lifshitz solutions are generated by a scalar 
along a flat direction of the potential that drives the U(l) coupling constant to zero in the 
IR. 

Analysing the phase structure, we find that, depending on (7, 5), continuous phase 
transitions of all orders appear as zero temperature is approached. Furthermore, the near 
extremal solutions have simple scaling thermodynamics functions. The entropy and heat 
capacity in particular scale as 



5" ~ Cq ~ 7 1 2(p-l)(l-5(5-7))+(5-7) 2 (1-19) 

near extremality. The entropy is finite at extremality only if 7 = 5. When the exponent 
in (1.19) is negative this indicates an unstable black hole. In this region this black hole is 



never thermodynamically dominant. 
When 

(p-l)(5 2 - 1 2 -4)(8(l-p) + (5- 7 )( 7 (l+p) + 5(5p-7))) 

A(2(l-p) + (5- 1 )( 1 + 5(2p-3))) 2 { • ' 
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(p_l)( 7 _0)a <q 



2(p - 1)(1 - 5(5 -7)) + (5 -7) 2 
then the system has charged excitations with a mass gap and a discrete spectrum. In the 
special case of p = 3 one more condition is needed for this: that the UV dimension of the 
scalar operator dual to <j) is less than 1. Such systems behave as Mott insulators. In such 
cases there will a finite temperature phase transition from an insulating to a conducting 
phase with a continuous spectrum and no mass gap. 
The extremal AC conductivity scales as 

I 6(l-p) + (5-7)(p7+(3p-4)i5) I -, 
a(0j) ~ J 2(l-p) + (5-7)(7+(2p-3)5) | _ (1-21) 

The exponent may become negative in a region of parameters, but in this same region the 
thermodynamics of such near extremal solutions is unstable. 
The DC conductivity for massive charge carriers scales as 



p ~ 7 n 2(p-i)(i-5(«-7»+(«-7) 2 . (1-22) 



■p-i 



As mentioned earlier, when the scalar is not the dilaton (k=0), we obtain p ~ S^~ , and 
in p = 3 in particular p ~ S. These low-temperature scaling laws are expected to be valid 
throughout the 7, 5 plane. When p = 3, there are several cases that exhibit linear resistivity 
at low temperature. They have 7 given by y± 



y± = 35 + 2k± 2yjl + (5 + k) 2 . (1.23) 

As the chemical potential scales in that case with a higher power of the temperature this 
behavior is valid in a regime where the chemical potential dominates the temperature as 
in real strange metals. 

The leading scaling behavior of thermodynamic functions and conductivities is ex- 
pected to be universal and be valid for all solutions for this class of theories near extremality. 
Therefore, they provide a profile that might help select EHTs as descriptions of concrete 
CM systems. 

We have also compared the entropy due to the charge carriers and compared it to the 
entropy in their absence. We found that generically the entropy at finite charge density 
dominates that in the absence of carriers as discussed in detail in the end of section |8.5| . 
Only in the gapped regions the both entropies are subleading and 0(1) and therefore a 
comparison at this stage cannot be done. 

Another interesting behavior concerns the zero charge density limit of the various 
solutions, in the context of the canonical ensemble. This can be taken explicitly for the 
7 = 5 and j5 = 1 solutions and compared to the zero charge density solutions. 

We find the following: In the j5 = 1 case there seems to be no transition if we study the 
equilibrium thermodynamic functions. In the 7 = 5 case however there is a phase transition 
in the zero charge limit: it is second-order if 5 2 < 1 and first-order if 1 < 5 2 < 3. The 
behavior of the conductivities in this zero charge limit is more interesting. The behavior of 



the AC conductivity at extremality for the gap less charged cases is summarized in (1.21) 
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and the exponent is independent of the charge density. A calculation of the conductivity 
in the zero charge case, done in section Q gives instead for the conductivity 

a~cj n , n = ' 7 ' - 1 (1.24) 

1 — o z \ 



and this is generically different from (1.21|). This indicates that for all ~f,5, there is a 



transition as Q — > that affects the scaling behavior of correlation functions. This is 
also true in the case of the standard RN black holes. Similar remarks apply to the DC 
conductivity at low temperature both at finite and zero charge density. 

Finally we have found solutions to the bulk equations in the limit where there are 
strong self- inter actions in the charge sector. This corresponds to the limit where the DBI 
action vanishes. In this limit the charge density has a limit which is independent of inte- 
gration constants indicating that in flavor brane systems there is a saturation limit for the 
charge density, imposed by the DBI action. We have found two non-trivial solutions to the 
equations, one that is extremal, and another for p = 4, that is near extremal but does not 
seem to fit the general pattern of Maxwell solutions described above. This system is an 
insulator at zero temperature. Its DC resistivity for massive carriers behaves as p ~ T 



i± 



2' 



indicating linear resistivity when the scalar is not the dilaton (k = 0). 

1.7 Outlook 

The set of holographic effective theories discussed here, parametrized by the exponents 7, 5, 
is severely constrained by the nature of the IR singularity. Imposing constraints so that 
the singularity is of the good kind removes part of the 7, 5 plane from consideration. 

The solutions fully analyzed, namely the uncharged solutions as well as 7^ = 1 and 
7 = 5 cross sections, give a rather representative profile of the generic case. The near 
extremal asymptotics are summarized by the scaling solutions that can be found for any 
value of 7, 5. 

Calculations of the transport coefficients, namely the conductivity, are important in- 
gredients for the characterization of the associated physics. We have given estimates for the 
scaling behavior of both the DC and AC conductivities. They suggest that there are EHTs 
that display the hallmark behavior of strange metals, namely linear resistivity. Moreover, 
the DC resistivity is related in the planar case to the entropy, correlating two interesting 
aspects of strange metal behavior, linear heat capacity to linear resistivity. 

There are however more features of the physics of the EHTs that need to be analyzed 
so that their suitability as theories of strange metal behavior can be assessed. An important 
general ingredient are the zero and low temperature spectra of fluctuations, that should be 
derived in order to completely characterize the low energy degrees of freedom, as well as 
the energy-energy and charge-charge correlators. 

In particular the interplay between insulating versus conducting behavior must be 
further analyzed. It is a generic property of the holographic systems described here, in the 
planar case p = 3, to favor conducting behavior. Indeed, with the exception of strongly 
relevant dynamics, the systems at small charge density seem to be conductors at any 
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small but non-zero temperature. There is generically a continuous phase transition to the 

extremal solution which is insulating. 

The phase structure also seems to be interestingly varying with the IR "strength" of 

the scalar operator captured by the exponent 5. There are indications from our analysis 

that the appearance of discrete spectra at low temperature is "delayed" by the finite charge 

density. At zero charge density, such spectra appear when 5 2 > 1 while at finite density in 

the 7 = 5 case they never appear, and for j5 = 1 their appearance is seemingly delayed 

until 5 2 = 1 -\ — 7=. This structure needs verification from a more detailed analysis of low 
V3 

energy spectra. 

The underlying fermionic nature of the holographic system is also a subject that needs 
clarification. Fermi surfaces and related observables are part of the phenomenology of 
strange metals and their cross-over behavior to Fermi liquids. A study of fermionic response 
functions along the lines of |jsj, ||, 10 1 is necessary in order to characterize the fermi surfaces 



and the presence of quasi-particles. 

Last but nor least, the superconducting instability must be analyzed. For this the EHT 
used here needs to be extended by the inclusion of an extra charged scalar operator, X. 
At the quadratic level, and in the most general case, this introduces two extra functions of 
the neutral scalar into the model. 

2. The zero charge density dynamics 



This is a warm up case that has been studied extensively lately in [16, 17, 19, 2C] in view of 
mapping the landscape of Einstein-dilaton gravity with a potential, with the intention to 
use such a model as a phenomenological model for large-N YM. Our aim here is to obtain 
constraints on 7 and 5 in order to map the allowed region of parameters for the non-zero 
charge density exact solutions that will follow in the later sections. The action is given by 



(1.1) where we consider zero charge density and to start with a Liouville potential, 

V{<j>) = Voe~ s t Z(cf>) = 0. 

As argued earlier we will assume that the potential will be modified eventually in the UV 
so as to lead to asymptotically AdS solutions. In the domain-wall frame ansatz for the 
metric 

Ar 2 
ds l = e 2A(r) [_ /(r)dt 2 + dx i dx ^ + £!_ (2-1) 

f{r) 
the equations are 

(p-l)A" + i^ 2 = 0, g'(pA' + g')+g" = 0, f = e° , (2.2) 

(p - l)A'(g'+ P A') - l -4>' 2 - e-vVicji) = , (2.3) 

where primes stand for derivatives with respect to r. 

Since A'{r) is monotonically decreasing asymptoting near the UV, so is e . There 
are two exclusive possibilities: either e asymptotes to an AdS scale factor or it vanishes 
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somewhere, and this is a naked singularity, [jlffl. This property will persist also at finite 
density. The general solution to this system has been found in [E^], 

For our purposes here it suffices to study zero temperature solutions and therefore we 
set g = in ([2~2|,|2.3|). The background solution is 



V S 4 



e A = r(^, e 6 * = , " , r 2 , f = l. (2.4) 

2 (^r - * 2 ) 

Defining general coordinates, 

ds 2 = -D(r)dt 2 + £(r)dr 2 + C{r) (dx 4 dx*) , (2.5) 

the background solution (|2.4|) becomes 



4 TA j>4 

£)(r) = C(r) = r ^^ , Bfr) = 1 , e S(t> = —. — r- r 2 . (2.6) 

2p x 



This can also be written in "AdS-like"-coordinates, 

ds 2 = r 2 ( _ df 2 + dx 2 + dy2) + r (p-l)^-2 dr 2 j (2 7a) 

-r^ 1 ^ 2 . (2.7b) 



y o^ 4 ^(P^D^ 2 



Jl> _«J2 
p-1 

This is not AdS p +i and explicitly breaks its SO(p,2) symmetry group. This breaking 
corresponds to a non-zero value for 5, and goes together with a non-trivial scalar field 
profile. On the other hand, setting 5 = restores the SO(p,2) invariance and yields a 
constant scalar field. The potential V(4>) then is simply a constant, as expected. However, 
none of the analytic solutions presented below have both a non-trivial dilaton and AdS p +i 
asymptotics in the case of a pure cosmological constant. 

By going through a conformal transformation of the metric g^ u — v e « g^ u in appropri- 
ate coordinates, one can show that this space-time is conformally flat, recovering Poincre 
invariance. In the string case, 7 = 5= */ ^5, the coordinate transformation induces 
a logarithmic branch and then the background in the string frame is simply Minkowski 
space-time. 

To carry out a general analysis, in the domain wall frame (|2.l| ), it is convenient to 
introduce a "superpotential" W{4>) as 

,, = dt_dW ,, _ ch4 _ W 

- dr ~ d0 ' dr ~ 2(p - 1) ■ l j 



Then the system of equations ( |2.2 , |2.3 ) is equivalent to fl2.8|) and 



2 „ = J^!__«\ (2.9) 

2(p-l) V<W l ' 



Moreover, (2.2,2.3) and ( ^.§|J2.9| ) have the same number of initial conditions, namely three. 



Note also that each solution W of (|2.9|) provides two solutions to the equation of motion, 



one satisfying (2.8) and another where the sign of both equations in ( |2.8j ) is reversed. 
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Here we analyze the general solution of the zero-temperature superpotential equation, 
eq. ffOD , (below we set A = e*). 

-A 2 (W'(A)) 2 + ^-^H/ 2 (A) = 2V(A) , (2.10) 

To exclude the well known case of IR fixed points, that have been discussed extensively 
before, we assume that the potential is monotonic. It must be positive near the UV fixed 
point, and we will assume that it is positive, V(X) > 0, everywhere. 
First we observe some general properties: 



1. The solution can only exist as long as |W(A)| > d ' V(X). 

2. The equation has a symmetry W — > —W, so we can limit the analysis to W > 0. 



3. For any value Aq 7^ for the scalar, there are two solutions of ( 2.10| ), W + (X), W-(X) 



passing through the point Ao, such that W+(Ao) = W-(Xq), and W+(Ao) = —W'_{Xq). 
In other words there are two branches of solutions: one where W and W have the 
same sign (i.e. WAX) > 0) , another where they have opposite sign (Wl (A) < 0). 



4. At any A* / where |W(A*)| = v/^y^V^A*), the derivative vanishes, W'(\*) = 0. 

5. A solution can go past such a point A* only if V(A*) = 0. Indeed, suppose that 
V'(A*) > 0. If the solution exists for A < A*, at the point A* we have: 



W(X*) = J 4(p ^ V(\*), W'(X*) = 0, F'(A*)>0 (2.11) 



=► W(X* + e) < J 4( - P X V (A, + e) , 

therefore the solution does not exists for A > A* . 

6. By the same argument, if V'(A*) < 0, the solution does not exist for A < A*. 

7. It follows from points 3,4 and 5 that, if V(A) is positive and monotonic, the two 
branches W+(A) and W-(X) (see point 4) are completely disconnected, since neither 
W' nor W can change sign. However two solutions belonging to different branches 
can be glued together at a point W' = 0. 

8. All solutions that reach A = have either W(0) = yJ^SzUvffi), or W'(0) = 00. 

In what follows we assume V(A) > and without loss of generality we take W(A) > 0. 
An analysis on how solutions approach critical points for p = 4 is described in detail in 



the appendix of [20 1 and can be generalized to arbitrary p > 2. This behavior is summarized 



in Fig|l| for the W- branch. 



22 



W(X) 



W(X) 




Figure 1: Superpotential on (a) the W- branch and (b) the W+ branch, close to a critical point. The 
black area is the "forbidden" region below the critical curve \ I ' p ~ ' V, where W' would become 



imaginary. The solution stops where it meets the critical curve. 



2.1 Solutions close to A ~ oo 



We can analyze the solution of ( |2.10| ) in the asymptotic region of large A. We assume for 
the potential a power-law behavior 



V(A) ~ V&A 1 * 1 

A— >+oo 



where we took 5 to be negative so that the potential is increasing/ 
There are two kinds of asymptotic solutions: 

1 . a continuous one-parameter family of the form: 



W C (X) = CX ±I{ + 



2V 



C(-^Tk(\5\Tk) 



x \S\Tk + 



where C is an arbitrary constant of integration; 
2. a single solution that asymptotes as 



w s {\) = w 00 \- + ••• , 



w n 



8F 



/' 



2(p - 1) 



^-* 2 ' 



(2.12) 



(2.13) 



(2.14) 



Notice that if Vq > 0, both types of solutions exist only if \5\ < —^: for \8\ > -^j the l.h.s. 
of the differential equation is asymptotically negative. In this case there is no solution that 
reaches arbitrarily large values of A, but rather all solutions to ( 2,1C| ) are of the bouncing 
type: they reach a maximum value A* where a W + and a W- solutions join. 



6 The potential has to increase at the AdS boundary. If it decreases in some intermediate range there 
there is an IR fixed point. We are not interested in this case here. 
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2.2 General Classification of the solutions 

The general classification of solutions for the superpotential is as follows: for any positive 
and monotonic potential V(A) with AdS asymptotics in the UV and 

V(X) ~ FooA^I, V >0, 

X— >+oo 

the zero-temperature superpotential equation has three types of solutions, that we name 
the Generic, the Special, and the Bouncing types: : 

1. A continuous one-parameter family that has a regular expansion near the boundary 
and reaches the asymptotic large-A region where it grows as 



W ~ C b W a *- 1 i + 



A — > oo . 



(2.15) 



where C\> is an arbitrary positive real number. These solutions lead to backgrounds 
with "bad" (i.e. non-screened) singularities at finite values of the radial coordinate 
ro, where e — > and A — > oo as 



e A (r)~ \r-r \p, A(r)V 2 ("- 1 ) 



2(p - 1) 



(2.16) 



pCfe (?"o - r) ' 
We call this solution generic. An example for p = 4 is shown in the left part of Fig.|2| 




30 40 



Figure 2: From left to right, Superpotentials of the "generic", "special" and "bouncing" type. The black 



area is the forbidden region below the curve 



4(p-l) 



V. The AdS boundary is at A = 0. 



2. A unique solution, which also reaches the large-A region, but slower: 

W S (X) = W^ + ■ ■ ■ , ^oo = 



8F 



S-* 2 ' 



(2.17) 



We call this the special solution. An example is shown in the center part of FigS 
We can subsequently find the dilaton and scale factor to be 



,1*1* 



2 us-^ 



^o<5 4 (r - r ) 2 



I r — j-q I (p-i)s' 2 



(2.18) 
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3. A second continuous one-parameter family where W(A) does not reach the asymptotic 
large-A region. These solutions have two branches that both reach A = (one in the 

UV, the other in the IR) and merge at a point A* where W(A*) = y ~ V(\*). 



The IR branch is again a "bad" singularity at a finite value ro, where W ~ A V 2 (f " ; 
and 

e A (r)~|r-r p, A(r)~V^> ~ J^ - ^ . . (2.19) 

pC b (r - r) 

We call this solution bouncing. An example is shown at the right of Fig.|2| 

Notice that, as two solutions with positive derivative cannot cross, the special solution 
(Fig.0) marks the boundary between the generic (bad) solutions, that reach the asymptotic 



large-A region as Av 2 <p _1 ) (Fig.g) and the bouncing ones, that don't reach it, Fig.g. Notice 
that, if \5\ > — zf, only bouncing solutions exist. 

Finally, all solutions end up in a naked singularity in the IR. This can be seen as 
follows. We can evaluate the scalar curvature from the equations as 

(2.20) 
p — x a p — x [ ^{P ~ -U 

For IT 

n r _ «-l- 1 - - 

(2.21) 

When a < such solutions violate Gubser's criterion. For positive a > there is an 
IR singularity. 

2.3 The spectra of spin-2 fluctuations 

At zero temperature there are three types of perturbations, transverse-traceless fluctuations 
of the metric giving spin two excitations, a scalar excitation that corresponds to a gauge 
invariant combination of the trace of the metric and the scalar p| and a vector fluctuation 
originating in the gauge field. The equations they satisfy were analyzed in detail for p = 4 in 



R = - 


p-l 2 


p 
p — 1 


k 2 J +1 iaH 

4(p-l) 


C\ a in the IR 








R ~ 

p — 1 


a 2 - 


p + 1 
4(p- 


1). 


C 2 X 2a + ■■■ 



1 48|, |17[ . We will only consider the graviton fluctuation h^ u here as the second order operator 
is the minimal Laplacian. By factorizing h^y^^x 1 ) = h{r)h l _ iU {x l ) with D^^j, = m 2 h^ u we 
obtain for the radial wave-function 

h" + p A'ti + m 2 e- 2A h = . (2.22) 

Changing variables and coordinates to the radial coordinates, we obtain a Schrodinger 
problem 

7 rr , 2 , d-Z _A , (P-l) A p-l v (P-l) 2 SO 

ctr 2 4 

(2.23) 

where dots stand for derivatives with respect to z. The Schrodinger potential always 

diverges V s — > +oo near the AdS boundary. 
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We parameterize again the superpotential as W = C X a in the IR (large values of 
A = e^). As negative values of a violate Gubser's criterion we restrict our discussion to 
a > 0. We obtain 

A = w^W> log "' " r °' • ""* = cW^r) ' (2 ' 24) 



The Schrodinger coordinate defined in ( 2.23 ) is 



( r — rn ) 2( P -i) a 2 i 

- Zo+ i-Ub- ' A= 2(P-D.--i to "'-'° l + - ( " 5) 

2(p — l)a 2 < 1 . In this case the IR limit r — > tq corresponds to z — > oo. The 
Schrodinger potential from ( 2.23|) is in this case 



C 2 -C 1 

Sr 1 ' ^ = — — <°> ( 2 - 26 ) 



2 Z 



4 a 2 



w^i 



and vanishes near the IR singularity. The spectrum is therefore gapless and continu- 
ous. This is an acceptable system. 

2(p — l)a 2 = 1 . In this case z — zq = log |r — ro| and A = z — zq The Schrodinger 



potential from (2.23) asymptotes to a constant 



(p-D 2 
V s a ^-^- , (2.27) 

and the system has a gap and continuous spectrum above the gap. The spectral 
problem is well-defined and the system passes the Gubser criterion. 

2{p — l)a 2 > 1 . In this case the IR limit r — > tq corresponds to z — > zq. The 
Schrodinger potential from ( 2,23; ) is 



{Z ~ Z0) ^ (a 2 - ^L_) 

and rises steeply in the IR. This suggests a discrete spectrum. Near z = zq the two 
linearly independent solutions behave as 

^i~(z-z ) C , i>2 ~ (z - zo) 1 ^ . (2.29) 

ifil is always normalizable near z = zq as j z ° dz t/j 2 < oo. 

3 

2' 

; ' <«<a/P^. (2-30) 



When £ > I, implying 



2(p - 1) " V 6(p - 1) 



the second solution -02 is non-normalizable. This is an acceptable state of affairs as 
only one of the two solutions should be acceptable in the IR. In the opposite case, we 
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must impose an extra boundary condition in the IR and this is against the dictums 
of the holographic correspondence. 

For the remaining values of a, namely 



«>^' (2 ' 31) 

both solutions are square integrable in the IR and such superpotentials are therefore 
unacceptable. 

With the above results in stock, we now return to the solutions for the IR superpotential 



for a given potential described in 2.2. 



For any 6, the bouncing solutions have a < and violate therefore Gubser's bound. 



When a generic solution exists (\5\ < ? j^ ), a = \/ 2 (n-i) an< ^ for p > 1, this is 



always bigger than a / if A • Therefore, such solutions are always unreliable since 



they have a bad Sturm-Liouville problem for spin-2 fluctuations. 



Therefore only the special solutions (that exist for \S\ < \/ ^3[ ) with a = '-j- have a 



p- 

~2(p+2) 



chance of being physical and reliable. They are acceptable for 5 < y 377-rxy • 
In the acceptable range the spectrum of spin-2 fluctuations is as follows: 



< \5\ < ./ , ^ . In this class the spectra are continuous without a mass gap. 



\$\ = \ (p-i) ■ This is a marginal case: the spectrum is continuous with a mass gap. 
The asymptotics in this case can be further refined [o] but we will not pursue this 
further here. 



2(p+2) 



7— — pr < 5 < \ gj _J . Here the spectrum is discrete with mass gap. The asymp- 
totic behavior of the masses of the n-th excited state is m„ ~ n. 



$ > y 3 L • Such cases are unreliable from a holographic point of view, as justified 
earlier. 



• 5 > J , jU . Such cases violate the Gubser criterion. 

2.4 The spectrum of current excitations 

We now consider a small gauge field perturbation of the form A fl (r,uj,k)e' lu ' t+ , \i 



0, 1,2, ...,p — 1, in the gauge A r = around the neutral background fl2.q). The effective 
action for the gauge field reads 



. , . dudP-^k ZVDBCp- 1 

,s '- 2 - ' dr ' ~W9 2 — 



L/iri r^-7 rC^-/l'i rZj/i.i\ \K^J\.^\ CaJ-/!^ 



DB BC C 2 DC 

(2.32) 
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where we have used a general diagonal metric of the form ( |2.5| ) . The perturbation equations 
are 



AS + d r 



i««Zi/|cp- 3 



4 + ^Mi - Mo) - §(^ 2 ^ - **(£ • A)) = , (2.33) 
4)' + 5, log ^JU - ^A + ^-1=0, (2.34) 



along with the Gauss law 



y/DB 



Dr 7, 



C 



C 



toA' = -k-A 



(2.35) 



From ( 2.35| ) we may solve for the longitudinal component of A in terms of A$. the two 
independent degrees of freedom are therefore A± and Aq and satisfy 



A", + d r 



)ogZ { /^CP-* 



A\ + I -^- - —k 2 ) A i 



7C 2~ 

4,' + $.io g -= r 4,+ 



L> C 

Bw 2 Z? 



^ A 



0. 



V^B u V D C 
We first consider the IR asymptotics using the extremal chargeless metrics (12. 



(2.36) 

(2.37) 



D = C = r(p- 1 ^ 2 , 5 = 1, 



(2.38) 



and 



VoS 4 



li^-b 2 

V p-1 



Z = e 



70 



Z r 



2- 



ZJ 



W 4 



V p-l 



with the boundary at r = oo and the IR at r — > 0. In the IR, we obtain 

A » + (?! + 2{P ~Jl ) ~A' ± = -m 2 r~7^A ± , m 2 = uj 2 - k 2 . 



5 (p — 1)5 2 J r 



(2.39) 



(2.40) 



Wc now ddjn<< § -= rT^W and A ± = V z~f> , with /3 = (^ + ^^) ( (p-^iCa ) 



<ij 



to 



obtain a Schrodinger problem 



V + ^3-^ = mV , 



(2.41) 



with the potential diverging to +00 in the far IR because z — > as r — > when (p—l)5 2 > 2 
or with the potential vanishing to in the far IR because z — > 00 as r — > when (p— 1)<5~ < 
2. 

When (p — 1)J 2 > 2, the two independent solutions are behaving near z = as 



^4. ~ 2r , ip_ ~ z 



-/3+1 



(2.42) 



If /3 > I , one of the solutions, ip+, is normalizable in the IR and the other is not, while for 
j3 < — h , "0- is normalizable and ip + is non-normalizable in the IR. For — i < j3 < | , both 
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solutions are normalizable and this is holographically unacceptable. We deduce therefore 
that 

7 , (p-3) \ / (p-l)8i \ 3 (2 , (p-3) \ / ( P -l)5* \ _1 

* (P-1)«5 2 JI(P-1)* 2 -2J 2 ° r V^ (P-1)«5 2 / VCp- 1 )* 2 - 2 / 2 ' 

(2.43) 

In the UV, we assume the solution to become asymptotically AdS with D = B = 
C = =-% + • • • , r — y and Z ~ 1 + • • • . The subleading terms marked with dots in the 
previous relations are associated to the scalar perturbation of UV scaling dimension A and 
are analyzed carefully for p = 3 (where they are important) in appendix |F|. 

We can calculate therefore the UV asymptotics of the Schrodinger potential first by 
considering the leading contributions due to the UV fixed point 

where the coordinate r vanishes at the UV boundary. We observe that the potential 
diverges V± — > oo at the boundary when p > 3. For p = 3we must calculate the subleading 
contributions as 



U±, p=3 = --A(2A-l)r^^ + --- (2.45) 

where the parameter k is defined from the UV asymptotics of the gauge function Z as 



described in ( F.l 



Unitarity implies that 2A — 1 > 0. When A < 1, the potential diverges in the UV. 
Regularity then implies that k < 0, and the potential asymptotes to Vj_, p =3 — > +oo. If the 
perturbation is softer A > 1 the potential vanishes in the UV. 

From the above we conclude that, 



Condition ( p. 43 ) must be satisfied in order for the spectrum to be holographically 



acceptable. In this range, (3 (/3 — 1) > always and the the effective potential diverges 
to +oo in the IR. 

For p > 3 the potential is diverging to +oo on both the UV and the IR. The spec- 
trum is therefore discrete and gapped. This describes an insulator (although not a 
conventional one, as there is no charge density). 

When p=3, and the UV dimension of the scalar A > 1 then the potential vanishes 
in the UV and the spectrum is continuous. This resembles a conductor. 

When p=3, and the UV dimension of the scalar A < 1 then the potential diverges 
both in the UV and the IR and the spectrum is discrete and gapped. This resembles 
to an insulator. 



7 The relevant quantum mechanical problem in treated by Landau and Lifshitz, and a regularization is 
presented that chooses one of the two solutions when the regulator is removed. However the end result seems 
to depend on the way to regularize, and is equivalent to choosing a boundary condition at the singularity. 
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When (p — 1)5 2 < 2, the two independent solutions are behaving near z = oo as 



ip+ ~ z 



P 



if)-. ~ z 



-/3+1 



(2.46) 



If P > 2 , one of the solutions, ?/;_, is normalizable in the IR and the other is not, while for 
j3 < — 2 , "0+ i s normalizable and ip+ is non-normalizable in the IR. For — , < j3 < ^ , both 
solutions are non-normalizable and this is holographically unacceptable. Thus the same 
condition ( [2.43 ) still holds for (p — 1)5 2 < 2. In this range, j3 (/3 — 1) > always and the 
the effective potential vanishes to in the IR. Thus the spectrum is continuous for p > 3. 
This resembles a conductor. 

The spectrum of Aq is similar. The constraints above combined with the constraint on 
5 derived in the previous section are displayed in Fig.0. 
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Figure 3: The region allowed on the 7 — 5 plane from the spin-two and spin-one spectra is displayed in 
blue for p=3 and p=4 dimensions. 



2.5 The phase diagram at finite temperature 

At finite temperature, the boundary conditions at infinity require (Euclidean) time to be 
a circle with temperature /3. There are two solutions to (|2.2| ), (|2.3| ) that are relevant. The 
first is the background solution with time compactified ( p.4[ ) , for which the temperature is 
arbitrary. This is the thermal gas solution. The other is the (small) black hole solution in 
domain-wall coordinates as in (12.11) 



r (p-l)6' 2 ) 



with temperature 



V n 5 4 



\p-i 



4vrT 



2p 
(p-l)5 2 



1 )r, 



/ = ! 



(p-iM 2 




rcA(^T]^-l 



(2.47) 



(2.48) 
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For 1 5 1 < \ h^i the temperature vanishes at extremality ro — > 0. In the intermediate 



case, |<5| = W^rx ; the temperature is constant and independent of ro- Finally in the case 

1^1 > \/^TT the black hole has infinite temperature at extremality. 

Such solutions provide the IR asymptotics of solutions that are asymptotically AdS in 
the UV. For p = 4 the complete phase diagram was drawn in [20| for all S. The situation 
for other dimensions is qualitatively similar and we describe it below. 



In the first case, \5\ < \/^i , there is a single branch of asymptotically-AdS black 
holes qualitatively similar to ( p. 47 ). Once T > they have lower free energy compared 



to the thermal gas solution. Their free energy scales with temperature as 



W ~ T 2 -(p~ 1 ) s ' 2 , 

and hence the system generically shows continuous phase transitions to the dilatonic 
background at zero temperature. For the case p = 3, which is interesting from the 
condensed matter point of view, there is a n -order phase transition for the range 

n — 4 9 n — 3 

o<<* < 7> n = 4,5,6,... 

n — 2 n — 1 

In particular, the transitions are fourth-order or higher. 



In the marginal case \8\ = a/^tj , there is again a single branch of asymptotically- 
AdS black holes that are qualitatively similar to ( 2.47| ). Now however they have a 



minimum temperature given by ( 2.48 ): T > T m j n = ■?-. For T < T m ; n only the 
thermal gas solution exists. Once T > T m ; n the black hole solution has lower free 
energy compared to the thermal gas solution and therefore dominates. The transition 
at T = T m ; n is continuous. Its order depends crucially on subleading terms in the 
potential, pj| 



In the case \5\ > yipi , there are at least two branches of asymptotically- AdS black 
holes. In the generic case there are exactly two branches. One branch corresponds 
to "small" black holes that are qualitatively similar to ( |2.47 ). They are the only ones 



that appear as solutions to the IR potential. The large black hole branch appears 
because of the AdS asymptotics and is not visible when the IR potential is used. 
The large black-hole branch merges with the small black hole branch at some critical 
value of the horizon. There is also a minimum temperature T m i n . At T = T m [ n there 
is a single black hole solution. At T > T m j n there are two black hole solutions. The 
value of T m i n cannot be estimated from IR quantities alone but contains information 
from the full potential. 

For T < T m [ Q only the thermal gas solution exists and it dominates the canonical 
ensemble. Even above T m ; n the thermal gas solution is dominant up to a temperature 
T c > T m ; n . At T c there is a first order phase transition, to the large black hole phase 
that dominates for all higher temperatures. 
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T min 



Ah 



Figure 4: Temperature as a function of the value of A at the horizon, A^ for the three different cases 
discussed in the text. Black holes exist only above T min whose precise value depend on the 
particular zero-T geometry. 



In Fig.|| the temperature of black holes as a function of the value of the scalar at the 
horizon is plotted. The three different behaviors described above are evident in the plot. 

In [^] it was argued that in non-generic cases more than two black hole branches 
can appear. A case is plotted in Fig.|| where there are two local minima in the relation 
between black hole temperature and horizon position. In such cases multiple first order 
phase transitions can appear, as indicated in Fig.|| where the free energy is plotted as a 
function of temperature. A concrete example of a potential leading to such a case was 
given recently in [49|. 



3. The dynamics at finite charge density 

At finite (charge) density the gauge field is non-trivial. We will study solutions where the 
backreaction of the finite density on the geometry is fully taken into account. In this section 
we will assume that the gauge field is emerging from the closed string (adjoint) sector and 
we will therefore use a Maxwell-type of action. This however will also be relevant for 
fundamental-type conserved charges in the limit in which the DBI action linearizes (weak 
self- interaction of fundamental charges). 

We consider the Einstein-Maxwell Dilaton action fll.l] ). Since we are interested in the 
IR behavior we will take the scalar potential and gauge coupling constant to have their IR 
asymptotics that we will parametrize as 



Z(tf>) 



i4> 



v{4>) 



-2Ae" 



(3-1) 



We have used the freedom to rescale the gauge fields in order to set the overall constant 
coefficient in Z to one. This however implies that there is an overall unknown (dimensionful) 
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Figure 5: T as a function of A^ in a case with more than two black hole solutions. 

1 



Fb 




Figure 6: The free energy F(T) for the case of Fig.g with multiple extrema. Si and S2 denote small BHs 
whereas B\ and B2 denote big BHs. The arrows represent direction of increasing Ah- There are 
two first order phase transitions in this case. 



constant multiplying quantities like the conductivity. This constant can be fixed once the 
EHT is embedded in a string theory. 

This class of EHTs are parameterized by the two exponents 7, 5 that we take to be 
real numbers. There is an extra parameter, A , which sets the length scale of the solutions. 
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It will always appear combined with the asymptotic value of the scalar field 4>o, in the 
following linear combination, Ae ^° . Of course, once the asymptotic UV asymptotics are 
corrected to AdS, this scale is an independent (IR) parameter. 

Because the scalar potential in ( |3,1| ) has a runaway minimum when the value of the 
scalar field is infinite, it will not admit solutions with regular Anti-de Sitter asymptotics. We 
will require however from our solutions that they are correctable to asymptotically AdS p +i 
solutions. This will be done by requiring that when evaluated in the solution, V(4>) — > in 
the UV boundary. In this way, simply by adding a constant to the exponential potential the 
solution will be transformed to an asymptotically AdS p +i. More complicated completions 
are also possible. For example we can replace the potential by Vrn = — Acosh(5<j)) with 
a similar effect, although in this case the dimension of the perturbing scalar operator is 
different. Moreover, in this case the RN solution in p + 1 dimensions is also allowed, as was 
explored recently in p8| . 

We will also be interested in toroidal metrics for the transverse p — 1 dimensions. We 
will keep p arbitrary in parts of this paper while we will take p = 3 when discussing the 
7 = 5 and j5 = 1 solutions, as this seems to be a value mostly relevant for applications in 
CM. 



The equations of motion stemming from ( LI ) are the following 



R 



R + V 



pu 



-9iiv 




V„(v 



F 2 + V'(<f>) 



2 ] z{4>) 



F p F 



itajr 2 



l e ^F 2 + 25Ae~ 
4 



(3.2a) 
(3.2b) 
(3.2c) 



It is then easy to deduce from (3.2) the value of the Ricci scalar on-shell 



R 



\ ( W 



p+l yi ( P -3)Z f2 



1 



and in the domain- 



p- l 4(p 

The relevant Einstein equations in the general coordinate system (|2 
wall coordinate system of ( |2.1| ) are depicted in Appendix |A]. 

In domain wall coordinates fl2.ip the equation of motion for the gauge field yields 

q 



(3.3) 



A' t 



(3.4) 



Z(cf))e^- 2 ) A ' 
where we have adopted Minkowski signature and primes stand for derivatives with respect 



to the radial coordinate. There is a conserved Noether charge, [22], not to be confused 
with the U(l) charge of the black hole solution 

Q = e P A g'e°-qA t , ^=0. (3.5) 

It is conserved in view of equation ( |A.5| ). Evaluating it at the horizon (where At = 0) we 
obtain that Q ~ TS. Therefore when it vanishes, it implies extremality. 

The finite charge density solutions we will examine have been obtained for the first time 
in [pq], see (|6.1|) , or are generalisations of solutions found in previous literature [50, 51, 23], 
see Q and (|8T| ). 
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For the sake of completeness and to give an overview of previous literature in four 
dimensions, let us mention that: 

• Solutions for zero potential have been obtained by [^] for general gauge coupling and 
by [J53| for the string case 7=1, while |>4j shows solutions for particular values of 7 
with a hyperbolic cosine gauge function. Near-extremal black holes for these theories 



are studied in [g^, |56|, [57]], and in particular are shown to be non-asymptotically flat 
or AdS. 

In the case of Liouville potential, the general uncharged planar solutions were reported 
in [47] ; charged spherical black holes for the generic action (|1 . 1[) were obtained in 



while topological solutions were found in [BQ] for hyperbolic and toroidal black 



holes. Wiltshire et al. have studied and classified their asymptotics in [59, 30 1. All 



the above dilatonic black hole solutions with a Liouville potential are included in |36|] . 
Cases of black holes with AdS asymptotics are known analytically for particular 



potentials in hyperbolic sines or cosines and were studied in |3l], |52|, 28, 53 1 



The definitions of the various thermodynamic ensembles considered are given in appendix 



4. The AC conductivity 

In this section we will describe the calculation of the AC conductivity in the spatially 
homogeneous case both for the Maxwell theory and the DBI case, in view of applications 
for solutions with string charge interactions. 

4.1 The AC conductivity in the Maxwell case 

To compute the appropriate correlator, as we are working in the fully back-reacting case, it 
will be necessary to also perturb appropriate off-diagonal components for the metric. We 
will do our calculations for the general coordinate system where the metric is given in fl2.5|) 
and <p and At are functions of r only. At is given by 

4M = ^I. (4.i) 

We consider the following (small) perturbation to this background 8 

A i = a i (r)e iut , g ti = ^(r)e^ . (4.2) 

The first variation of the rx Einstein equation gives 

g'u - %9ti = -ZAiA't , (4.3) 



3 We can also perturb g r t instead of gti with similar results. 
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while the gauge field fluctuation equation gives 



., p-3 J D A f \ i „__ t p-3 /.B 2 A , 1 i ' C 



d r [ zc— x /-A' i + zc^J-co 2 A i +± [g'u - -^gu 



(4.4) 



Substituting ( [4.3[ ) into (4.4) we finally obtain a second order equation for 

a, = , 



/ P-3 /Z) , \ B-3 

9 r IzcVW—oJ] +ZCV 



^B 2 q 2 y/DB 



(4.5) 



where as usual primes stand for derivatives with respect to r. 

This equation can be eventually mapped to a Schrodinger-like equation by changing 
the radial coordinate to z and redefining the radial wave-function as 



dz 
dr 



<B(r) 

W) 



* 



p-3 

Z = ZC~ 



to obtain 



d^_ 

' dz 2 



+ Vy = LO Z V, 



V 



Vz' 

q 2 D t \{8 Z Z) 2 , l fl (8 Z Z 



ZCP~ ! '4 Z 2 + 2 ; 



z 



(4.6) 



(4.7) 



In many cases, for plotting the potential it is useful to express it in the original r coordinate, 
reading 



V(r) 



q 2 D D 
ZCP- 1 ~ B 



±(d r log zf - x -d 2 r log z - l -{d r log z)d r log (zJ^ 



■ (4-8) 



Near an AdS boundary, z is the standard radial coordinate with z = being the 



boundary, D — > p-, B 
near the AdS boundary 



Ti 



and Z asymptotes to a constant Z — > Z^. Then, 



V{z) 



Z b 



zx 2(p-2) (p_l)(p_3) 

+ 4^2 



+ 



(4.9) 



For p > 3 the potential diverges quadratically at the boundary and if there is a similar 
behavior in the IR the spectrum will be gapped and discrete. For p = 3 however we must 
calculate the subleading contributions as the potential above vanishes in the UV. These 
come from the scalar perturbation and are 



V± 



k 



p=3 



2A-2 



A(2A - \)r 1 ^ 1 + 



(4.10) 



where the parameter k is defined from the UV asymptotics of the gauge function Z as 



described in (|F.l 



Unitarity implies that 2A — 1 > 0. When A < 1, the potential diverges in the UV. 
Regularity then implies that k < 0, and the potential asymptotes to Vj_, p =3 — > +oo. If the 
perturbation is softer A > 1 the potential vanishes in the UV. In that case the spectrum is 
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gapless and therefore the system is a conductor. We conclude that strongly coupled doped 
systems in p=3 are generically conducting. 

If the UV fixed point is a Lifshitz point instead with Lifshitz exponent w > 1, then 
the UV asymptotics of the schrodinger potential are modified to 

The relation above is valid when w 7^ 2. It implies that UV Lifshitz asymptotics with 
w = 3/2 make the second term vanish for p = 2, and w = 5/2 make it vanish for p = 4. 
Therefore Lifshitz behavior affects the presence of a gap in conductivity in the IR. 

Finally when w = 2, the schrodinger coordinate is z = logr, and the second term in 
the potential ( 4.11] ) becomes a constant proportional to (p — 3) 2 and therefore vanishes at 



p = 3. The first term as usual asymptotes to zero at the boundary. 

Near a regular horizon at r = r^ we have the following expansions for the various 

functions 

D = c D {r - 7-0) + 0[(r - r ) 2 } , B = ° B - + regular , (4.12) 

[r - ro) 

while C and Z',Z asymptote to constants. From ( }4.6| ) the horizon is at z = —00 with 

/IE , 



r — ro — ev c s , and the potential vanishes exponentially there, 

q 2 c D 



V{z)-ih^I eVcB +••• ( 413 ) 

A)^o 

Finally near an extremal horizon with finite entropy, we have the following asymptotics 

9 cr 

D = c D {r-r y + ■■■ , B = - -r + regular , (4.14) 

(r - r ) z 

while C and Z 1 ', Z asymptote to constants. Then z — */— _ and the potential vanishes 
as 

z o <• 
In the case where Z is a constant a general formula for the AC conductivity was given 



m 

1 + 7c(cj) 

where TZ{uj) is the (energy-dependent) reflection coefficient of the Schrodinger problem 
above. In the case of a non-trivial gauge coupling constant Z, the associated conductivity 



was derived in [23] with the result 



/ \ 1 — TZ(uj) i Z ,. *-\ 

1 + 7c(wJ zwz () 

where the dot on the extra factor implies derivative with respect to 2 and this is evaluated 
at the AdS boundary. For regular boundary behavior of Z, this factor vanishes, and the 



conductivity is given by (4.1 
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Finally, in [34, 23 1, the zero-temperature asymptotics of the AC conductivity were 



calculated for a Schrodinger potential that asymptotes to 

V(z) = ^*, (4.18) 

near the extremal horizon, with the result 

<t(w) ~ uj 2 ^- 1 . (4.19) 

We may now use equations above to evaluate the extremal conductivity in the charge- 
less background solution ( p. 47 ). The Schrodinger coordinate in the IR is 



r (p-l)S^ 

-T— (4-20) 



1 — f 
and the effective potential ( |4.8| ) in the IR becomes 

V <" = # ' C = ((p - 1)^ - 2)2 (421) 

We observe that c > — ■j always. When 5 2 < — ^-, z — > oo in the IR and the potential 
vanishes there. This signals a continuous spectrum and we obtain the scaling of the AC 
conductivity as 

a ~ u n , n = v / 4cTl - 1 (4.22) 

There is a whole range of parameters where the exponent becomes negative. It becomes 

n = — 3 when 

1-5 2 2{l-5 2 ) /A . 

7 = -^^ or 7 = -^^ (4-23) 

In the opposite case, 5 2 > — j, z — > in the IR. In this case the potential diverges. 
When it diverges to positive infinity it indicates a discrete spectrum. In p = 3 an extra 
conditions is needed, A < 1. 

As we will see in the rest of the paper, neither the extremal AC conductivity nor the 
DC conductivity are continuous as the charge density q — > 0. 

4.2 The AC conductivity in the DBI case 

We will now consider the DBI action in the Einstein frame and turn on a A x (r, t) as well 
as a gtx{ r -,t) perturbation. 

The relevant effective action expanded to quadratic order in the perturbation is 



S DB i = -fdP +1 i 



e 2H >Z{4>)C^r x l <DB - e~ 2k ^F 2 t + 



1 zc iE^l DF 2 X - BF 2 X + 2g tx F rx F rt + Be^gj 

2 y/DB - e~ 2k <i>F' r \ 

Considering A x (r,t) = A x (r)e lult we obtain the fluctuation equation 
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(4.24) 



~d r ( -gQtx 



D 



d r e~ 2k<t> x —-, ,yje 2k <t>q 2 + ^ZHJp^A' + 



+u 2 e~ 2k ^ 



BC 2 

B 
DC 2 



^ q 2 e 2k^ + e ^ Z 2 C p~l Ax 



while from Einstein's equations we obtain 



g tx y/DB 
Combining the two we obtain for the gauge field the equation 



(4.25) 



(4.26) 



D 



o = d r [ e -^, ^^vTwc^i; + 



BC 2 



+ 



w 



B ^/q 2 e 2k ^> + e^ZtCP- 1 q 2 ^DB 



A,. 



DC 2 e 2k <>> Cp- 1 

In the strong coupling limit q S> e ^ZC~^~ the fluctuation equation becomes 



(J,- 



—A>] + 

BC 2 x 



aj 2 e 



B q 2 VDB 



DC 2 CP~ 



A x = 0. 



(4.27) 



(4-5 



The Schrodinger problem relevant for conductivity calculations that emerges from ( 4.2S| ) 
then is 

,7-, / D 



— - — 

dr~ ~ Vl>' 



Ce^ 



with the Schrodinger potential 



„ 3C' 2 k 2 , 2 kC', C" k., q 2 De k * 
^=4^ + T + 2C* ~ 2C ~ 2* + -CP= 2 ~ 



(4.29) 



(4.30) 



Near an AdS boundary the DBI action linearizes and the fluctuation equation becomes 

j 2 VDB 



-ZC— A' x + 



Or 



matching the equation derived in (|4.5|) 



o B p-3 



CP' 



-1 



Ar. = . 



(4.31) 



5. The DC conductivity 

For the DC conductivity we will explore several possibilities. In the Maxwell case we will 
do a drag calculation of the DC conductivity for massive carriers. In the DBI case without 
backreaction the calculation was first done in [21]. We will also amend this calculation 
below and apply it to the chargeless solutions of section [2L 
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5.1 The DC conductivity in the probe DBI case 

We consider again the DBI action ( p,2| ) where we will drop the second piece as it will not 
be relevant for the conductivity calculation. Keeping a general coordinate system, working 
in p + 1 dimensions and considering the ansatz for the gauge field At(r), A x (r, t) we obtain 
the action 



S D bi = - f <P +1 xe 2k,t> Z^- det(g + e~ k +F) 



I) 



d 4 xe 2k<t) Z C v ^r A BBC - e~ 2k <t> (cA' 2 + BA 2 - BA' 2 ) , (5. . 

and the gauge field equations of motion 

ir,% A'. 

= q, (5.2) 



ZCiAL 



DBC - e~ 2 W 



(CA' 2 



( 



t 



Z C—BA, 



V 



BBC - e- 2k * [CA' 2 + BA 2 - BA' 2 



+ BAl ~ DA x 
\ 
) 



( 



-0 r 



Z C^BA' 



\ 



\ J BBC - e- 2k <f> [CA' 2 + BA 2 X ~ DA' 2 



0. 



(5.3) 



/ 



Making the further ansatz A x = —Et + a(r) we obtain from ( |5.3D that the first term is 
identically zero and 



Z C^Ba' 



w . 



'BBC - e~ 2k< f> {CA' 2 + BE 2 - Ba' 2 ) 
with w, q constants. We may solve 

y/BCe 2 ^ - E 2 



(5.4) 



A'=±q 



a = ±w 



BB 



C y/(Bq 2 - Cw 2 ) + BCP^e^Z 2 ' 
BC ^BCe 2k ^ - E 2 



D ^/(Bq 2 - Cw 2 ) + BCP- l e 2k *Z 2 



and obtain the on-shell action as 



o 



on— shell 
DBI 



d A x Z 2 cVBBCe 2k ^- 



^jBCe 2k <t> - E 2 



^(Bq 2 - Cw 2 ) + e 2k <$>BCv- x Z 2 



(5.5) 
(5.6) 

(5.7) 



As first observed in 21] the square roots in numerator and denominator change sign 
between the horizon and the boundary and this must happen at the same radial point 

r = r*, 

B(n)C{n)e 2k ^ = E 2 , C(n)w 2 - B(n)q 2 = B{n)C(n)P- l Z 2 e 2k ^ , (5.8) 
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from where we obtain 

w *=(^. + C P-3 Z 2 e 2k^\ E 2 e -2h4>, (59) 



Near the AdS boundary, D = B = C ~ *y and 

q£ /r\P- 2 w£ /r\v~2 

At =» + ^rwSv + -"' a=ao+ j^wSv + " (5 - 10) 

We also identify 

(J t )=q, (J x )=w, (5.11) 

so that the (non-linear) conductivity is given by 



a 



(J»> 



9 2 + Cf _1 ^ 2 e 2 ^. (5.12) 



As first pointed out in |21| the first term under the square root is proportional to 
the charge density q and is therefore the contribution to the conductivity from the charge 
carriers in the system. The second term is independent of the charge density and is due to 
pair creation of charges in the medium. As pair production is expected to be Boltzmann 
suppressed at high carrier masses, this term is expected to be dominant for almost massless 
carriers, 

&pair — Ly* Zj* [p.LOj 

For heavier carriers the first term is dominant 

^a g - cdt (5 ' 14) 

and can also be calculated from a drag calculation as we will show in the next section. 
For the solution in ( 2.471) , equation (5J3) becomes 



8 i 4fc 

(p-l)*2^ 6 



2p 

r \ (p-i)s' 2 
r* 



2 te - ^ 2 



and the conductivity is given by 



4 | 4( 7 +fc) 



a = e -k<t>O u s (P-^y(Jt)2 +e 2( 7 +fc)0o r f ' * , (5.16) 

where we have taken Z = e 1 ^ . Consistency requires that the electric field is not very strong 
so that r* is fully into the IR region. 

When the electric field vanishes r* — > r$ and the conductivity becomes 



. 4 | 4(7+*;) 



and using 



ff = e -*«0r d ^ 1 ^y(jt)2 + e 2( 7 +fc)^o r ^ * j (5.17) 



47T 

r = {kTT , k =-- -j—^ r , a = — , (5.18) 



4^ - l) ' T^t^7 I 



WW~ l ) Jj ^ IW 
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we finally obtain 



2fc(p-l)<S+4 / 4(p-l)[l + (7+fc)i] 

a = e -^o( KT ) (p-d^-2 y(J*) 2 + e 2(7+fc)<^o( K T) 2-Cp-d^ , (5.19) 

For the resistivity at low density it is the second term in the square root that dominates 

I 2(p-l) 7 a+2(p-3) 

p = - ~ T (p-!)^- 2 , (5.20) 

cr 

while at higher densities it is the first 

1 2fc(p-l)<S+4 

p = - ~T 2-(p-d^ (J*) . (5.21) 

cr 

This result agrees with the drag-force computation in the next section as was first shown 



in ||1|. 



In the second case, taking into account the bound, 5 2 < J P _J derived in section ^^ it 
can be shown that the resistivity can never be linear. In the first case we can attain linear 
resistivity when 

6 2 _ 2(p-2) 

7 = Tlincar = J^ • (5.22) 

2d 
Note that in this dilute limit the resistivity scaling is independent of the nature of the 
scalar field, i.e. whether it is dilatonic in nature or not, which was parametrised by the 
number k. 

It should be stressed that a DC resistivity formula of the form (5.19) offer the possibility 



of interpolating between two different powers of the temperature as has been observed in 
cuprate materials where the resistivity is linear in the strange metal region and quadratic 
in the overdoped region. This will happen here if 5 = — |. We will also see in section 8.7 
that this can also arise in that case as well. 

5.2 The drag calculation of the DC conductivity 

Massive fundamental carriers can be thought of as endpoint of strings, as they end on flavor 
branes. For massive carriers the associated flavor branes are located away from the IR bot- 
tom of the bulk space. Their motion can be modeled as the motion of a fundamental string 
in the semiclassical background described by the dilatonic charged black hole solutions. 

We will consider here a massive charged particle moving in the appropriate background 
with a world-sheet action given in the Nambu-Goto form 



Sng 



T f f d 2 ^g + f drA^ , g a p = g^8 a x^x u , (5.23) 



where Tf is the fundamental string tension. The string has two end-points that each couples 
to the bulk An. For the drag calculation we follow one of the endpoints, as the second is 
very far away. £ a are the world-sheet coordinates while x^ are the space-time coordinates. 
We have also taken the Chan-Paton charge of the endpoint to be unity. 
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If the background is invariant under translations we have the following world-sheet 
Poincare currents 

vr£ = 7f« + A^rT = T f ^~gg^g vix dpx v + A^ , (5.24) 

with conserved momenta 

P„= I dairJ l + A lx = T f f da^gg^g^d^ + A„ , (5.25) 

where rj a/3 is the flat two-dimensional Minkowski metric. 

The bulk equations of motion stemming from the effective action for the coordinates 
for which the background is translational invariant are 

d a Tt« = , (5.26) 

while at the open string boundary we obtain 

T f ^gg aP g^dpx u + qF^x" = . (5.27) 

We now consider a space-time metric in a generic coordinate system 

ds 2 = -g tt (r)dt 2 + g rr (r)dr 2 + g xx (r)dx i dx i , (5.28) 

while the bulk gauge fields are A t (r), A x \ = —Et + h(r). We choose a static gauge with 
a = r and t = t and make the ansatz 

x l =X = vt + £(r), (5.29) 

which is motivated by the expectation that the motion of the string will make it have a 
profile that is dragging on one side as it lowers inside the bulk space. 
We may now calculate the components of the induced metric 

9tt = ~9tt + 9xxX 2 = -g tt + g X xV 2 , g T a = 9xxX'X = g xx v£,' , (5.30) 

9aa = 9rr + 9xxX 12 = g„ + g X x£,' 2 , - det g = g tt g r r + 9xx{gtt£,' 2 - 9rrV 2 ) , (5.31) 

and 

-T rji 9<JCF -cr m 9t<T i r or) \ 

ir t =Tf-=gtt, Tv t =-T f -=g tt , (5.32) 

Tt T x =T f (g <Tff X-g tTT X')^=, < = T f (g TT X' - g aT X)-^= . (5.33) 

As all 7f are time independent, the only non-trivial bulk equations are d^it® = 0. It implies 

m 9tt9xxk , , /-- .\ 

ir x = —If — z = constant . (5.34) 

The boundary equation for fj, = t and \x = x are equivalent and become 

T f -^g tt + Ev = -)• tt x = E. (5.35) 



9 
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We can solve ( 5.34| ) to obtain 



y/g- - " '• :> 



r = w— — : 9tt 9xxV * x . (5.36) 

V 9tt9xx ^Tjg tt g xx - ^ 

To ensure we have a real solution, there must be a turning point at r = r s where the 
numerator vanishes and the denominator should also vanish at the same point giving 

2 = 9tt{r s ) 

QxxVs) 

Finally, equation ( |5.35 ) becomes 



Kx = -TfV9tt(r s )g X x(r s ) ■ (5.37) 



Tfy/gtt(r,)gxz(r,) = -E, (5.38) 

which gives the electrostatic force on the endpoint. As v is constant, the equation of motion 

in the x-direction is 

% = -Tt x + qE, (5.39) 

at 

and the steady state solution is tx x = E. 

We denote by J 1 the density of charge as computed from the asymptotic form of At. 



From Q5.37 ), tt x is an implicit function of v. J\ is essentially the blackness function, and 



therefore vanishes at the horizon r = r^. Therefore, the velocity vanishes when r s — > r^. 
For small velocities we obtain 

7t x m -T f g xx (r h ) v + 0(v 2 ) . (5.40) 

Writing J x = J l v we obtain 

jx = jt v ~ Jt ^ ~ J \ E , (5.41) 

Tfg xx (r h ) Tfg xx [r h ) 

and we obtain the DC conductivity and related resistivity as 

Jt T f g xx (r h ) 

a ~ m T^ ' P ~ ^~Tt • ( 5 - 42 ) 

T f g xx {r h ) J 1 

The end result depends on temperature only via the dependence of g xx (rh) on temperature. 
Note that the g xx factor should be evaluated in the string metric. As discussed earlier in 

ADD 

9' xx = 9!Le k +, 0<fc<y^, (5.43) 

and k depends on the relation of <f> to the dilaton. 

Therefore the final result for the temperature-dependent DC resistivity is given by 

J 1 
for massive charge careers and is in agreement with ( |5.14| ). 
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In the case that k = (the scalar field does not involve the dilaton) we can directly 
relate the temperature dependence of the resistivity to the entropy. The reason is that 
from general principles the entropy is proportional to the area of the horizon 



(p-i) 



S~9Zx(r h ) 2 



Relation (|5.45| ) implies 



P(T) 



constant . 



(5.45) 



(5.46) 



In particular for planar systems this indicates a resistivity linear in the entropy. This is a 
very interesting correlation that is emerging from holography. 



6. The 7# = 1 solutions 

For this relation of parameters the general solution is obtained integrating the equations 
of motion fl3.2|) 9 . This range of parameters is most interesting since it encompasses both 
the string theory case 10 , 7 
solution reads, 



5 =- 1, and the Kaluza-Klein case, 7 == \/3, 6 = -4= . The 



d^ 



V(r)dt 2 



<r- \ 
\ r I 



3-<5 2 



4(l-<5 2 ) 
(3-4 2 )(l+« 2 ) 



+ e° 



, dr 2 
V(r] 



+ 



+r' 



V(r) = g)' 



(t) 

} mtf_ 



3-8 2 



2(<S 2 -1) 2 
(3-<5 2 )(l+,5 2 ) 



+ 



dx 2 + dy 2 

(1 + J 2 )^ 2 ' 2 ^ 2 
45 2 (3-5 2 ) 2 r*- 2S2 ' 



(r±) 



3-«5 2 



*2-5 2 



f \ 25 

7. 



m =t 4 / m 2 



V r 

„/?2-<5 2 



(l + £ 2 )q 2 
45 2 (3-<5 2 ) 2 

4,5(<5 2 -l) 



(3-i 2 )(l+« 2 ) 



,4 



$ 



$ 



(3 _ j2) r 3-5 2 

«2-<5 2 



di. 



(3 - 5 2 )r 



3-<5 2 



(6.1a) 

(6.1b) 

(6.1c) 

(6.1d) 
(6.1e) 

(6.1f) 



7,5 = 1, 

where the parameters m and q are integration constants linked to the gravitational mass 



and the electric charge (note that the electric potential (6.1) has been fixed accordingly). 
An initial condition for the scalar, <pQ, can aleays be absorbed in a redefinition of A — > Ae ^° . 

9 We refer the reader to [fed] for particulars on the method of resolution. 
10 This corresponds to 4d string theory with the vector arising from the NS sector. 
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There is an overall scale £ which can be fixed freely as the metric is scale-invariant (up to 
redefinitions of m and q, of course). This scale can be identified as an IR scale via 

3- 5 2 

e = 6 -^- . (6.2) 

-A 
We could also consider "de Sitter"-like solutions by setting 

' 3 " <2 (0.3) 



A 



in the usual fashion, and these solutions can be obtained from fl6,l|) by making the change 
£ —> ia. However, this forces to Wick-rotate both the time coordinate t and the radial 
coordinate r, yielding time-dependent solutions rather than black hole space-times. We 
will thus restrict our analysis to cases where £ 2 is real and positive. It is worth emphasizing 
that what controls the scalar curvature of the solutions, ie., de Sitter or Anti-de Sitter like, 
is the sign of the product (3 — 5 2 )A. We will now focus on the static case, where the radial 
coordinate r is spacelike: 

1. A < and 5 2 < 3, which will yield black hole solutions; 

2. A > and 5 2 > 3, which will yield singular solutions. 



The Ricci scalar, computed with (|3.3j ) is: 



K = 4Ae" 



25 2 A 
3^2' 



m 



3-<5 2 



1- ^ 



,3-<5 2 



3-5 2 (T- 



1 + 5 2 



3-<5 2 



(6.4) 



and it is straightforward using ( |6.1| ) to observe that for 5 2 < 3, both r = and r_ are 
singular. This confirms that in the generic case, the only event horizon is at r = r+, and 
space-time only extends up to r = r_. Quite interestingly, when <5 = 1, r_ ceases to be 
a singular point and turns into an inner horizon, and space-time extends all the way to 
r = 0. 

The solution has an extremal limit when 



m. 



o2 _ S 2 _ VT+Pqi 2 - 



„3-5 2 



25(3 - 5 2 ) 



(6.5) 



The two points r^ then merge and become singular, except again in the limit 5 = 1 where 
the extremal black hole is regular. 

From Q6.4D , we observe that in the case where 5 2 > 3 (which is perfectly admissible 
from previous considerations, provided we take a positive A), the Ricci scalar blows up 
when r — > +oo and vanishes when r — ¥ 0. This is a hint that our space- like coordinate is 
ill-chosen, and that we should change r — > -. However, upon doing this, we find that the 
new radial coordinate is actually spacelike inside the horizon radius, so in the end we obtain 
a space-time describing a cosmological horizon covering a naked singularity. Although this 
is reminiscent of de Sitter space, the inner singularity spoils the regularity of space-time. 
A more appropriate picture would rather be that this depicts the interior of the black hole. 
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The 'string' solution j = 6 = 1 can be shown to be the limit of the equivalent solution 
in |58|, |36| with a horizon of non-zero curvature k = ±1 by taking the limit k — > 0, 5 — > 1 
while keeping the ratio of the two fixed. The question raised in this previous work [58] 
that the limit 5 — > 1 was ill-defined in the string case is then resolved here, a regular limit 
exists but the horizon's topology is severely constrained and becomes toroidal. 

For 5 taken in the range 5 2 < 3, we have Anti-de Sitter like solutions. The solution 
( |6.1| ) is then truly a black hole solution with an event horizon situated at r+ and a curvature 
singularity at r_ . The space-time is static outside the horizon and the singularity is timelike. 
The black hole space-time has to be cut off at r_ and does not extend all the way to r = 0. 
The background space-time obtained either by taking m = q = or the asymptotic limit 
r — > +oo of ( |6.1| ) goes all the way to r = 0, where a curvature singularity sits. This does 
not come as a surprise, as the background solution still contains a non-trivial scalar field 
which is singular asymptotically and thus cannot be treated as a perturbation. It is also 
worth pointing out that the two-torii collapse at r = r_, that is for finite "radius". This 
kind of feature seems to be characteristic of general dilatonic solutions, see [53] for the flat 
potential case in string theory 11 or references |)5|, |66J for related work. 

Finally, the limit 8 — > allows to recover the planar Schwarzschild-Anti-de Sitter 
solution. Indeed, from (qJ), it is seen that the only consistent way to take the 5 — > 
limit is also to set q — > 0, which effectively makes the Maxwell field vanish. Therefore, the 
Anti-de Sitter- Reissner- Nordstrom solution does not seem to be a limit of fl6.1|) . 

The temperature is given by 
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(6.6) 



Let us note already that the behaviour of the temperature in the extremal limit is not 

uniform in the 5 2 < 3 range. Indeed, for <5 2 < 1 -\ — 1=, the temperature vanishes in the 

v3 

extremal limit, whereas it diverges for 5 2 > I+-75 and is finite and non-zero for 5 2 = 1 + -^=. 
Thus, the extremal black hole seems ill-defined in the range 1 + —7= < 5 2 < 3, as it does 
not seem that they can be end states of the evaporation through thermal radiation of the 
black hole. This behaviour has been noted for a long time in string theory black holes and 
is commented upon in |)7], |6q] . From similar cases of uncharged black holes investigated 
in pQ ], this behaviour corresponds to the small black hole solution that is unstable, and 
behaves like the flat space Schwarzschild black holes. In such cases we expect that the 
AdS asymptotics will generate at least one branch, the large black hole branch, that will 
dominate eventually the ensemble at larger temperatures. 

In order to determine the gravitational energy of the solution, one can integrate the 



contribution at infinity of the boundary term ( B.ll ), 

2<5 6 - 15<5 4 + 2(W 2 



M„ 



8vr 



\3-6 2 



+ 



(3-8 2 )(l + 5 2 ) 



-(r_) 



3-5 2 



(6.7) 



11 This solution, or more precisely its generalisation to free 7, was proven in | p2J , |3q ] to be the general 
solution of the system of equations of motion without a scalar potential. 
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where lo 2 is defined as the volume of the compact two-dimensional horizon. This reduces 
in particular for 5 2 = 0, 1 to M g = ^tu. However, the boundary term at infinity stemming 
from the scalar field ( B.12| ) gives a non-zero contribution, 
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which must be added to find the total energy, 
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The conserved electric charge ( B.19| ) is 
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We will make the following rescalings in order to absorb volume factors in the thermo- 
dynamic quantities, throughout the rest of this paper: 
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E 
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Q 

which will not modify the expression of the first law. 

6.1 Regions of validity 

The UV region of the solutions is r — > 00 where the scale factor increases without bound. 
For the solutions to be IR regions of asymptotically AdS solutions the scalar potential must 
become small in the UV region, 
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(6.12) 



Therefore these solutions can indeed be IR regions of asymptotically AdS solutions. 

In the case where the conserved charge is originating on flavor branes with a DBI 
action, according to ( p,16j ) we can linearize the action when 
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As 5 < y3 this approximation is always good in the UV. In the IR we must distinguish 
the following cases 
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For 5 2 < 1 the right hand side of ( |6.13| ) is a function that becomes smaller and smaller 
in the IR. A necessary condition is that at the horizon it is small and this can be 
adjusted by adjusting the UV (boundary) value of the scalar (j) or effectively the IR 
scale £. The condition of validity of ( |6.13[ ) at the horizon is 
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Yh as a function of r + is a function that diverges at extremality. As it also diverges 
at the boundary, it is enough to require 

r _ \ 4+2fc<5 a 

-) »!' (° 5 ) 

for the Maxwell approximation to be everywhere valid in this case. 

For 3 > 5 2 > 1, the factor Yh vanishes at extremality and therefore the Maxwell 
approximation breaks down there. In this case the charge carriers are in a dense 
regime. To test whether in this regime the DBI action can be treated like a probe, 
we must have that V > e 2k ^Z or e(- 5+ "s +2k )^ < 1. Indeed, in this case as 5 2 > 1 
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We can therefore use the probe techniques to calculate the conductivities with the result 



derived already in section 5A 

6.2 Thermodynamics 

Throughout this section and both for the grand-canonical and canonical case, several cases 
will be distinguished depending on the value taken by 5 2 . However, we will only plot one 
curve per range, representative of the behaviour. Where useful, we will also display the 
behaviour at the limiting values for these ranges. 

6.2.1 Grand-canonical ensemble 

The grand-canonical ensemble is defined by keeping the temperature and the electric chem- 
ical potential fixed: this means that we can take as the thermal background the space-time 
with the black hole switched off (m = q = 0), which also coincides with the asymptotic 
limit of the solution if 5 2 < 3: 

ds 2 o = r 2 (-dt 2 + dx 2 + dy 2 ) + (j) dr 2 , (6.17a) 

A = $dt, (6.17b) 

^ - Gf ■ 

For a flat potential 5 = 0, one recognizes AdS^ as expected. 



(6.17c) 
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Figure 7: Slices of the temperature and the electric potential versus the horizon radius r+. 
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Figure 8: Three-dimensional equation of state r+(T, $) for the lower (0 < 5 2 < 1), intermediate (1 < S 2 < 



1 + -4=) and upper (1 + -4= < 6 <3) range from left to rig 
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The equation of state ( pJ3| 
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(6.18) 



gives the horizon radius in terms of the grand-canonical thermodynamic variables (T, $). 
From the previous expression, it becomes apparent why in the string limit the Gibbs poten- 
tial is identically zero for all values of the thermodynamic variables: this is a signal that the 
grand-canonical ensemble breaks down in this limit, as the temperature and the chemical 
potential cannot be considered as independent variables. Thus, in the 5 2 — > 1 limit, the 
grand-canonical ensemble is ill-defined and we shall have to restrict the discussion to the 
canonical ensemble. 

For all temperatures, there is a maximum value to the chemical potential, correspond- 
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ing to 

and the region <£ > <3? e is forbidden because of the non-trivial exponent in ( |6.18| ) , contrary 
to what happens for AdS-Reissner-Nordstrom black holes ||. The extremal black holes 
r_ = r + = r e are attained only for this particular value of the chemical potential, and it 



is straightforward to observe that the equation of state ( 6.18 ) becomes ill-defined in this 
limit, as one cannot know both the temperature and the black hole radius. Thus, one has 
to conclude that the extremal black holes do not exist in the grand-canonical ensemble. 
However, for values of <£ > 3> e , the neutral black holes can still exist and compete with 
the dilatonic background. Indeed, there can be a non-zero chemical potential for a neutral 
black hole while keeping a zero charge, as only the derivative of the vector potential has 
physical meaning. In this case, the equation of state ( |6.18 ) reduces to 

r\- 52 =^^, (6-20) 

and gives the horizon radius as a function of the temperature only. 

We plot slices of the equation of state (16.18; ) at fixed chemical potential and temper- 
ature in Fig.0, while the full three-dimensional plots are in Fig.g. Three ranges can be 
distinguished, 

• Lower range < 5 2 < 1 : The black holes resemble toroidal AdS-Reissner-Nordstrom 
and behave according to the standard physical intuition. Indeed, for finite $, the 
black hole (charged or neutral) disappears (zero horizon radius) as the temperature 
goes to zero, and grows to cover the whole of space-time as the temperature grows 
to infinity. This is to be compared with Fig. 2 in Q: the term responsible for the 
vanishing of the radius at non-zero finite temperature is absent because the spatial 
curvature of the horizon is zero. The extremal limit <& = <3? e can never be attained 
on the left plot of Fig|7] since $ ^ <£ e there, so it comes as no surprise that the zero 
temperature limit switches off the black hole, instead of going to the extremal limit. 
The right plot of Fig.|7] is at finite temperature and so does not display extremal black 
holes. Instead, the extremal limit $ —> $~ leads to a black hole the engulfs the full 
space-time as r + — > +oo. The horizontal axis $ = shows the uncharged dilatonic 
black holes. 

• Intermediate range 1 < 5 2 < l-\ — \=\ This range seems quite different from the usual 

V3 

behaviour as the black hole radius decreases with the temperature! The background 
space-time, given by r + = is attained for infinite temperature and in the zero 
temperature limit the black hole covers the whole of space-time. The right plot of 
Fig.|7| shows a finite radius for the uncharged dilatonic black holes and no black hole 
for the value $ = <3? e , all of this expected as the temperature is finite and does not 
allow extremal black holes to be reached. 

• Upper range 1 -\ — %= < 5 2 < 3 : The same analysis of the left plot of Fig. [7] as in the 
intermediate range applies, while the right plot corresponds to the same analysis as 
in the lower range. 
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Figure 9: Slices at constant chemical potential and constant temperature of the Gibbs potential. The 
T = 2 value for 5 2 — 0.5 allows to see graphically the non-zero limit as $ — ► 





Figure 10: Three-dimensional representation of the Gibbs potential for the lower (0 < 5 2 < 1), interme- 
diate (1 < 5 2 < 1 + -4=) and upper (1 + -^ < 8 2 < 3) range from left to right. 



Then, calculating the value of the Euclidean action (jBjJ) and subtracting the back- 
ground contribution, we find 
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where we have identified the temperatures of the black hole and the thermal background 
on the outer boundary in order to do the subtraction. The Gibbs potential of the black 



hole in the thermal background (6.17) is 
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and here we have to be careful when evaluating the Gibbs potential for extremal black holes. 



Inspecting the limit $ — > $ e (e.g. approached from below) in (6.22), we find that it yields 



52 



zero in the intermediate and upper range, but in the lower is minus infinity: in this range 
we cannot define the extremal limit from ( |6.22 ) as we have an undetermined expression. 
Thus, the proper, unambiguous way to calculate the Gibbs potential of the extremal black 
holes is by evaluating ( |6,21 ), from which we find it is always identically zero. We show 
slices at constant chemical potential and temperature of the Gibbs potential in Fig]9|, while 
the full three-dimensional representation is in FigJlO|. 

We also note that the Gibbs potential is always zero if 5 2 = 1, that is in the string 
case. The explanation is simple: in this limit, the equation of state ( |6.18| ) becomes a 
relation between T and $ and thus they cannot be taken as independent variables; the 
grand-canonical ensemble is ill-defined in this case. 

As it was done in |2C] in the S 2 = 1 case a finer parametrization of the potential 
asymptotics is necessary, V ~ e ^(jP m order to resolve the IR behavior of this regime. We 
will not pursue this here. 



It is also easy to observe, either from ( |6.21| ) or ( 6.22] ), that the neutral black holes 
always have greater Gibbs potential than the charged black holes, and so will not be 
favoured globally as long as the charged black holes exist. 

We may now calculate the entropy (|B.8|), 
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which is the quarter of the area of the horizon as expected, the electric charge 

Q = (3-*V-2g r 3-* 
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which is equal to its usual value ( |6.10 ), and the energy 
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which coincides with the previous expression (|6.9| ) 

Turning to the thermodynamic analysis, we find that the same three ranges as above 
must be distinguished. 

• Lower range: Similarly to the AdS-Reissner-Nordstrom black holes, the dilatonic 
black holes dominate at all values of the temperature, since their Gibbs potential is 
negative, see Fig.J9| The black holes are stable, as can be assessed both from the 
positivity of the heat capacity flB.26|) and negativity of the second derivative of the 



Gibbs potential at constant chemical potential, see Fig. 11. They are also stable to 



electric fluctuations: the slope of the Gibbs potential, as the chemical potential varies, 
is always negative and decreasing, which corresponds to positive electric permittivity. 
However, the charged black holes dominate the ensemble only for <I> < <l> e , whereas it 
is the neutral black holes which dominate for greater values of the chemical potential. 

The remarkable feature that at zero temperature and for large enough <3?, the extremal 
black holes were dominating || is no longer true, as the extremal black holes are ill- 
defined except if <£ = <3? e . As the temperature approaches zero, the black hole shrinks 
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Figure 11: Heat capacity and second derivative of the Gibbs potential at constant chemical potential. 



until it disappears completely and only the dilatonic background is left. This is not 
a proper phase transition as there are no competing solutions: only the background 
exists. If both T = and <]? = <l> e , then the dilatonic background and the extremal 
black holes are competing, but none of them dominate the other. 

» Intermediate and upper range: The black holes, charged or uncharged, are always 
globally and locally unstable and decay to the dilatonic background at all values of 
the temperature and chemical potential. 

The heat capacity at constant chemical potential is given below, 

9. / T \ T=«* / <I> 2 \ M i-a 



$2 

1 ~~ $2 
A e 



(6.26) 



1 -5 2 \?>-5 2 
and is positive for 5 2 < 1 and negative otherwise. 

6.2.2 Canonical ensemble 

The canonical ensemble is defined by keeping the temperature and the electric charge 
fixed. Then (|6.6j ) can be considered as an implicit equation of state giving r + in terms of T 
and Q: 
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We plot it for various values of the coupling constant 5 in Fig, 12 and Fig.[l^. We notice 
three different types of behaviour: 

• In the lower range, 5 2 < 1, there is a single black hole branch for each doublet (T, Q). 
This is similar to the 5 = AdS-Reissner-Nordstrom black holes for q < q CT it, see 
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Figure 12: Temperature and charge versus horizon radius r+. 




Figure 13: Three-dimensional equation of state r+(T,Q) for the lower (0 < (5 2 < 1) and intermediate 
(1 < S 2 < 1 + -2=) range from left to right. 



Fig. 3 in ||. The extremal limit r_ — > r + has zero temperature, and the radius of 
the black hole grows with the temperature. At Q = 0, the endpoint of the curve on 



the right in Fig. 12 is the neutral black hole. For the limiting case 5 = 1, there is 
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a maximal temperature for large black holes, prefiguring the behaviour in the next 
range. 

In the intermediate range, 1 < 5 2 < 1 + -7=, there are two branches, small black holes 

V3 

(SBH) and large black holes (LBH), which merge at a transition point, 
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This radius corresponds either to a maximal temperature at fixed electric charge, or 



to a maximum charge at fixed temperature, see Fig.[12j. The critical point exists only 
for 1 < 5 2 < 1 + —7s, and the critical temperature is given by: 
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Again, this relates to Fig. 3 in ||, except that the third branch is gone due to the 
modified UV asymptotics of the solutions. The small black holes branch end on the 
other side to the extremal black holes at T = and r + = r e , while the radius of 
the large black holes diverges in the zero temperature limit. This alone would cast 
doubt on their physical relevance, see the energetic analysis below. At Q = 0, the 



endpoint of the curve on the right in Fig. 12 is either the neutral black hole (non-zero 
radius) or the background extremal black hole, which in this case has zero radius and 
so coincides with the dilatonic background. 

In the upper range 1 -\ — 7= < S 2 < 3, there is again a single branch. The temperature 

V3 

diverges in the extremal limit r_ —tr+, which consitutes a lower bound for the black 
hole size. This is however expected for dilatonic black holes |67, 68] and signals the 
breakdown of their statistical description. On the other side, the large black holes 
have arbitrarily small temperature. For the limiting case 5 2 = 1 + -7=, the horizon 
size is actually independent of the electric charge and therefore has no minimum : 
the temperature diverges for zero radius. 



Calculating the value of the action ( |B.17| ) in the canonical ensemble (where we have 
subtracted the value of the Euclidean action for the extremal black hole, which we use as 
the thermal background Q), 
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yields the Helmholtz potential of the black hole, 
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which we plot at fixed charge density or at fixed temperature in Fig.|l4| for the two ranges 
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Figure 14: Slices of the Helmholtz potential at fixed charge and temperature. 



The energy of the solution is, flB,23| ) , 
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the chemical potential is, ( |B,25| ), 



3-<5 2 



4(2<5 4 - h5 2 + 1)Q 2 
(5 2( 3 _ (5 2)2 r ,3-^ 



8Q 



5 2 



1 + 5 2 



16Q£' 



2-<5 2 



P 



( 3 _J 2 ) r 3-<5 2 Vl + <5 2 ' 

the entropy (B.24) is equal to one quarter of the area of the horizon, 
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since the entropy of the extremal background is zero, and goes to zero at extremality, 



indicating the presence of a non-degenerate ground state, 67, 68]. We derive the scaling 
of the entropy at very low temperatures and find, 
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We observe that for small temperatures, the entropy scales linearly with the temperature 
for the special value 5 2 = 1 + |v5- The validity of this approximation is for T <C Tq 12 . 



2 the numerical factor in front of the subleading term is of order one for 8 2 — 1 + lyo- 
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Figure 15: Heat capacity and double derivative of the Helmholtz potential against the temperature as a 
function of the temperature. 



Moreover, the entropy will vanish at extremality (zero temperature), for any value of 
5 2 belonging to the lower and intermediate ranges where the black holes are well-behaved. 
This indicates the presence of a non-degenerate ground state and is of direct interest for 
condensed matter systems. We can compare with the neutral black hole entropy, which 
scales like the temperature as, 
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The neutral black holes are well-behaved for 5 2 < 1, so it makes sense only to compare 
them with the charged black holes in this range. The difference between the two leading 
powers in the temperature is, 
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and this shows that the neutral black holes contribution is always subdominant in the ranges 
of interest. This implies that the dominant part of the entropy at very low temperatures 
is due to the charge density. As we will see also in the sequel this seems to be a generic 
property of the theories we are analyzing. 



We also derive the heat capacity at constant electric charge, Fig.|15, 
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which scales linearly with the temperature at first-order for the same special value as the 
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Figure 16: Electric permittivity and double derivative of the Helmholtz potential against the charge. 



entropy, 5 2 = 1 + t\/5, and for T <^Tq 
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The linear approximation has the same validity as for the entropy, when T « Tq. 
The electric permittivity at constant temperature is, Fig,[lq, 
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This ensures that the following first law is verified, 

d{E - E e ) = TdS + ($ - $ e )dQ ^^ d£ w = TdS w + $ w dQ , 
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as expected since we use the extremal black hole as a thermal background. 

Let us study the energetic competition between the black holes and the thermal back- 
ground, e.g. the extremal limit. 

• Lower range: the single black hole branch is energetically favoured over the thermal 
background as the Helmholtz potential is always negative, all the way to the Q = 
axis for fixed temperature, where the neutral black holes are seen to dominate instead 
of the extremal background. Moreover, the black holes are stable both to thermal 



Fig. 15 and electric fluctuations Fig.|lq. 



13 the numerical factor in front of the subleading term is of order one for 5 2 — 1 + %y/E. 
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(see text) 



non-extremal 
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Figure 17: On the left, phase diagram in (T, $) space for the lower range, 8 2 — 0.5. CBH = Charged 
Black Holes, NBH = Neutral Black Holes, DB = Dilatonic Background, EBH = Extremal 
Black Holes. On the right, for comparison, we reproduce Fig.12 of [n] displaying the phase 
space for spherical AdS-RN black holes 



Intermediate range : the Helmholtz potential is always more negative for the small 
black holes branch than for the large black holes branch. The latter crosses the F = 
plane at some (T, Q) but this is irrelevant to the thermodynamics. Both branches 
exist only in a certain region of the (T, Q) plane limited by the line (Tph, Qph) and 
beyond which only the thermal background exists. 

At Q = and for finite temperature, the neutral black holes have positive free energy 
(they are the endpoint of the large black holes branch) and so it is the extremal 
background that dominates (which is the endpoint of the small black holes branch). 
We know from the analyses of p(J that a large black-hole appears in that case above 
a certain temperature, T m i n . For lower temperatures it is the extremal background 
(with discrete spectrum and a mass gap that dominates). At a higher temperature 
T c > T m i n there is a first order phase transition to the large black-hole phase. The 
transition becomes continuous if T m i n = T c . 

The small black holes always dominate the phase space, and are stable to thermal 
fluctuations Fig. 15 but not to electric fluctuations Fig.|i~6|. The large black holes 
display the opposite behaviour. 



• Upper range: Since the Helmholtz potential is positive for all values (T,Q), the 
regular black holes are always disfavoured compared to the thermal background. 
The black-holes are unstable in this range, and the situation is very 

In the case of the lower and intermediate range, these plots are to be compared with Fig. 5 
in H and Figs. 4, 6 in Q in which the so-called branch 3 does not exist as argued above. 

6.3 The phase structure 

6.3.1 The Grand-Canonical ensemble 

We present the phase diagram for the lower range in Fig.[l7], which is the only one that is 
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non-trivial. 

In the region <I> < <I> e , the charged black holes dominate at finite temperature, and so 
do the neutral black holes on the vertical axis $ = 0. The transition to the neutral black 
holes is smooth, in particular no phase transition happens, as is obvious from ( |6,22| ), 

On the horizontal axis T = 0, only the dilatonic background exists (no black holes) 
as the extremal black holes can only be reached at the point <3? = <3? e . The approach to 
zero temperature exhibits critical behaviour, there is an n -order phase transition to the 
dilatonic background for values 

^i<5 2 <^|, n = 4,5,6... (6.42) 

n — 2 n — \ 

Note that the upper bound goes to 1 as n — > +oo, allowing for abitrarily high-order phase 
transitions. 

If <3? — > $~ and T is finite, then the Gibbs potential diverges and there is a zeroth- 
order phase transition to the neutral black holes. For <J> > $ e and finite temperature, the 
neutral black holes dominate the ensemble over the dilatonic background. It makes sense 



on physical grounds to expect that the naked singularity of ( 6.17 ) should preferably be 
cloaked by an event horizon. At zero temperature and $ = $ e , the equation of state ( 6.1§| ) 
for 5 2 < 1 can still be satisfied, but fl6.21|) yields G = identically, and so we cannot 
discriminate at this particular point between the dilatonic background and the extremal 
limit. 

We can compare with the phase space in the case of AdS-RN spherical black holes, 
see Fig. 6 of ||. The effect of missing the asymptotically AdS region (and therefore the 
large black hole branch) is to destroy the Hawking-Page transition, and so the black holes 
dominate the whole phase space. Then, the effect of the scalar field is that charged black 
holes exist only up to some critical value for the chemical potential, and then there is a 
transition to the neutral black holes. Furthermore, the extremal black holes now exist 
only for a specific value of the chemical potential, and do not dominate the phase space 



anywhere. This can be seen comparing the equation of state ( 6.18 ) with equation (20) of 



H , which we reproduce below (notwithstanding some numerical factors) , 



$ = yjrl - 2r + T + n , (6.43) 

where we have explicitly introduced the normalised spatial curvature of the horizon, k = 
0, ±1 for planar, spherical or hyperbolic geometry. Specifically, one sees that for zero 
temperature zero k, one can always obtain in the planar RN case an extremal black hole, 
whereas ( |6.18 ) forces either $ = <I> e (fixed value) to keep r + = r e finite, or r + = (no 



black hole). 

6.3.2 The Canonical Ensemble 



The order of phase transitions is as follows, inspecting the Helmholtz potential Fig. 14 and 
eq. ( pip . 



Lower range: There is a continuous phase transition to the extremal background 
at finite charge in the zero-temperature limit. Inspecting ( |6.31| ), one finds that the 
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Figure 18: (T, Q) phase diagram for the lower (left) and the intermediate (right) range. EBH = Extremal 
Black Holes, CBH = Charged Black Holes, SBH = Small Black Holes. 
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Figure 19: (T, Q) phase diagram in the case of spherical AdS-RN black holes on the bottom line, taken 
from Iql. 



transition is third-order for < 5 2 < 1 %= and second-order for 1 7= < <5 2 < 1. 

There is no discontinuity in any derivative of the free energy at the points 5 2 = 
(0, 1 ± ~7g)> but subleading corrections to the scalar potential might induce critical 
behaviour here. 

For zero charge and small temperature, the neutral black holes dominate. This is the 
range in which (for p = 3) these black holes have a reasonable thermodynamic limit 
even without AdS-completion, see sec. || for details. The phase transition from the 
charged black holes is at least of fourth-order, and of n -order for 



n — 4 ,. 9 n 
r < 8 < — 



n 



n 



1 



n 



4,5,6. 



(6.44) 



Note that the upper bound goes to 1 as n — > +00, thus spanning the whole lower 
range and admitting phase transitions of arbitrary order. Again, there is no transition 
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at the boundary of these intervals, but transitions can be induced through subleading 
terms in the scalar potential. 

Finally, approaching the zero charge axis of the phase diagram at finite temperature 



never shows critical behaviour, since the free energy ( 6.31 ) is a series expansion 
in integer powers of Q. The charged black holes settle continuously and without 
transition in their stable endpoints, the neutral black holes. 

Intermediate range: There is a zeroth-order (discontinuous) phase transition to the 
thermal background at the point (7ph,Qph) since the Helmholtz potential is discon- 
tinuous there and jumps to zero . 

At finite charge, in the zero-temperature limit, the Helmholtz potential and its first 
derivative are continuous, but higher-derivatives diverges for the small black holes 
branch: There is a continuous phase transition to the thermal background of re- 
order for the parameter values 



In — 2 r, In— 1 , 

1+2 Vsrri <s <1 + 2 \fcr (6 ' 45) 

Note that as n — > +oo, one reaches the endpoint of the Intermediate Range, 5 2 = 1 + 

-7=, so phase transitions of any order are possible. At the endpoints of these intervals 
V3 

there are no phase transitions with our choice of potential, but again subleading 
corrections to the scalar potential might change this. 

At zero charge and finite temperature, the (small) neutral black holes are unstable 
and the thermal background dominates. Above some finite temperature a first order 
transition exists to a large black-hole phase. 

At finite temperature, in the zero charge limit, one follows the stable branch (the 
small black-holes) to its endpoint, which has zero free energy and coincides with the 
thermal background. This is consistent with the fact that in the neutral case, the 
black holes are stable only for 5 2 < 1. There is no phase transition because the free 
energy is an integer-power expansion in Q. 

• Upper range: The non extremal (small black holes) are always unstable, so the ex- 
tremal background dominates. There will be large black holes in the higher energy 
completion of the theory, that will dominate above some finite temperature. 

We have plotted the phase diagram (T, Q) for the lower and intermediate ranges in 



Fig.|18|, and reproduced the phase diagram for spherical RNAdS black holes from Fig. 6 in 
Q. As branch 3 is absent in our case, there is no first-order phase transition with branch 
1 in the lower and intermediate ranges : the Hawking-Page transition is absent. In the 
latter, this also explains the zeroth-order phase transition, which in the vocabulary of || 
corresponds to the maximal temperature where branch 1 (small black holes in our case) 



14 Once one corrects to full asymptotically AdS solutions then this will be generically becoming a first 
order phase transition to a RN-like large black hole phase. With special tuning this transition may become 



a second- or higher-order transition as in fed, WH, |46j 
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and branch 2 (large black holes in our case) merge. However, the fact that the extremal 
black hole background dominates on the T = axis is of course unchanged. 

6.4 The AC conductivity 

In this section we calculate the optical (AC) conductivity for small frequencies at zero 



temperature, i.e. for the extremal case of (6.1). We proceed along the lines of section 



The understanding of the small frequency, small temperature behaviour is particularly 
interesting, since the theory of quantum critical behaviour predicts universal scaling laws 
of the form 

M(t(lj) ~ u n , (6.46) 

in the strange metallic region around a quantum critical point. 
The extremal limit of the solution (B.ll) is reached for 



= J = JW^W) =q - (6 - 47) 

Since the low-frequency behaviour of the AC conductivity is only sensitive to the near- 
horizon geometry, it suffices to analyse the fluctuation problem Q4.2J ) in the extremal near 
horizon geometry of ( J6.1| ) , which reads 

l+4<5 2 -<5 4 2S 2 

n (3-S' 2 )(l+S 2 ) f^rrl+S' 2 n 36 4 -4,5 2 -3 

d, 2 = -e p( J 2 ( U p + ? |^ W P 2 ^^d^ + 

qi+s 2 ^ i 

2<5 2 2(a 2 -l) 2 

+fqT+& pW^W+^i (dx 2 + dy 2 ) , (6.48a) 

e* = qT+F pW^W+s 2 ) . (6.48b) 

Here the near-horizon coordinate 

/r\3-5 2 
P={- i ) -Q, (6-49) 

has been introduced, and the time coordinate t = it has been rescaled in order to work 
with dimensionless coordinates only. Note that the AC conductivity is calculated in the 
canonical ensemble, in which the temperature and the charge density of the solution are 
fixed. Following the thermodynamic discussion of (|0]), we may restrict our attention 
for now to the case 5 2 < 1 -\ — y=, which admits a sensible extremal limit in which the 
temperature of the black hole horizon tends to zero. 15 The relevant coordinate change to 



transform the fluctuation problem (4J3) into the Schrodinger problem (4/7) then is given 

by 



26 



■1 



[l + 8 z g 1+a 



:W^ 



<P) = ^1^2_ l P { ' i - S ' Z){1+SA) ■ ( 6 - 5 °) 



15 Note that the restriction coming from the consideration of boundary conditions at the naked singularity 
(see section 2.3) of the uncharged background, 6 < |, is even stronger than this bound. From the point of 



view of the construction of "natural" effective holographic theories the parameter range S 2 > | is thus ruled 
out both by the breakdown of the thermodynamic description for S 2 > 1 H — % and by the "bad" nature of 
the IR singularity at S 2 > |. 
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In this coordinate system the extremal near-horizon metric reads 



ds 2 = D(z)(-dt 2 + dz 2 ) + C(z)(dx 2 + dy 2 



D(z) 



e2~~ 



C(z) = fq^z 



26 



2*4 



J?4 3<5 4 - 6<5 2 - 1 ' 

1+52 z 

1 + 5 2 



> 1+4.5 
'3<5 4 -65 2 -l 



_^43<5 4 -6<5 2 -l 
q T +^ : z 



1 + 5 2 



M 4 -6<5 2 -l 



' _^l35 4 -65 2 -l ' 
g "^ 1 + fi ' 



45(g^ — 1) 

3<5 4 -6,5 2 -l 



(6.51) 
(6.52) 

(6.53) 

(6.54) 



The effective Schrodinger potential (4.7) in the near-horizon region (the horizon is at zh 
—00 for 5 2 < 1 + -t=) then reads 



V = -ff 



2(£ 2 (5 2 -4)-l)(£ 2 (2£ 2 -5)-l) 
(3<5 4 - 6<5 2 - l) 2 



(6.55) 



Note in particular that the extremal near-horizon potential does not depend on the charge 
density q. Since the Schrodinger potential has the same form as observed in [23], we 



would like to apply the matching procedure developed there in order to obtain the scaling 
behaviour of the AC conductivity. To justify this procedure, the full Schrodinger prob- 
lem needs to be frequency-dominated near the horizon and near the boundary in an AdS 
completion, and potential-dominated in between. 

Fig.|2^ shows the full Schrodinger potential for various values of 5 2 < 1 + - 2 =, both at 
extremality (blue solid curves) and slightly away from extremality (red dashed curves). We 
find that away from extremality the potential changes qualitatively only near the horizon, 
according to the results in sec. ||, and in particular that the near-horizon region is separated 
from the asymptotic region by a potential wall. 16 Applying the results of ]^3|, we can thus 
immediately deduce the small frequency scaling of the optical conductivity to be 



<t(u) 



UJ ' 



n 



V4c + 1 - 1 



(3-5 2 )(55 2 + l) 
\35 4 - 65 2 - II 



1. 



(6.56) 



The exponent n{5 2 ) is plotted in Fig.21 for < 5 2 < 3. It diverges at 5 2 = 1 + -t=, and 

otherwise is positive in the region < 5 2 < 1 -\ — %=. At the point 6 2 = 5+ | 33 where 
n vanishes (which is also the point at which the coefficient c in the potential vanishes), 
subleading corrections to the Schrodinger potential will become important. At 5 2 = 1/3 
the minimum n = 5/3 is attained. 

We therefore conclude that in the region where the IR black-holes are stable 5 2 < 
1 + -7=, the solution ( |Q| ) reproduces low- frequency zero temperature scaling laws for the 



16 In the case when the full Schrodinger potential seems to diverge at large r one needs to keep in mind 
that in an AdS completed solution the potential will drop to zero for p — 3 at the boundary according to 

(£§. 
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Figure 20: The Schrodinger potential V(r) for the A x fluctuation around the 7$ = 1 solution at extremality 
(blue solid curves) and slightly away from extremality (red dashed curves) in the range S 2 < 
1 + 2/V3. The non-extremality is parametrized by the non-extremality parameter e via m — 

26(^-6'^) "t" V- F rom l e ft to right, the plots show the respective potentials for (S, e) = (0.7, 0.3), 

(1.3,0.1) and (i/l + 2/y / 3 — 0.03,0.1). The horizon is at the left of each plot at the point 
where the the potential hits zero. The boundary is at the right of each plot. However the 
behavior of the potential there depends on the UV asymptotics not included here, but which 
were described in section 2.4. 



AC conductivity with positive scaling exponents n only. Furthermore, the scaling exponent 
is bounded from below, 



5 

n > — 

~ 3 



.57) 



The spectrum in this class is continuous without a mass gap and the system remains 
a conductor at arbitrarily low temperatures. 




Figure 21: The scaling exponent n of the optical conductivity (j(uj) ~ uj n as a function of 5 2 (0 < 5 2 < 3). 



On the other hand, in the region 3 > 5 2 > l + 2/\/3, the situation is different, as can be 
seen from Fig.^. In this region of 5, the horizon is at z = 0+, and the asymptotically AdS 
region is also at z = 0+, i.e. the z-coordinate is multi- valued. There are two qualitatively 
different regimes, depending on whether the near-horizon potential at extremality ( |6.55| ) is 
diverging to +00 (c > 0), or to —00 (c < 0). In both cases V ~ -% with z — >• and c still 
given by ( 6.55| ). 

The division line between these two regimes is 



5 + V33, 
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at which c(5*) = 0. The left plot in Fig. 22 shows the situation for c > 0: At extremality 
the potential diverges as +l/z 2 . Taking into account the behavior at the UV, detailed 



in section 2.4 we conclude that when the dimension of 



in the UV is \ < A < 1 then 

this describes a finite density system whose charged excitations have a gap and a discrete 

spectrum. 17 It is therefore an insulator. It resembles very much to a 2+1-dimensional Mott 

insulator, where insulating behavior is expected due to strong interactions. 18 

As the temperature is raised from zero, the IR black-hole we found is not thermo- 

dynamically dominant. In analogy with the chargeless case analyzed in a UV complete 

context in [BO] , we expect that the extremal solution will dominate the thermodynamics up 

to a temperature T c > 0. At some lower temperature T m j n < T c a new black hole solution 

will appear due to the behavior of the potential at higher energies. There will generically 

be a first-order phase transition from the gapped solution to the large black-hole solution 

that will be a conductor. This will be an example of a insulator/conductor transition. As 

5 2 — > 1 H — 7= this transition is expected to become continuous. 
v3 

On the other hand for c < (right plot in Fig. 22), the coefficient c > —j as long as 
5 2 < 3. Although the potential diverges, the quantum mechanical problem is well defined 
and we expect not to have any negative energy levels. Therefore in this case the system 
seems to have again a continuous spectrum and is probably gapless. 




Vlrl 
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Figure 22: The Schrodinger potential V(r) for the A x fluctuation around the -yS — 1 solution at ex- 
tremality (blue solid curves) and slightly away from extremality (red dashed curves) in the 
range 8 2 > 1 + 2/y/H. The non-extremality is parametrized by the non-extremality parame- 

From left to right, the plots show the respective potentials for 



ter e via m 



'1+6 



* v ^ ; " I 
2S(3~S' 2 ) "r 



(8, e) = (Jl + 2/\/3, 0.1) and (\\/h+ V^ + 0.002, O.Ol) . The horizon and asymptopia is as 
in figure M, with the difference that the non-extremal case here is without thermodynamic 
significance. 



6.5 The DC conductivity 

We calculate the DC resistivity, for yS = 1 solution using the drag force formula derived 

17 The spectrum is continuous and gapless if A > 1. 

18 John McGreevy informed us that they also found a qualitatively similar solution in a theory with 5 — 0. 
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in section (7.2). The DC resistivity is given in ( 5.44| ). For jS = 1 case, it reads 



/ ( r + 



P^jtl- 



2+25fc 






r 2 



3-<5 2 



2(i' ! -l)(<5' ! -l+2fcd) 
(3-« 2 )(l+« a ) 



(6.58) 



where r e is the horizon radius in the extremal limit and we use the relation r + r_ = r 2 . 
For the canonical ensemble, the total charge density is fixed and thus J* is constant. It 
can be readily checked that the resistivity has dimension [p] = 2 — p using the convention 
of |J27| , because [Tf] = 2, [,P] = p and all the others are in dimensionless combinations. It 
is therefore dimensionless in the planar case. 

To investigate a temperature-dependent resistivity, we use the temperature expression 
(|6.6j) in a slightly different form using r e 



4tt£ 



I ) 



l-<5 2 



1- -4r 



3-S 2 



1-2 



(3-i 2 )(l+<5^) 



(6.59) 



Putting these two equations together we generate three plots in Fig.|2^ for three different 

ranges of 5, namely, lower range with < 5 2 < 1, intermediate range with 1 < 5 2 < 1 -\ — \= 

V3 

and upper range with 1 + -7= < 5 2 < 3. 

v3 



«=0.5, q=l/5000 

5=1.3, q=0.021 
4=1.52, q=0.251 




Figure 23: DC resistivity for the j5 — 1 case. Plots are for the three different ranges of the parameter 
8. The lower range < 5 2 < 1 is denoted with the black solid line, the intermediate range 
1 < 5 2 < 1 H — % with the blue dotted line and the upper range 1 H — % < S 2 < 3 with the 
red dashed line. All the plots display dimensionless quantities, the ir-axis being IT, and the 
y-axis being ^-p with k — 0. The lower part of the 1 < 8 2 < 1 -\ — 1= line corresponds to the 
thermodynamically dominant small black hole. 



Resistivity reveals different behaviors for different ranges of 5, while inside a given range 
the qualitative features are the same. For example, it increases as we increase temperature 
for < 5 2 < 1. It flattens out for decreasing k, decreasing £, but does not change much 
when we vary q. In the intermediate range, it has two different branches corresponding to 
the two black hole phases. The resistivity of one branch increases and that of the other 
decreases as the temperature increases. 

A special case, 5 2 = 1, reveals a limiting behavior of the lower range case, while another 
special case, 5 2 = 1 + -t=, is a limiting case of the upper range case. Those two cases are 
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special because the temperature expression simplifies, and we can convert the temperature 
expression to obtain r + in terms of temperature T. Thus, for those cases we obtain the 
resistivity in a closed form 



T f /r e \2+2fc 



W=l~ T t\- 



(1 - 27r£T) _i ^ 



(6.60) 



while for 5 2 = 1 + -7= we obtain 
V3 



T f 






-V3-\ / 3+2^k 



£T 



f 



V3 



2y/3lT 

\/3-l 



4(^3-1) 



£T 



l+V / 3+2v / 3fc 



(6.61) 
For 5 2 = 1 and k = 1, there exists a temperature where the resistivity diverges, T = -^i- 
This suggests some sort of critical behaviour that deserves further attention. We have also 
checked that the result for 5 2 = 1 is compatible with the one calculated for 7 = 5 solution 
in the next section. 

We will now consider the small temperature limit of the resistivity for the lower and 
intermediate range of 5. This limit corresponds to the extremal limit r + — > r e . Therefore, 
the leading behavior of the resistivity is 



Pleading 



~. 11 (H) 



,p 
J* 



2+2<5fc 



4tt 



-3-<5 2 \r. 

! )+ 
i+6S 2 -:iS 4 



1-8 2 2(S 2 -l)(S 2 -l+2kS) 

{IT) 



l+6i 2 -3i 4 



(IT) 



2(S 2 -l)(S 2 -l+2k5) 



l+6<5^-35 4 



(6.62) 
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Figure 24: Exponent of temperature in resistivity, p ama uT, for two extreme cases, k — 0, 1. There are 
always two distinct points for the exponent to be 1, which reveals the linear temperature 
dependence of resistivity. 



To leading order at low temperatures, there are two different values of 5, one in the 
lower and the other in the intermediate range, for which the resistivity has a linear temper- 
ature dependence. This is the case for all < k < 1. This is seen in Fig.p4|. In p articular, 
for k = 1, the two values are 5 ~ 0.17782 and 5 ~ 1.2359, while 5 = J\ (5 - 2\/5) ~ 0.3249 

and 5 = a/| (5 + 2y5) ~ 1.3763 for k = 0. In Fig. 25, we compared the exact result and the 
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leading order approximation with a set of values, /c = 0,<5 = -iA(5 + 2\/5) , corresponding 
to the linear temperature dependence of the resistivity. 



P 



p = i+2/ys 

Exact 
Leading 



Figure 25: The exact and leading order resistivity in the intermediate range with particular values S 2 
1 + -T-, which gives the linear temperature dependence to leading order. 



In order to assess the range of validity of the leading approximation we will compare 
the leading term to the subleading one by varying r + as r + = r e + 5r with 5r given by 



5r 



2(3 - 5 2 ) 



1+66' 2 -3<5 4 



4-7T 



5 2 



IT 



i+6« 2 -;w 4 



(6.63) 



We obtain 



Pnext— to— leading 
Pleading 



c 



i+as' 2 -3S i 



An 
3^2 



(:T 



\+&S- L -iS' i 



(6.64) 



where the constant factor is C 



-l+206{6+k)-14:5 :i k+5 i (-13+28(8+k)) 



2(-3+<5 2 r(l+«5 2 ) 
We require Pnext-to-leading < Pleading, which is equivalent to 



£T< 



1-8 2 



I, 



(6.65) 



Therefore the leading approximation is valid when the dimensionless temperature scale is 
small compared to an appropriate power of the dimensionless charge density. In particular 
for S < 1, the linear range can be substantial at larger charge densities while for 5 > 1, for 
small charge densities. 
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7. The 7 = 5 solutions 

Taking 7 = 5, one finds the following scaling solution 19 : 

dr 2 



ds 2 
V(r) 

A 



■V(r)dr + e° 



V(r) 



+ r 2 (dx 2 + dy 2 



2mt j2 /- 1 + 



4(l + 5 2 )r 2 



25 



$ 



(1 + 5 2 )r l + 62 

q e s2 



dt, 



7 = 5. 



(7.1a) 
(7.1b) 
(7.1c) 

(7.1d) 

(7.1e) 



We observe that this solution reduces when 5 = 1 to the solution ( |6.lQ for the same value 
of 5. Let us also note that for 5 = (and so 7 = 0, too), we recover the planar version of 
the Anti-de Sitter-Reissner-Nordstrom black hole, as expected from the action. 

The IR scale £ 2 = _a is the same as in the previous solution, and, for the same 
reasons, we shall restrict our attention to black hole solutions, that is A < and 5 2 < 3. 



The Ricci scalar on-shell is: 

K = 4Ae" 



1 



25 2 eV{r) 
(3-5 2 ) 2 r 2 



(7.2) 



so there is a curvature singularity at r = for 5 2 < 3, radial infinity is regular if 5 2 < 3 and 
there are event horizons at the zeros of the potential (|7.1| ) . The size of the torus vanishes 
now at r = 0. 

We now analyze the conditions for V(r) to vanish. Indeed, solving simultaneously for 
V(r) = and ^ = 0, we find the extremal value, 



2n2 



q 2 £ 



2r: 



3-5 2 



4(3 - 5 2 ) 



m, 



(1 + 5 2 )e- 



2-S 2 ' 



(7.3) 



This implies that for r + > r e , or equivalently m > m e , there are two event horizons, one 
inner and another outer, which are degenerate in the extremal case (where the bound is 
saturated). Below the bound, a naked singularity exists. The extremal black hole is defined 
only for 5 2 < 3. 

The temperature of the solution is, 



(3-5 2 ) 5 2_ 2 X _ S 2 



47T 



-1° 



I - I — 



(7.4) 



19 We refer to [B6[ for details on the method of resolution. In the language of this paper, the solution ( |7.l| ) 
is the general solution of the equations of motion provided we set the integration constant h = 0. 
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and vanishes in the extremal case. The integrated gravitational mass of the solution from 
flRll) is 



M g = ^-m, (7.5) 

and in this case the scalar contribution ( |B.12| ) is zero. This is expected since the scalar 
field always has its background value, it is not backreacted on by the black hole, contrarily 
to (|6.1|). The electric charge (B.19) is 



Q = ^-Q- (7-6) 

l07T 

Thus, there are two independent integration constants specifying the solution, as well 
as an independent overall scale that we have fixed to its 'natural' value, e.g. the maximal 
number allowed by the equations of motion. They are m, q and £ as before and there is 
no relation between them, m and q can be considered as the 'reduced' mass and charge, 
while I is the IR radius. However, it is possible that this is not the most general solution 
for 7 = 5, because of the method by which we obtained ( |7.1| ). 

For the remainder of this section, we shall use the same rescalings as for the previous 
7<5 = 1 solution, as in ( |6.12| ). 

7.1 Regions of validity 

The UV region of the solutions is r — > oo where the scale factor increases without bound. 
For the solutions to be IR regions of asymptotically AdS solutions the scalar potential must 
become small in the UV region, 

V~ e - S + ~ (~) — ¥ 0. (7.7) 

Therefore these solutions are indeed IR regions of asymptotically AdS solutions. 

In the case where the conserved charge is originating on flavor branes with a DBI 
action, according to ( |C,16| ) we can linearize the action when 

ZCe k<t > = £ 2 (-J > |g| . (7.8) 

As 5 < v3 this approximation is always good in the UV and it always breaks down in the 
IR at extremality, (r — > 0). 

To test whether in the IR extremal regime the DBI action can be treated like a probe, 
we must have that V 3> e 2k< ^Z or &- °+ k )<t> <g 1 . Indeed, in this case since k > 0, 

e^=n 25{5+k) ^0. (7.9) 

\£J r^O y ' 

We can therefore use the probe techniques to calculate the conductivities in that case with 



the result derived already in section 5.1 . 

7.2 Thermodynamics 

Throughout this section and both for the grand-canonical and canonical case, several cases 
will be distinguished depending on the value of 5 2 . However, we will only plot one curve 
per range, representative of the behaviour in this range. Where useful, we will also display 
the behaviour at the limiting values for these ranges. 
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7.2.1 Grand-Canonical ensemble 
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Figure 26: Slices of the equation of state T(r+, $) at fixed potential (left) and &(r+,T) at fixed temper- 
ature (right) 
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Figure 27: Equation of state r+(T, <E>)in the lower (0 < S 2 < 1) and upper (1 < S 2 < 3) range from left to 
right. 



The equation of state ([74] 



T=(3-5 2 )£ 52 - 2 r 1 + - s2 



2\2p2-25 2 ^2 



(i + pyp-n'Q 

64(3 - 6 2 )r 2 ~ 2S2 



(7.10) 
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Figure 28: Slices of the Gibbs potential at fixed potential (left) and fixed temperature (right) 



can be rewritten so as to express the horizon radius as a function of the thermodynamic 
variables (T, <&) 



32(3 - S 2 ) 
(1 + <5 2 ) 2 $ 2 



T 



+ 



^2 



(l + 5 2 ) 2 $ 2 

3 - <5 2 ' V (3 - 5 2 ) 2 ' 16(3 - 5 2 ) 



+ 



(7-11) 



As in the previous solution, the grand-canonical ensemble breaks down in the string limit 
(5 2 = 1) since the temperature and the chemical potential cease to be independent variables. 
We plot slices at fixed chemical potential and temperature of the equation of state in Fig. 26 



and the full three-dimensional representations in Fig. 27. 
We have to distinguish between two ranges, 

• Lower range S 2 < 1 : There is a single black hole branch. At fixed potential, the black 
hole grows with the temperature and the endpoint of the curve is the extremal black 
hole. Contrary to the previous solution, the extremal black hole can now be reached 
for any value of the chemical potential. At fixed temperature, the vertical axis $ = 
corresponds to the neutral black holes. 

• Upper range 1 < 5 2 < 3 : There is a single black hole branch. At fixed chemical 
potential, the curve starts at zero temperature at the extremal black hole, and then 
the radius actually decreases for the non-extremal black holes as the temperature 
increases. At fixed temperature, the vertical axis $ = displays the neutral black 
holes, and here also the radius diminishes as the chemical potential grows. 

We now turn to the calculation of the Gibbs thermodynamical potential. Evaluating 
the value of the Euclidean action (|B.l| ) and subtracting the background contribution, we 

find 

4Q 2 ^ 2 



-PG = I- I 



„8 2 -2 



pt- 2 {5 2 - 1) 



rV 2 + 



(1 + 5 2 )r 



1+5 2 



(7.12) 
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Figure 29: Three-dimensional representation of the Gibbs potential in the lower (0 < S 2 < 1) and upper 
(1 < 8 2 < 3) range from left to right. 



where we have identified the temperatures of the black hole and the thermal background 
on the outer boundary in order to do the subtraction. The Gibbs potential of the black 
hole in the thermal background (6.17) is then 



g[t,<!>] = I s2 ~ 2 (5 2 - ly'l- 8 ' 2 



i+ 64 r+ 



(7.13) 



slices of which at fixed chemical potential or fixed temperature are plotted in Fig.Eq, while 



the three-dimensional representation is in Fig.2£. 

From this expression, we derive the entropy (|B.8| ), 



S = r\ , (7.14) 

which is the quarter of the area of the horizon as expected, the charge density (|B.9| ), 

Q= i -^S 1 r^%, (7.15) 



16^ 2 + 
which is equal to its usual value (ffl|) , and the energy (| 

„_ 2 .,_,, , AQ 2 e 



E = 2£ 



rV 2 + 



(l + <5 2 )r 



l+<5 2 



Am = M, 



9 ■ 



(7.16) 



which is equal to the gravitational mass ( |7.5[ ) as expected for a scaling dilaton. The 
quantities calculated by use of flB.7|) , ( |B.8| ), (|B.g| ) satisfy the first law flB.21 ). 



Studying the thermodynamics of the solution, two behaviours have to be distinguished: 
• Lower range: At fixed potential, the charged black holes dominate the ensemble 



(Fig.[28| and Fig. 29) , at all temperatures, even in the extremal limit. At fixed poten- 
tial, the black holes also dominate at all values of the chemical potential, including 
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Figure 30: Heat capacity (left) and second derivative of the Gibbs potential with respect to the tempera- 
ture (right) at fixed potential. 
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Figure 31: Electric permittivity at fixed temperature. 



the neutral black holes at $ = 0. The black holes are stable thermally and electrically, 
see Fig.EOl and Fi K M 



Upper range: The ensemble is dominated by the dilatonic background for all (T, 3>) 
values. 
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Figure 32: Slices of the equation of state T(r+, Q) at fixed charge (left) and Q(r+, T) at fixed temperature 
(right). 




Figure 33: Equation of state r+(T,Q) for the lower (0 < S 2 < 1) and upper (1 < S 2 < 3) range from left 
to right. 



7.2.2 Canonical ensemble 

In the canonical ensemble, the equation of state is 



T = ( 3 _ 52^-21-^ 



2r>2 



1 



U 2 Q 



(3 - <5 2 )H 



(7.17) 



and allows to determine implicitly the horizon radius in term of the temperature and the 
charge density. Slices at fixed charge density or fixed temperature are plotted in Fig. 32, 
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Figure 34: Helmholtz potential versus the temperature at fixed charge (left), or versus charge at fixed 
temperature (right), for 5 2 — 0.5 (solid line), S 2 — 1 (dotted line) and S 2 = 1.5 (dashed line). 



while the three-dimensional version is in Fig.|33[ We notice two behaviours depending on 
the value of 5 2 : 

• In the lower range, S 2 < 1 , there is a single black hole branch for each doublet 
(T,Q), and the limiting case 5 2 = 1 has a maximal temperature at large radius, 
again signalling the change of behaviour in the upper range. 

• In the upper range, 1 < 5 2 < 3 , there are two branches, small black holes and large 
black holes, which merge at radius, 

r ^ = W=rf l (7 - 18) 

This corresponds either to a maximal temperature at fixed charge density, or to a 
maximum charge at fixed temperature, see Fig.^. The critical point exists only for 
5 2 > 1, and the line in the (T, Q) phase space so defined is given by 



■Ph 



u d 



S2 _ 2 3-S 2 {4(3 + 5 2 )t 2 Q 2 Ph 



3 + S 2 \(8 2 - 1)(3 - 6 2 ) 



(7.19) 



We now calculate the value of the Euclidean action ( |B.17| ) , subtracting the appropriate 
extremal black hole background, 



I-L 



PI 6 



(1 - 5 2 )r 



2\-iS 2 



2\«2n2 



4(3 + 5 2 )l 2 Q 



+ 



„3-<5 2 



(7.20) 



(l + 5 2 )rl +s2 (1 + S 2 ) e 
from which we can deduce the Helmholtz potential in terms of the thermodynamic variables, 



W[T,Q] 



„& 2 -2 



2\i)2rfl 



(5 2 - l)r%-° + 



_ S 2 4{3 + 5 2 )£ 2 Q 



(1 + 5 2 )r. 



1+5 2 



4£ 2 Q 
(1 + 5 2 ) \3~lP 



3-6* 
2/^2 \ —4- 



(7.21) 
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Figure 35: These plots show the heat capacity (left) and the second derivative of the Helmholtz potential 
with respect to the temperature (right) versus the temperature, for 8 2 = 0.5 (solid) and S 2 = 1.5 
(dot). 



Slices at fixed charge density or at fixed temperature are plotted in Fig.J34|- 
We can check that the first law is satisfied using (|B.7|) , (|B.8| ) and (|B.8|) : 



E\v = E — E e = 4(m — m e ) 
16Q£ s2 



$ 



w 



$-$ f 



[1 + P) 

A h 



-1-8 2 -1-8 2 

r , — r a 



Sw = S = r, = — — . 



(7.22a) 
(7.22b) 

(7.22c) 



The entropy, similarly to regular AdS-RN black holes, is finite at extremality, suggesting a 
large degenerate ground state. Its low temperature behaviour is as follows, 



S w = rl 



I + 2 l + <« 2 - 9 »© + - 



£2-52 'e 



(7.23) 



valid for T < T 20 . 

Here again, we can compare with the uncharged black holes, for which the entropy 
scales with the temperature as 



Sjv ~ T 1 - 



(7.24) 



which for < 5 2 < 1 is always at least quadratic and so is subdominant both with respect 
to the finite and linear piece in ( 7.23| ), in the range of interest. 
The heat capacity is explicitly: 



Cq 



2 r 2 r 4 _ r 4 



1 - 5 2 ri 



4 ' 
Ph 



(7.25) 



20 the numerical coefficient in front of the subleading term will be of order one in the range of interest. 
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Figure 36: Electric permittivity (left) and second derivative of the Helmholtz potential (right) versus 
charge, for S 2 — 0.5 (solid line) and S 2 = 1.5 (dotted line). 



and is displayed in Fig. 35. It scales linearly with the temperature near extremality, 



C, 



Q 



2r 



e n 



1 + (2<5 2 - 3) 



5T 



+ ■■■ 



(7.26) 



valid for T < T 21 . 

The electric permittivity is 



4,(3+S 2 )Q 2 



e T 



1 + 8 2 1 ^A2 1+ (l-« 2 )(3-«52)r* 



,l+5 2 



J 3-c5 2 ' 



1 + 



4Q 2 



(7.27) 



(3-5 2 )r 



see Fig.|36 

We will now study the energetic competition between the black holes and the thermal 
background, e.g. the extremal limit. 

• Lower range: The single black hole branch is energetically favoured over the thermal 
background. The 5 2 = 1 displays an interesting maximal temperature, even though 
the separation into two branches has not yet occurred. The black holes are stable 



both against thermal and charge fluctuations, as seen in Fig.B5j and Fig. 36 



Upper range: The small black holes are energetically favoured both with respect to 
the background and to the large black holes, up until the critical point where the 
two branches merge and cease to exist. When they are small enough, the large black 
holes become energetically favoured compared to the background, but still have a 
greater free energy than the small black holes. The small black holes are stable 
against thermal fluctuations, see Fig.RBJ, but not against electric fluctuations, see 



Fig.|36|. The large black holes have the opposite properties. 



Hhe numerical coefficient in front of the subleading term will be of order one in the range of interest. 
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7.3 The phase structure 

7.3.1 Grand-Canonical ensemble 
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Figure 37: Phase diagram in the (T, $) phase space for the lower range. 



The phase space is non-trivial only in the lower range 5 2 < 1, and is depicted in Fig,|37J 
In that case, the charged black holes dominate for all (T, $), except when T = (extremal 
black hole) or $ = (neutral black holes). There are no phase transitions to the dilatonic 
background in the interior of the phase diagram, the extremal black holes can exist for any 
values (T, <&), contrary to the case 7<5 = 1. These black holes behave more like AdS-RN 
black holes. The only critical behaviour appearing is when approaching zero temperature 
at $ = 0, or zero chemical potential on the T = axis. In both cases there are phase 
transitions of n -order to the corresponding dominating solution in the parameter range 



n 



<5 2 < 



n 



n 



n 



1 



n 



4,5,6... 



(7.28) 



In particular, these transitions are fourth-order or higher. 

In the upper range, the dilatonic background dominates everywhere and there are no 
phase transitions. 

7.3.2 Canonical ensemble 

The order of phase transitions is as follows, inspecting the free energy eq. ( |7.21| ) 

• Lower range: There is no phase transition to the extremal background at zero tem- 
perature for non-zero charge of the black hole. There appear however diverging mixed 
derivatives of the free energy as one approaches Q = T = simultaneously, which are 
of at least third-order. 

In the zero-charge limit at finite temperature, the neutral black holes dominate and 
there is a continuous phase transition from the charged black holes in this parameter 
range. Finally, higher-order phase transitions appear as T — >• 0+ at exactly zero 
charge, as discussed around eq. (|6.44| ). 
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Figure 38: (T,Q) phase diagram for the lower range (left) and the intermediate range (right). EBH = 
Extremal Black Holes, CBH = Charged Black Holes, NBH=Neutral Black Holes, SBH = Small 
Black Holes. 



• Intermediate range: There is a zeroth-order phase transition to the thermal back- 
ground at the point (Tph,Qph) since the Helmholtz potential is discontinuous there 
and jumps to zero. This is again to be remedied by an AdS-completion. 

At zero-temperature and finite charge, the small black holes again settle continuously 
in their extremal state, as seen by the continuity of the Helmholtz potential and its 
second derivative. 

Aproaching zero charge at finite temperature, the first derivative of the Helmholtz 
potential w.r.t. Q diverges, so there is a first-order phase transition to the thermal 
background (which is the endpoint of the stable small black holes branch). In par- 
ticular, on the vertical axis Q = 0, the neutral black holes do not dominate (they are 
the endpoints of the unstable large black holes branch). 

We have plotted the phase diagram (T, Q) for the lower and intermediate ranges in 



Fig. 38 



7.4 The AC conductivity 

We now turn to the optical (AC) conductivity for small frequencies at zero temperature, 
i.e. for the extremal case of (|7.1[ ). We proceed along the lines of section ||. The extremal 
limit of the solution (7.1) is reached for 



3-<5 2 



ra(l + S 2 

2£~ 



l + 5\. 



4£ 2 (l + 5 2 )' 



(7.29) 



3-5^ 

Since the low-frequency behaviour of the AC conductivity is only sensitive to the near- 
horizon geometry, it suffices to analyse the fluctuation problem (|4.2|) in the extremal near 



horizon geometry of (7.1), which reads 

ds 2 = -2£ 2 (3-5 2 )p 2 dt 2 + £ 2 pf 



Px = const . 



dp 2 



2(3 - 5 2 )p 2 



+ / a pg(dz a + <V 



(7.30a) 
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Here the near-horizon coordinate 

r 4 1 + 5 2 n . . 

P = j ~ P* , P* = 3 _p 9 , (7.31) 

has been introduced, and the time coordinate t i— >■ It has been rescaled in order to work 
with dimensionless coordinates only. Note that the AC conductivity is calculated in the 
canonical ensemble, in which temperature and charge of the solution is fixed. The relevant 
coordinate change to transform the fluctuation problem ( [4.5[ ) into the Schrodinger problem 
( |4.7| ) then is given by 

■m—sc^W- (7 ' 32) 



The effective Schrodinger potential (4.7) in the near-horizon region (the horizon is at yu = 
— oo) then reads 

V = ^ ■ (7.33) 



Since the Schrodinger potential has the same form as observed in [23], we immediately 
deduce the small frequency scaling of the optical conductivity to be 

a(co)~uj n , n = 2. (7.34) 



We thus find that the solution (|7. 1| ) gives rise to a low-frequency scaling law for the AC 
conductivity independent of the parameter 5. A similar case of "universality", i.e. inde- 
pendence of the scaling exponent from the parameters in the action, was observed in [23]. 



As a consistency check we can use the fact that the solution (6.1) coincides with the 7 = 5 



solution ( |7.ip for 6 = 1, and hence the scaling exponent ( |6.56 ) agrees with ( |7.34j ) for this 
particular choice of 5. 

In particular it indicates that in this case the charged carriers are gapless with a 
continuous spectrum everywhere. 

7.5 The DC conductivity 

Similar to the DC conductivity calculation for 7<5 = 1, we consider a DC resistivity for 
7 = 5 solution. The resistivity reads 

„ r ^y'^|(!±) 2+Mt , (7.35) 

which is much simpler than 7 (5 = 1 case and the resistivity is a simple power of r + . 

To consider a temperature-dependent resistivity, we use the temperature expression 
fl7.4j ). Resistivity in the 7 = 5 case has slightly different features compared to the j5 = 1 
case. The upper range in the j5 = 1 case disappears and merges to the intermediate range. 
Thus there are two distinct regimes, which are < 5 2 < 1 and 1 < 5 2 < 3. For given 
ranges of 5, the qualitative behavior of the resistivity does not change as we can observe 
in the plots in Fig.p9|. 

For small temperatures, r + approaches the extremal limit r e = I 170JW; ) • Therefore 



.4(3-^) _ 
we expand r + around the extremal r e , r + = r e + 6r, to obtain 5r as a function of T. Then 
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8=1.3, q= 1/200 



Figure 39: DC resistivity for 7 = 5 case. At the left side is a plot for 7 = <5 = 0.7 in the lower range and 
at the right one for 7 = 5 — 1.3 in the upper range, both with k — 0. Plots are generated for 
dimensionless quantities, :r-axis being £T, y-axis being ^-p. The lower part of the upper range 
line corresponds to the thermodynamically dominant small black hole. 



we substitute 5r(T) in the resistivity expression to obtain the small temperature expression 
for the resistivity, 



5r 



IT 



2 T) 



3-5 2 
The low temperature resistivity, including the leading non-trivial order of T, is 



n ~ T f ( r e\ 2+2Sk (, 1 <^ + 25k) 

p ~ jt\e) V 3-5 2 



r \<5 -1 



f. 



IT + 0{T 2 



(7.36) 



(7.37) 



Thus the first non-trivial temperature dependence is always linear unless we have 5k = — 1, 
which always gives a constant resistivity. Note that there is a non-zero resistivity, /?o, at 
zero temperature. This is different from yS = 1 case. These features are depicted in the 



Fig. 40, 
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Figure 40: Resistivity for both leading order (solid line) and exact case (dotted line), which is in the lower 
branch of the upper range with k — 0. Leading order resistivity grows faster than the exact 
case at first, and then they meet each other near the turning point. Plots share the same axis 
conventions as the previous one. 



We estimate the range of the validity of the leading order approximation by requiring 
the next-leading order contribution to the resistivity to be small compared to the leading 
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order term. Note that this condition automatically makes the leading terms much larger 

than the linear term in the temperature. 

We obtain 

, T P \l-S 2 



£T < -r 



2 



(7.38) 



8. Near-extremal scaling solutions 



A solution can be obtained without any constraints on the coupling constants for general 
values of 7, 5. 



ds< 



A ~fh-s) , o e" v dr 2 2 (7-5) 2 , , 2 , 2N 
-V{r)r~ A ~^~dt 2 + ' IA - 1 A ■ ■ 



-wAV(r] 



+ r «™ (dx + dy ) , 



00 -4<^i 




.4 = 2W — e _ 2 00 [r-2m]dt 



wu 



V(r) = r{r — 2m) . 
wu = 37 2 — ( 



27^ + 4 , u = 7 2 - 7^ + 2 , v = 5 2 - 7(5 - 2 , 



(8.1a) 

(8.1b) 
(8.1c) 
(8.1d) 



For future use, we also give the solution in the so-called domain-wall coordinates (fO|). 



This is achieved with the change of coordinates r = p 4 , t = A\/ — je ^°t, and yields 



ds 2 = p- 



(i-sy 



[dx 2 + dy 2 - f{p)dt 2 ] + j?L 



/(p) 



16 (-A) 



ww- 



e -*Ay-iOr-*) a ( p T _ 2m 



e ^o p -(7-<5) 

8 /uA (7+g) , 

iim V u 



> 4 — 2m 



dt. 



3.2a) 

3.2b) 

3.2c) 

J.2d) 



The solution is well defined only if uv < ( since we take A < to have a potential 
bounded from below). 
The Ricci scalar is 



K = AAe~ s<h r ss — 



1 



2(7 -<5) 5 



W1i' 



r + 4m 



(7 - <5) 5 



wu' 



(8.3) 



There is an event horizon at r^ = p^ = 2m and a curvature singularity at 

r = p = for all 5 2 < 3. Indeed, it is not hard to see that if wu > (which is true in 
particular for all 5 2 < 3), 

wu = (3 - 5 2 ){1 + 7 2 ) + (1 - 7 5) 2 > , V5 2 < 3 , 

then the Ricci scalar diverges at r = and is regular at r^. On the other hand, if wu < 
0, then we need to change variables r — > - in order for the Ricci scalar not to diverge 
asymptotically, but then this changes the nature of the solution and we have a cosmological 
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space-time. This is precisely the Gubser bound for such solutions. The region in the 
(7, 5) plane that gives a cosmological solution violates the Gubser bound and so should 
be considered as unphysical: the conditions for having an acceptable naked singularity at 
extremality are 

wu > , uv < (8.4) 

These respective regions in the upper half (7, 8) plane are plotted in the left part of 



figure Fig. 41, Note the symmetry by a rotation of tt around the origin, reflecting that any 
simultaneous change of sign in 7 and 5 can always be absorbed by a change of sign in </>. 

Here again, m is related to the mass of the solution, and <f>o to the overall scale (and 
thus can be fixed at will), but we lack a second dimensionful integration constant. As such, 
the charge density built through the usual Gaussian integral is finite but universal. This 
hints that the solution is not the generic one. 

We will now proceed and take both 7^ = 1 and 7 = 5 in ( |8.1| ), and compare with both 
full solutions we have in our possession in these two limits, ( |6.1| ) and ([Q|), by going to the 
so-called near- extremal limit for both. 

Let us determine the near-extremal behaviour of the full solution (64.) that is, set 22 



„3-<5 2 



„3-<5 2 



+ £, with £ small compared to r_. We obtain: 

2<5^£l 4.5(<S 2 -1) 

e* = r _ ( 1+d )£ (3-^)(i+<^) 



A 

ds 2 



$ + 



s 2 



1 + 5 2 



4(S 2 -1) 



dt. 



3.5a) 

15b) 



- [£(£ - 2m)] £(3-^)(i+^)d£ 2 + 

2(<5 2 -l) 2 

+f(3-* i )(i+**) (dx 2 + dy 2 ) , 



dC 5 



(3 - <5 2 ) 2 £(£ - 2m) 



+ 



(8.5c) 



where we have set r_ ~ r + at zeroth order and m = r 
seen to be exactly m = ^=^ r+ = r_ 



3-<5 2 



o_ r2 

r_ . The extremal limit is 



4h> 



r e . Then, identifying in (|8.1| ) 

Mr :j 



(1+* 2 ) 



<>-2 



3-5 2 



1. 



the near-extremal limit of the 7$ = 1 solution fl6.1|) and the 7(5 = 1 limit of (|3.1|) agree. 



We now turn to the second full solution for 7 = 5 (7.1), and set this time 



,l+<5 2 



_ + £. Then, being careful to keep all the terms to high-enough order, this yields: 



.,25 



A 



$ + V2rI s2 £ 



dt, 



„2S 2 



ds 2 = -f (f - 2m)dt 2 + 



2(3 - <5 2 )£(£ - 2m) 



+ (dx 2 + dy 2 ) . 



J.6a) 
!.6b) 

16c) 



2 We could as well work with r+, this is just a matter of convenience. 
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This is just a product of A0IS2 X R2, with a constant dilaton but with a non-zero gauge 
field, rh has been redefined: 



1 frt 



rn 



1 



2 \n 

where here again rh = corresponds precisely to the extremal limit r_ = r_| 
then compare with ( |8.lD , and find perfect agreement, provided we set in ( |8.1| ) 

-(3-5 2 )A = r 2 = i. 



r P . We 



(8.7) 



Here, we are able to ascribe a precise interpretation to the generic 7, 5 solution (j§,l|). 
Indeed, it turns out that, both in the limit j5 = 1 and 7 = 5, it coincides exactly with the 
near-extremal limit of the full solution. Therefore we interpret ( |S.1| ) as the near-extremal 
limit of the full generic solution. Yet, such a general solution remains elusive and we keep 
a more detailed search for a later work. 

8.1 Extension to arbitrary dimension 

From the p + 1-dimensional equations of motion, one can find the p + 1-dimensional gen- 
eralisation of the solution (|3,2j), using a scaling Ansatz: 

dr 2 



f(r) 



ds 2 = e 2A [-/(r)dt 2 + dxidx*] + 

e^=(-) 2(p - 1) , e* = e«> 
8(p - l)£ 2 (-A) e"^ 



/(r)' 
r\ (<S-7) 



Ww n 



r\ 1 2( P -i) 



/I* 



4(p - 1) /^ 2 Av (7±g) 



<w 



w; p u V it 
w p u = 2(p - 1) + P7 2 - 2j5 - (p - 2)5 2 
We have also reintroduced the IR scale £. 



r^ 2 (p-i) 



r\ 2iF^TJ 



2m 



2m 



n = 7 — 7(5 + 2 , u 



7<5 



(8.9) 
(8.10) 

(8.11) 
(8.12) 



The extremality Noether charge (3.5) is given by 

e-^°£ 2 (-A) 



Q 



m . 



u 



.13) 



x (p-i)(~f-8V 
the temperature T ~ m WpU while the chemical potential as a function of tempera- 



ture 



<I> 



4(p-l) 2+jfr 



« '■/. 



TTUT 



£ 2 Au \£(-A)e- 9 <t>° 



ip—p(y- 



(8.14) 



There are constraints on the parameters: for the solution to exist and have AdS character 

-u>0 , v<0 , w p >0. (8.15) 

These capture the Gubser criterion for the associated naked singularity at extremality. 



They are portrayed in figure 41 for p = 3 and p = 4. 
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Figure 41: This graph shows the Gubser bounds on the near extremal solution on the whole of the (7, 5) 
plane for p = 3 and p = 4. The blue regions are the allowed regions where the near extremal 
solutions are black-hole like. The white regions are solutions of a cosmological type and there- 



fore fail the Gubser bound. The constraints for a black-hole type solution are given in (S.15) 
The dashed blue line is the 7 = S solutions while the solid black line corresponds to the 7<5 = 1 
solutions. 



8.2 Regions of validity 

The UV region of the solutions is p — > 00 where the scale factor increases without bound. 
We have to be careful with the interpretation of this class of solutions in the UV. Although 
they represent exact solutions of the system, the question as to whether their UV region 
can be successfully modified to different asymptotics is subtle. 

For example the solution ( |8.l|) has a potential that behaves as 



V 



~ e 



P 



5( 7 -«5) 



(8.16) 



Therefore, for £(7 — 5) > 0, the potential diverges in the UV, and implies that there is a 
direct problem of considering such solutions as IR limits of asymptotically AdS solutions 
for example. A hint of what really happens is given by the j5 = 1 solutions, which have 
been shown to asymptote to the extremal scaling solutions in the IR, and to the extremal 
uncharged solutions in the UV, (see ( JS.5D and appendix E.5 ). 



In appendix [E] we give a detailed analysis of the system of equations with the extremal 
uncharged and charged solutions as fixed points. We then study the stability of these 
fixed points and show that the uncharged fixed point solution is always unstable, while 
the extremal charged solution is stable once an initial condition is chosen to avoid a bad 
singularity. Therefore, with an appropriate tuning of UV conditions, there is a solution 
that asymptotes to the uncharged solution ( |2.47| ) in the UV and the charged solution flS.ip 
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2 4-4-2 2 



Figure 42: These figures summarize the reliability constraints on the spin-2 fluctuations for p — 3 and 
p — 4. The blue region depicts the part of the (7, 8) plane which satisfies the Gubser bound 
constraints in (8. IE). The yellow-brown region is determined by 6 > 1 and is a reliable region. 



The second reliable region is b < 1 and a > 3 from (D.8) and is portrayed as the purple region. 
Note that this barely touches the Gubser allowed region for p = 3 and has a small overlap for 
p = 4. The solid black line corresponds to 7(5 = 1 and the dashed blue line to 7 = 8. 



or ( |8.8| )-( pTT2| ) in the IR. This is supported by the stability analysis as well as the explicit 
7<5 = 1 and 7 = 5 solutions 23 . 

Because we have shown in section || that the uncharged solutions have vanishing po- 
tential in the UV, this implies that all charged scaling solutions are UV completable. 

In the case where the conserved charge is originating on flavor branes with a DBI 
action, according to ( |C.16| ) we can linearize the action when 



Ze (p-l)A+k<p ^ l( 7 _5)[_4(fe+7)+(7-tf)(p-l)] > 



5.17) 



When (7 — 5) [— 4(fc + 7) + (7 — 8)(p — 1)] > 0, the Maxwell approximation breaks down 
near the UV. In the opposite case it breaks down in the deep IR. 
Finally the probe approximation is valid 



ZVe 2( P -l)A 



~ p^to-V » q 2 , 



(8.18) 



in the Maxwell limit. The probe approximation thus always breaks down in the IR for 
p > 3. In the strong self-interaction limit this same condition becomes p( s ^^>( s+ ^ +2k > <g 1 



according to section C.2 



23 In reference |69] that appeared after the present paper, this same statement was argued and shown 
numerically to be true in several concrete examples. The discussion here has been added after the appearance 
of |69], and in order to justify the UV completability of the solutions. 
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8.3 Reliability of the IR fluctuation problem at extremality 

The extremal black- hole solutions are generically naked singularities (except when 7 = <5). 
The region in the (7, 5) plane where these satisfy the Gubser bound is portrayed in figure 



41. To asses the IR reliability of such naked solutions, we must investigate the Sturm- 
Liouville problems for the fluctuations. This is done in appendix |D| and we review our 
findings here. 

The spin-2 Sturm-Liouville problem is reliable at extremality when 

b > 1 or (b < 1 and a > 3) , (8.19) 



where a, b were defined in ( |D.2| ) and ( p,4j ). The regions are portrayed in Fig. 42. 

On the other hand, the spin-1 (current) Sturm-Liouville problem is reliable at ex- 
tremality when 

• b > 1 or 



b < 1 and a < — 1 where a is defined in ( D.14| ) or 



b < 1 and 2 > a > — 1 and q 2 > q 2 , where go is defined in (D.18). 



Therefore, as is shown in Fig. 43, there are no constraints on the (7, 5) plane from the 
spin-1 fluctuations (though the third constraint imposes that the charge density should be 
large enough). For p = 3, only the first and the third constraints have some overlap with 
the Gubser constraints, while for p = 4 the second spin-1 constraint has a (quite) small 
overlap with the Gubser ones, too. Overall, the spin-1 problem is always reliable in the 
Gubser region for p = 3, while for p = 4 there is a small part of the Gubser region where 
there is a lower bound on the charge density. 

A further comment concerns the correlation of the constraints above and the behavior 



of black holes near extremality. As discussed in section SJi the near extremal black holes 
are thermodynamically unstable in the region 

b-W-W > 1 (8.20) 

w p u 

where the temperature at extremality diverges. In analogy with the uncharged case, in such 
situations the system is expected to have no stable black hole solutions up to a temperature 
T m i n above which new large black holes will appear due to a modified UV potential. This 



region is described in Fig. 44. 



This region is exactly the same as the b < 1 region above. Therefore, in the thermody- 
namically unstable region, one also needs a > 3, so that the spin-2 problem is reliable. For 
p = 3, this is never really the case as this region touches the allowed region at two points, 
(7, <5) = (1,-1) and (7,5) = (—1,1) where v = 0. We conclude, that in the p = 3 case, 
for the thermodynamically unstable region the spin-2 problem is not reliable. The spin-1 
problem in the same region is on the other hand reliable except the case b < 1,5 < —1, 
which has a a non-zero overlap with the allowed region and where the spin-1 problem is 
reliable only at sufficiently high charge density. 

Finally, in the thermodynamically stable region, the spin-1 and -2 problems are always 
reliable since it corresponds to b > 1. 
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Figure 43: These figures summarize the reliability constraints on the spin-1 fluctuations for p — 3 and 
p — 4. The blue region depicts the part of the (7, 8) plane which satisfies the Gubser bound 
constraints in (8.15). The yellow region is determined by b > 1 and is identical to both spin-2 



and thermal stability constraints. The purple region is determined by b < 1 and — 1 < a < 2 
in (D.17), with a lower bound on charge density q > q in (D.18), while the green region is 
b < 1 and a < — 1. The overlap of the latter region with the Gubser one is trivial in p — 3, 
and small but non-trivial for p — 4. The solid line corresponds to 7<5 = 1 and the dashed blue 
line to 7 = 5. 



8.4 Lifshitz solutions 

A generalization of AdS solutions involves asymmetric scaling symmetries of the Lifshitz 



type, [70|. A generalized Lifshitz metric can be written as 

-odr 2 



ds z 



dt 2 dxMx 4 
+ 









By changing variables to u 



As' 



the metric can be rewritten as 



,dn 2 



5V 



;2z 



+ 



dxMx* 



To bring this metric to the domain wall frame we must identify 

A(r)_}_ ....■■. 1 - ~ du 

U 

so that 



a9{r) 



,1-z 



,2(2-1) ' 



dr 
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10 -((«-«3 
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o9(u) 



{z-lf 



5.21) 



$.22) 



.23) 



(8.24) 



Comparing with the scaling solutions ( |3.2j ) , we find that they are of the Lifshitz form 
only when 5 = and 7 < 0. In that case the non-extremal solutions read 

„2 r ,-.x £±^ii m=§ 



A(r) 



((' 



1 



1 

z-1 



o9( r ) 



1 



f 2 



1 
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r \^^i- 



e <Ao r , 



.25) 
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Figure 44: These figures summarize the region where the (small) near-extremal black holes are thermo- 
dynamically unstable. The blue region depicts the part of the (7, S) plane which satisfies the 
Gubser bound constraints in (BT3). The purple region is the thermodynamically unstable 
region determined by 1 — 2 "~ ' < 0, as indicated in equation (5.31). The overlap region 



is the region where there will be a gap in the black hole temperatures. The black solid line 
corresponds to f8 — 1 and the dashed blue line to 7 = 6. 
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Specializing to the planar case p = 3 the parameters are 
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Such solutions were interpreted as Lifshitz solutions before in [51]. 
The non-extremal version provides explicit Lifshitz black-holes 
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The fluctuation spectra for such solutions pose no constraints. 

We conclude that such solutions are generated by a scalar along a flat direction of a 
potential, which drives the coupling constant of the U(l) to weak coupling. 

8.5 Thermodynamics 

We will discuss here the thermodynamics for p = 3. Other dimensions show similar behav- 
ior, and the equivalent expressions are straightforward to compute. 
The temperature is, 

T = — \/^Ae~^»(2m) 1 - 2il = L , (8.31) 

4-7T 

and we observe that it vanishes in the m = (extremal) limit if and only if the exponent 
is positive, that is 

wu - 2( 7 - 5f = (3 - 5 2 ){l + 7 2 ) + (1 - 7 <5) 2 - 2( 7 - 5) 2 > 0. (8.32) 

This precisely reduces in the limit we already observed in the 7<5 = 1 solution, 5 2 < 1 + -3=. 
We recover the two standard dilatonic cases, 

• If 1 — 2 ^ 7 ~ - > (lower range), the extremal temperature is zero. 

• If 1 — 2 ^ 7 — < (upper range), the temperature diverges in the extremal limit. 
The gravitational mass is 

M g = ^J^-( 7 -5fm, (8.33) 

4-7T V TO 

and reduces to the corresponding near-extremal limits of the 7<5 = 1 and 7 = 5 solutions. 
Notice that if 7 = 6 the mass is identically zero, which is expected since in that case the 
geometry is not a black hole but simply the direct product AdS^ X S^- The electric charge 

U} 2 fuK 7^0 O 



Q = 7T\— e— *°, (8-34) 

8-7T V U 

is universal and does not contain an independent integration constant. Finally, the entropy 

is 2 

5 = ^(2^)2^ (8 . 35) 

and we observe that it vanishes at extremality provided that 7 7^ 5 (otherwise it is finite) 
and 1 — 2 ^ 7 ~ - > 0. This is consistent with the behaviour of the two previous solutions 
for 7# = 1 and 7 = 5. 

Given all of this, we work in the canonical ensemble, where only the temperature is 
allowed to vary (and the charge density is fixed to its value (8.34). Then the Helmholtz 
potential is 

,(7-<5) 2 ' 



1 



m, (8.36) 



^2 / r -ic 

* - " 1 — z^ 

wu 
and again reduces both to ( |6.31| ) and ( 7.21| ) in the appropriate limits. Unsurprisingly, 



8vr 



inspection of its sign reveals that the black holes are globally stable in the lower range but 



93 



Thermal 
Instability 



H 



Thermal 
Stability 




17 



Figure 45: On the left panel, we plot the region of local stability of the near extremal black hole. The right 
panel shows a variety of phase transitions of the near extremal black hole to the background 
at zero temperature. In the blue region continuous transitions occur, in the purple region 
adjacent to the blue one the transitions are of third-order. The stripes starting with yellow to 
the left of the blue and purple regions depicts transitions of fourth- (yellow) up to tenth-order. 
Above them all higher-order transitions also occur. 



not in the upper range, similarly to the 7^ 
capacity, we find 



1 solution. Moreover, computing the heat 



C, 
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^2 (7 - £) 2 



wu 
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2(-y-6) Z 
4:irT \ um-2(7-<5) 2 



.37) 



which will be positive in a range a little bit smaller than the lower range and negative in 
the upper range, thus rendering the black holes locally stable in the former and locally 
unstable in the latter, as displayed in the left panel of Fig.fl5|, 

Finally, second- or higher-order phase transitions will occur in the zero temperature 
limit if the appropriate derivative of the Helmholtz potential ( |S,36| ) diverges. In particular 
an n -order transition occurs if 



1 < 



wu 



--n<0. (8.38) 

wu — 2(7 — by 

This is depicted in the right panel of Fig.^, showing the regions in which transitions of 
second to tenth order occur. Two remarkable features are that, in most of the parameter 
space that is thermodynamically stable (which amounts to requiring a well-defined extremal 
limit and vu < 0), second-order transitions appear, and there are no first-order transitions. 
We note that it does not make sense to analyse the thermodynamically unstable region 
in this way as long as the AdS-completion has not been taken into account, since our 
experience with the solutions for j5 = 1 show us that this is crucial for the correct structure 
of the phase transitions. 
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Figure 46: In the shaded region, the entropy at finite charge density dominates the one for zero density 
at very low temperatures. In the rest of the (7, 5) diagram, the comparison cannot be made 
as both entropies are expected to be of 0(1) instead of 0(N 2 ). The vertical axis represents 
the value of 8, while the horizontal axis the value of 7. 

Finally we would like to compare the entropy at finite charge density to that for zero 
charge density. We have already done this for the 7 = 5 and the j5 = 1 solutions and 
found that in both cases, the finite density entropy is dominant compared to that of zero 
density for low temperatures. For the extremal solutions we are discussing (at finite charge 
density) the entropy behaves at low T as 



2 (7-*r n 

S ~ T n,u-(~f-5)^ ^ 



5.39) 



while at zero density from ( 6.36J ) S ~ T 1 - 6 ' 2 . When 5 2 < 1, both powers in ( |8.39 ) and ( |6.36j ) 
are positive and the entropy at finite charge density dominates the one for zero density. 
When 5 2 > 1, there are two regions. In one, both ( 8.3S| ) and ( |6.36D are negative. In this 
region the theory has a gap and the theory remains in the extremal solution as T — > 0. In 
the rest, the neutral extremal solution dominates at low temperature with a temperature 
of O(l), while at finite density the entropy is non-zero and 0(N 2 ). Therefore in this last 
case again the entropy at finite charge density dominates the one for zero density. This is 



summarized in figure 46 



We stress again that all these considerations should be taken as valid in the near- 
zero temperature regime: as the temperature grows, one will need the (yet unknown) 
full solution and not its near-extremal approximation in order to study finite temperature 
thermodynamics. 

8.6 The AC conductivity 

We now turn to the calculation of the AC conductivity at zero temperature, i.e. in the 
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extremal limit, and low frequencies. We will first work in general dimension p + 1, and then 
analyse the special cases of planar systems (p = 3) and bulk systems (p = 4) . The scaling 



solutions at arbitrary dimension are given in subsection 8.1 



The low-frequency scaling of the AC conductivity for the solutions ( |8.1| ) is again cal- 
culated from the general formulae (4.7), fl4.18| ) and ( 4.1SJ ), wherever these are applicable. 
It is convenient to work with the solution written in domain wall coordinates ( |8.8| )-( 8~T^ ). 
From (8.2) it is evident that the extremal limit is m — > 0, after which the coordinate p (or 
equivalently £ in ( |8.8| )) is already the near- horizon coordinate - the horizon is at p = 0. 
Requiring the temperature of the solution (^3) to go to zero in the extremal limit yields a 
constraint 

(wu) p > (p - 1)( 7 - 5) 2 . 



The relevant coordinate change to arrive at the Schrodinger problem is given by 



dz u(wu) p e 



dp 8A(p - 1 



<50O (p-l)("[-6) 2 -(v]v.)p 

p 



2(p-l) 



?.40) 



5.41) 



Thus if the condition for a extremal zero-temperature limit flB.40 ) is fulfilled, the power in 



.41) is always smaller than minus one, and hence the coordinate change is polynomial. 



The Schrodinger potential found from (4.7) is again 



V(z) = -, 



(p-l)(5 2 - 1 2 -4)(8(l-p) + (5- 7 )( 7 (l+p) + 5(5p-7))) 



4(2(1 - p ) + (5 - 7 )( 7 + 5(2p - 3))) 2 
The AC conductivity at zero temperature scales then for small frequencies as 



a ~ uj 
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.42) 
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As a consistency check we observe that both (|8.42|) and ( 8.43| ) reduce to corresponding 
results for the special cases 7 = 5 (eqs. ( |7.33 ) and ( |7.34|) ) and 7^ = 1 (eqs. ( |6.55 ) and (|6.56| )), 

and reproduce the results of |23|]. For the Lifshitz solutions (5 = 0, 7 = • 



1 ) we find 

2—1 / 



n 



(3(z-l)+p) s 



J.44) 



which reduces to a ~ uJ 1 for p = 3. This is different than the co~ 2 ' z result of [27]. The 
difference can be traced to including the full backreaction as well as the non-trivial scalar 
factor in front of the gauge field action. 

We now turn to a more detailed analysis of the special cases of planar systems, p = 3, 
and bulk systems, p = 4. The above calculation again made use of the matching procedure 
of |2J| which, from experience with the 7<5 = 1 solution (|6.1| ) and the fact that ( |8.1[ ) 
reduces to fltO]) for 7 = 1/5, is applicable in the region of the (7, 5) plane given by ( 8.4CQ , 
i.e. whenever the solution has an extremal limit in which the temperature drops to zero. 
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Figure 47: Contour plot of the scaling exponent n in the (7, 8) upper half plane for p = 3 (left figure 

< 5 < ./§) and p = 4 (right figure, < 8 < ,/§)• Left figure: Contours correspond to 
n — 1.52, . . . , 8.36, starting with n = 1.52 in the upper right corner and increasing in steps of 
0.76. The black solid line 7 = 8 is n — 2, and brighter colors correspond to larger n. Right 
figure: Contours correspond to n — 2.2, . . . , 12.1, starting with n = 2.2 in the lower right 
corner and increasing in steps of 1.1. The black solid line 7 = 8 is again at n = 2. In the 
yellow regions the computation of n cannot be trusted, since an explicit AdS completion of 
the space-time is needed to render the thermodynamics well-defined. The scaling exponent 
diverges to +00 along the dashed black line in both cases. 



Planar systems (p = 3) The left plot in Fig. 47 shows the scaling exponent n in the 
upper half (7, 5 < a/5/3) plane for p = 3. The maximum value of 5 is set by the spin 2 



fluctuation constraint, see sec. 2.3. The yellow shaded region is either forbidden by thermal 
instability (above the black dashed line, see Figj45|), 24 by the spin 1 fluctuation constraint 
Q2.43 ) (the wedge-shaped region touching 7 = 5 = 0), or by requiring that the r-coordinate 



is space- like (yielding u > 0, v < 0, which forbids the yellow region on the left below the 



6(j) 



7 + 2v / 3~+f 



$.45) 



black dashed line). See figure [4l| for a discussion of all the constraints on (7,(5). In the 
allowed region the scaling exponent is always positive, n > 0. We find n = 2 along the 
solid black 7 = 5 curve, with n increasing until it diverges along the dashed black curve 

1 

3 

For comparison, the j5 = 1 line is also depicted as the dotted black curve. We thus 
conclude that whenever the underlying scaling solution is thermodynamically stable and 
completable to AdS space, our system describes a AC conductivity scaling law at zero 
temperature and small frequencies with positive scaling exponent, a ~ uj n>0 . 

In particular this implies that in this regime the spectrum of charged fluctuations is 
gapless and continuous and therefore describes a conductor. 

24 Note that the thermal instability constraint is not to be taken too serious since, as discussed before, an 
AdS completed solution with these parameters will have proper thermodynamics and hence a meaningful 
holographic interpretation. 
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Figure 48: p=3. Left figure: The region on the (-y<5) plane where the IR black holes are unstable, and 
the coefficient c is positive. In this region the extremal finite density system has a mass gap 
and a discrete spectrum of charged excitations, when the UV dimension of <j> is A < f . This 
resembles a Mott insulator and the figure provides the Mott insulator "islands" in the (7, S) 
plane. Right figure: The region on the (71$) plane where the IR black holes are unstable, and 
the coefficient c is negative. In this region the extremal finite density system has a gapless 
continuous spectrum at zero temperature. In both figures the horizontal axis parametrizes 7, 
whereas the vertical axis 8. 



The region for which 



(wu) p < (p-l)(j-Sf 



$.46) 



is the one where the "small" black holes found in the IR EHT are thermodynamically 
unstable and the region where the extremal background is dominant thermodynamically. 
In this case we have two possible behaviors for c defined in ( 8.42 ): 

(a) c > 0. This region is depicted at the left of figure |4S|. Here the potential for charged 
excitations at T=0 diverges to infinity in the IR, and also in the UV if A < 1. These systems 
have a gap and a discrete spectrum of charged excitations at zero temperature and therefore 
resemble Mott insulators. 



(b) c < 0. This region is depicted at the right of figure 48. In this case the system is 
gapless with continuous spectrum at zero temperature. It is therefore a conductor. 
The behavior above generalizes what was found in the j5 = 1 and 7 = 5 solutions. 



Bulk systems (p = 4) This situation is shown in the right plot in Fig.|47[ The conclusions 
are similar to the planar case, in particular the positivity of the scaling exponent in the 
physically reasonable region, and the divergence along the dashed black curve, in this case 
given by 

5( 7 ) = - (27 + v^a/iO + 3 7 2 ) . (8.47) 



The spin 1 fluctuation constraint is now violated in the horn-shaped region, and the r- 
coordinate is timelike in the yellow region on the left below the black dashed line. One 



should also note that for general dimension p, the scaling exponent always behaves as 

n p ( 7 , 7 )=2 (8.48) 

along the line 7 = 5. 

Concerning the low energy spectrum some general properties are different when p > 4. 
The Schrodinger potential always diverges near the AdS boundary, and therefore when it 
also diverges positively in the IR the spectrum is gapped and discrete. This indeed happens 
on "islands" in the (7, 5) plane analogous to those of the p = 3 case. Moreover one again 
expects here the insulating phase to transform via a finite-temperature phase transition to 
a conducting phase. 

All the rest of the allowed regions have a continuous and gapless spectrum at extremal- 
ity. 

8.7 The DC conductivity 

From ( 5.44| ) we evaluate the DC conductivity for massive carriers using 

o ■ ~ a s (on) ~ e 2^(Po) e ^(Po) „ T m p 2k(p-l)(5- 1 ) + 2(5- 1 ) 2 

Pmasswe 9xApQ) e e i ' mp - 2 (p - 1)(1 - 6(6 -j)) + (5 -j) 2 ' 

(8.49) 

where we have extended the calculation to p+1 dimensions. This contribution is inversely 

proportional to the charge density. 

In the planar case p = 3 the relevant exponent reads 



4k(5 - 7) + 2(5 - jf 
4(1 - 5(5 - 7)) + (5 - 7) 5 



<±rvyu - 1 -r ±\u - 1 , c ,m 



The exponent becomes unity for two values of 7 



7± = 35 + 2k ± 2^/1 + (5 + k) 2 . (8.51) 

For a non-dilatonic scalar, k = and the temperature dependence of the entropy and the 
resistivity are the same. As the entropy scales with the same factor as the heat capacity 
we find that linearity of the first implies linearity of the second. 

For the Lifshitz solutions, we must take 5 = and 7 = ~\ ~nEjv- hi this case we 
obtain that 



z 
in agreement with [27] when k = 0. 



! + l ^-«, (8.52) 



We can also estimate the conductivity from pair production from ( [5,13| ), 

. , n ^ 7 -i T n p n (7 ~ S) [2(p - 1)7 + (p - 3) (7 - 5 - 2k(p - 1))] 
PP*~9M>Z ~Tv , n p = 2(p _ 1)(1 _ (5(5 _ 7)) + (< j_ 7)2 

(8.53) 
For the Lifshitz solutions the exponent becomes 



jM^fe-iW-iX,-!) 
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In the planar case p = 3 and 



47(7 - 5) 
U3 = 4(l-^- 7 )) + (5-7) 2 (8 ' 55) 



This exponent will become 2, if k = — f together with (p. 51). Therefore when 



k = -- , 7± = 25 ± V4 + 5 2 (8.56) 

the DC conductivity interpolates between a T behavior for massive carriers and a T 2 
behavior for light carriers. As the linear behavior is inversely proportional to the charge 
density we expect the T 2 behavior to dominate eventually at higher doping as indeed 
happens in cuprates. 

There are two special values in ( S.56| ), for k = 0, 5 = and 7 = ±2. One of them 



corresponds to the z = 2 Lifshitz geometry. The other is k = 1, S = —2 and 7 = — 4± 2V2. 



9. Strong self-interaction solutions 



In section |C.2j we have derived the equations for solutions that describe strong self-interactions. 
We reproduce them here in the domain wall frame 

1 n" \n\p-9+k<t> 

(p-l)A"+ l -r = 0, pA ' + 9 ' + 9 7 - l j^ =w = 0, (9.1) 

(p - l)A'(g' + P A') - \S 2 - e~ 9 V + \ q \ e -9-<p-i)A+W = , (9.2) 

with the gauge field being in this limit 

A' t = e A+k * . (9.3) 

We have found the following three charged solutions of the equations above: 
9.1 Strong coupling solution I 



e A = rW ^j) ; e <P = e <Po r c P ; e g = e 9o r ci 

with 



1-(^ C: 



(9.4) 



c p = 0, ci = 2, c 3 = l,2, e S0 = |e fc ^, q = -2Ae-( 5+fe ^ . (9.5) 

For this solution 

T ~ ro , S = constant . (9-6) 

These are AdS2 black holes driven by strong charge self-interactions. 
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4 

2(p-l) ; 


e <t> = e <Po r c P ; 


e 9 = e ffo r ci ; 


(9.7) 


C\ = 


= 1 2p AS + k)\ 


g = -2Ae-( 5+fc ^° . 


(9.8) 



4 



9.2 Strong coupling solution II 

e A -. 

with 

c p = 2(5 + k), 

p — 1 

This is an extremal solution. 

9.3 Strong coupling solution III 

For 5 = ±1 and p = 4, there is another solution 

e =r4, e ^ = r^4, / = _|1 ( r4 - r 4 ) , | 9 | = -V , A t = -{r 

(9.9) 
We can obtain the temperature and entropy as 

AttT = A^ , S ~ e 3A = r| , (9.10) 

giving 

4 4 /4vrT\t 3 

«=-, * = -,(— ) . S~T,. (9.11) 

Consistency with the strong coupling limit implies that 

q»e^ +k)4,+ZA . (9.12) 

This can be valid in the IR depending on the sign of 7 + k. 
The DC resistivity for massive carriers is 

p~T 1± l fe . (9.13) 

It is linear if the scalar is not the dilaton. 
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Note added 

After the appearance of this paper, |6^] appeared which independently analyzed the scal- 
ing solutions of section ||. Moreover it has connected them with the uncharged extremal 
solutions in the UV, and defined the physics by using the non-conformal holography. Our 
appendix |E|, added later to respond to questions posed by the referee is in agreement with 



|6S] although the techniques used here are somewhat different. We thank E. Perlmutter 



for discussions on his work. 
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APPENDIX 

A. Equations of motion for the EMD system 

The relevant Einstein equations stemming from Ql.l|) in the general coordinate system Q2.5| 



arc 



r + {p-i) 



c" p-i fry , b' , c\ c 

~C 2~ 



D + B + c) C ° ' 



-r7-^ + 2{c-D- U + ^ + (3 ~ p) " 



g 2 £ 



- inr + £zi!£f!^S + 2 



1 

2^ 2 C V 2 C D 

along with the solution for the gauge field 

F rt = A' t = ^y . 

ZC* 



c ) zcp- 1 

q 2 B 



0, 



2ZCP- 



0. 



In the domain- wall coordinate system of (2.1) the equations of motion read 



= (p - l)A" + \r , 

n = £_ 4- „' + r>A' - 

g< ^ 9 ^ P e 2(p-lM+»Z(0)(/ ' 

= (p - 1)AV + pA') - e-^(0) - i 



'2 



+ 



2e 2{p-l)A+g Z l 



(A.l) 

(A.2) 
(A.3) 

(A.4) 

(A.5) 
(A.6) 

(A.7) 



B. Definition of thermodynamic ensembles 

B.l Grand-canonical ensemble 

In this appendix we will identify the various thermodynamical ensembles we will work 
with. To do this, we consider the Euclidean version of action Ql.l| ) in p + 1 dimensions 
and supplement it with appropriate p-dimensional boundary terms, in order to make the 
variational problem well-defined. Although these boundary terms might be neglected in 
order to obtain classical solutions of the equations of motion, they are crucial once one tries 
to calculate the Hamiltonian associated to these solutions |[71J| . We will start by adding 
only the Gibbons-Hawking boundary term W2h to the usual Lorentzian action: 



16tt J m 



R 



1 



(dtf 



1 



e 70 F 2 + 2Ae - 



+ 



1 



-ill 



dM 



'h p K 



(B.l) 



where dAi is the boundary of our space at radial infinity, hfu, the induced metric on the 
p-boundary and K p = hp U V ^riy its extrinsic curvature. n M is a vector orthogonal to the 
boundary and pointing outwards. A convenient way to calculate the boundary term is to 
use the following relationship: 

° ' ^, (B.2) 



dM 



<h p E K p 



dn 



dM 
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see [72]. We then consider the variation of this action flB.l| ) [73]: 



51 e = (Terms giving the equations of motion) [5g^ v , 5A^ , 



+ 



+ (Gravitational boundary term)<% 
ib7r JdM V 

'ITT f \AW)n^- 



(B.3) 



The Gibbons-Hawking boundary term allows to fix only the metric components on the 



boundary and leave their first derivatives free. Variation of (B.l) will yield the equations 
of motion if the variation is at fixed chemical potential on the boundary, which we will 
denote <3?, as well as at zero variation of the scalar field at the boundary (we will see 
however that this variation is not independent of the other quantities varied). 

This will define the grand-canonical ensemble as we use the gauge potential at infinity 
and the temperature as parameters and they are intensive variables. Again, it will not be 
necessary to define a thermodynamic variable related to the scalar field. The other extensive 
variables S and Q are determined as functions of them. The Gibbs thermodynamical 
potential in the grand-canonical ensemble is 



G[T, $} = E-$Q-TS, 



(B.4) 



and is related for a particular solution to the Euclidean continuation of the Lorentzian 
action by Ie — 1% = —f3G where 1° is the background contribution. The first law is 
expressed as: 

dE = <S>dQ + TdS , (B.5) 



or equivalently in terms of the Gibbs potential: 



dG = dE- TdS - SdT - $dQ - Qd$ = -SdT - Qd$ . 



(B.6) 



We can then deduce the other thermodynamic quantities of the solution by taking the 
appropriate derivatives of the Gibbs potential with respect to the thermodynamic variables, 



E 

S 

Q 



G-T 

_dG_ 

~df 
dG 

<93> T 



dG 
df 



$ 



dG 



T 



(B.7) 
(B.8) 
(B.9) 



Specialising to p = 3 (four-dimensional case), we can also determine the various ther- 
modynamic quantities as conserved charges stemming from the symmetries realised asymp- 
totically on the boundary, following the Regge-Teitelboim-Henneaux procedure, 0, 74]. 
We detail this procedure in the case of the conserved charge associated to translations in 
time on the boundary, that is, the energy. Indeed, it has been shown in other contexts 
that a scalar field could backreact strongly on the metric in such a way as to modify the 
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asymptotic boundary conditions, and thus generating extra contributions to the variation 
of the conserved charge associated to the energy, |75], |65|, |66| . 

Although our setup generically breaks the 50(4,2) symmetry of the Anti-de Sitter 
space-time on the asymptotic boundary, see ( |2.7D , considering the timelike Killing vector 
4r generating the time isometry associated to the energy will be enough for our purposes 
and allows to follow the procedure defined in [71, 74, 75, 35, 66 1. Referring to these papers 



for details of the computation, one uses the following minisuperspace metric 

r{r) 

which yields the following boundary terms (the simplified analysis of |7t| is still valid in 
our case, with appropriate boundary conditions) 

SB g = ^[N (RR'5f - (f)'RSR) + 2/ 2 (N5R' - N'SR)}™ , (B.ll) 

with primes denoting derivatives with respect to r, f3 the inverse Euclidean temperature, ui2 
the volume of the compact planar horizon. r + means the inner boundary of the space-time, 
which for a black hole is the outer event horizon. 



To this gravitational contribution, one must add a scalar term, |7^, 65, 66 1 



«J, = j^[^/V**]~. (B.12) 

The contribution on the asymptotic boundary of these two terms will yield the variation 
of the energy, that is 

5E = -SB™ - SBf , (B.13) 

which can then be integrated. One should also take into account the boundary term 
stemming from the integration by parts of the Maxwell Lagrangian, but this does not 
present any particular difficulty. Summing all these boundary terms, both at the horizon 
and at infinity, will give the Gibbs potential. 
The electric energy is 



$Q = -L / ^N |e^F"V4)l , (B.14) 

167T J goo L J 



which allows to identify the conserved electric charge, the Aq term on the boundary being 
the chemical potential at infinity. The chemical potential for a black hole is usually fixed 
such that the Maxwell vector A^ is always regular on the horizon of the black hole. This 
implies fixing $ so that Ao(r+) = 0. 
This yields finally, 

E = M g + M lj> , (B.15) 

G = Mg + M^-QQ-TS, (B.16) 

for the energy and the Gibbs potential. 



105 



B.2 Canonical ensemble 

On the other hand, consider adding to the action the following Maxwell boundary term 

(B.17) 



Ie = Ie + ^t [ Jtf s e*F'»'n li A v , 

16vr J M V 



which will combine with the previous electric two-dimensional boundary term when the 
action is varied and yield: 

8Ie = (Terms giving the equations of motion) [Sg^ , SA^, 6(f)] 
+Gravitational boundary term + 



+- 



(B.18) 



16vr 

1 
!6tt j dM 

Here, the equations of motion will be obtained when we keep the quantity [y—hF^n^ 
fixed on the boundary at infinity, which is precisely the usual electric conserved charge, 



167r JsL 



(B.19) 



This is the canonical ensemble, as thermodynamical variables are the temperature (in- 
tensive) and the electric charge (extensive parameter). The free energy in the canonical 
ensemble is the Helmholtz potential, 

W[T,Q] =E-TS, (B.20) 

related to the Euclidean action by — J3F = I — Iq. The first law is expressed as: 

&E = TdS + $dQ , (B.21) 

or equivalently as 

dW = dE - TdS - SdT = -SdT + $dQ (B.22) 

We can then deduce the other thermodynamic quantities of the solution by taking the 
appropriate derivatives of the Helmholtz potential with respect to the thermodynamic 
variables, 

aw 



E w = W -T 



Sw 



dW 



dT 
dW 



dT 



Q 



Q 



dQ 



(B.23) 
(B.24) 
(B.25) 



We also define various quantities which will be useful to evaluate the thermodynamical 
stability of the solutions 0], namely the heat capacity at constant electric charge and 
dilaton charge, 



Cr 



T 



dS_ 
dT 



(B.26) 
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and the electric permittivity, 



-H'S) t . (B.27) 



Stability of the solutions requires that both quantities are positive. This corresponds 
respectively to the fact that larger black holes should heat up and radiate more, while 
smaller black holes should go colder and radiate less; and that the chemical potential 
should increase when more charge is added to black hole || , making it harder to move away 
from equilibrium, as expected from classical physics or more generally, from Le Chatellier's 
principle. 

Equivalently, one can require that 

d 2 W\ fd 2 W\ 

as can be seen by varying the free energy in the canonical ensemble ( [B.20D . 

Finally, let us note that all these boundary terms in the case of black hole space-times 
should be considered both on the boundary at spatial infinity, the outer boundary of space- 
time, but also on the horizon, the inner boundary. In the case of the Maxwell boundary 
term, this is not a problem: indeed, there is a gauge freedom in choosing the value of the 
chemical potential on the outer boundary, <J>. It cannot be set to zero for a black hole space- 
time as the coordinate transformation required to do that would be singular on the horizon 
1 72], so we fix it so that it is zero on the horizon: this way, the Maxwell boundary term 



never contributes on the inner boundary. On the other hand, for the Gibbons-Hawking 
boundary term, including the inner boundary contribution would precisely mean that we 
do not fix the temperature of our thermodynamical system, but rather its entropy: this 
means going over to the microcanonical ensemble. However, in order to carry calculations 
in this ensemble, we would need to be able to count the available microstates needed to fix 
the entropy to a particular value. Unfortunately, this requires a quantum theory of gravity: 
some progress in the case of extremal (or so-called BPS) black holes have been made in 



some specific realizations of String Theory [76], but this is as far as our knowledge extends 



for the present. Thus, we will leave aside the microcanonical ensemble in the following. 

C. DBI dynamics and the strong interaction limit 

In the case where the U(l) symmetry arises from fundamentals, the gauge field originates 
from a stack of flavor D-branes. In this case the proper gauge field action is the DBI action 
that includes "self-interactions" of charge. 



Sdbi 



f d p+1 x Z(<f)) V- det( 5 CT + F) - v 7 - det g° , (C.l) 



where we have adjusted the action so that the potential is still given by V in ([O]). Therefore 
in this case the full action is S g + Sdbi instead of S g + Sf- The metric that enters the DBI 
action is the induced metric in the string frame. As we do not consider extra embedding 
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coordinates, we can rewrite it using the Einstein metric using ( |l.6|) . 

S DB i = - fd p+1 x e 2k(f, Z^)[^- det{g + e~ k <t>F) - y'-dets] , Z = M p ~ l e^ H 'Z . 

(C.2) 
Here the coupling Z has been defined such as to appear in the kanonically normalised 
Maxwell term, -Z(4>)F 2 /4. 

The equation of motion for the gauge field is 



e + **ZJ-det(g + e'^F) (g + e -^Ff v - (g + e~ k(p F) 



0„ 



while the variation with respect to the metric gives 



k(j)-p\un 



0. 



(C.3) 



1 5S DB i Ze 2 ^ det{g + e" fc fF) 



-g 5g^ u 4 y detg 

, ab . 



g» a (g + e **f)° gb, + (ji++v) 



Ze 2k4> 



where (<? + e fc( ^F) is the inverse matric of (<? + e fc ^F) , . 

For a non-trivial F r t profile in the coordinate system (|2.5|) we obtain 



- det(g + e-^F) = C p -\DB - e'^F 2 ^ , 
and the equation of motion (C.3) becomes 



Or 



Solving we obtain 



ZC^Frt \ n ZC^Frt 



^DB-e- 2k ^F 2 t j 







^DB - e- 2k ^F 2 t 



■rt 



4 



q^/DBt 



yjq 2 + ^fc^C^" 1 ) 



The on-shell action is 



rion she// 
a DBI 



d p+l x e^Z 



det(flf + e~ fc fF)| - ^detg 



K / dre 2fc *ZC— \/M 



r 



e^^ZcV 



\ vVe^ + e^Z 2 ^" 1 



V p / dr 



\f~DB 



(p-i) 



2 ZC^ 



TT + W 



The last line is the result in the linearised (Maxwell) approximation. 
The Einstein equations ( |A.5| )- ([A.7|) now become 



r + (p-i) 



~c 



D' B' , o ,C"\ a n 



D" C" 4-p (C D'\ (Bf D' 
~D~~7T + 2 \~C~U) \~B + ~D 



q 2 e 2k( i'BC- R 2 1 
^/ q 2 e 2k<j> + e 4k4> Z 2 C p-i 



0, 



(C.4) 



(C.5) 



(C.6) 



(C.7) 



(C.8) 



(C.9) 



(CIO) 
(C.ll) 
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(C.12) 
In particular, in the domain wall frame, (D = e +9 , B = e~ 9 , C = e ), the equations 
above translate to 

1 „" 2 -2(p-l)A -g+2k<j> 

(p - 1)A" + -<f>' 2 = , M' + </ + ^- q =0, (C.13) 

2 g g'Je Ak( l>Z 2 + g2 e -2(p-l)4+2fc^ 

(p - 1)A'( 5 ' + pA') - -0' 2 - e~ 9 (y + e 2fc< ^) + e'^e^Z 2 + g 2 e -2( P -i)A+2fc^ = . ( C14 ) 
There is still a conserved Noether charge 

Q = ( J>A' e9 _ A ^ = 0. (C.15) 

dr 

and it vanishes at extremality. 

C.l The weak self- interact ion limit 

In a region where 

Z(0)C i£ 2 iI e^>|g|, (C.16) 



we can linearize the DBI action and the equations ( |C.11| ) and (C.12) reduce to the Maxwell 
equations in ( |A.5[ )-( |A?7| ) as can be directly verified. 

C.2 The strong self-interaction limit 

If the strong self-interaction limit it is the opposite inequality that is satisfied 

Z(^)C^e fc *«|g|. (C.17) 

In the domain wall coordinates, the DBI square root is approximated as 

Z 2 e (p-l)A+3kcj> 



e 4fc0 Z 2 + q 2 e 2k<t>-2{p-l)A ^ I | e k«-(p-lM + e + _ _ _ 

2 M 

and the rest of the equations become to leading order 

(p-l)A" + -^ = 0, M' + <7' + ^-L X = 0, (C.19) 

(p - 1)AV + MO - tU' 2 - e~ 9 V + \ q \ e -3-(P- 1 ') A + k ^ = . (C.20) 

Note that in this limit the gauge coupling constant function Z(<f>) drops out of the 
equations to leading order after taking into account the inequality ( p.lT] ). Moreover the 
gauge field solution becomes 



F rt = A' t = sign(q)e A+k * 
Therefore to leading order the gauge field solution is q- independent. 



Z 2 e 2(p~l)A+2h(t> 

1 v — + 



(C.21) 
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C.3 The probe limit 

In this limit the presence of the gauge field solution is subleading in the gravity plus scalar 
equations, with equivalent condition 

V » e 2k <*>Z ( - e2k * Z ° P ~ _ 1^1 (c.22) 

In the linearized limit this becomes 

V»^T, (C23) 

while in the strong self-interaction limit it is 

V^e 2k ^Z. (C.24) 

D. The Sturm-Liouville problem for spin-1 and 2 in general charged back- 
grounds at extremality. 

We will study here two special fluctuations of the system, namely spin-2 excitations and 
the spin-2 excitations arising from the gauge field at extremality. This will give us infor- 
mation on the reliability of the extremal naked singularity as to the spectral properties 
of the system. If the related Sturm-Liouville problem is well-defined with a normalizable 
and a non-normalizable solution, then the associated physics is largely independent of the 
resolution of the singularity. In the opposite case it depends in a crucial way. 

D.l The spectrum of the Laplacian 

The spectra of the spin-two excitations are readily computable as they satisfy the standard 
Laplace equation. In the background above and for a plane wave basis for the perturbation, 
h(r)e lujt+ it becomes (in Euclidean space) 



J +pA ) h -\P^ + T? 



fc'' + (^ +M ')*'-( " + * U = (D.l) 



At extremality (m = in ( 8.10D ), / ~ r v and 



r +rJt = t , a = 2 + f^l + ^. l + ^_ (D . 2) 

/ r 2 2(jp - 1) 2(p - 1) 

while the terms depending on frequencies or momenta scale as some inverse powers of r, 

Pe" ~ r 2(2 + s il -S)) + ^ ' fe 2A ~ r2+5(7 - 5)+ rf • l ^ 

The Gubser constraint ( |8.15|) implies that a > 1. 

Therefore, near the singularity r = 0, equation ( p.l| ) becomes 

h"+*h'-^- b h = , b = -v+ { ^l > 0. (D.4) 

r r A0 2{p — 1) 
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with w ~ co. Equation can be solved by the modified Bessel functions, and the two 
independent solutions are 



h± = r 2 < 



IV 



l-b 



'M I TTTif 



*-» i nV* 



(D.5) 



1) 



When 6 < 1, the argument of the Bessel function is small near r = and the solutions 
behave as 

hi ~ r 1_a , hi ~ constant (D.6) 

Using the standard norm 

dr r a h(r) 2 (D.7) 



/o 

we observe that hi is normalizable at r = when a < 3 while /12 is always normalizable as 
a > 1. The Sturm-Liouville problem is therefore unacceptable when 1 < a < 3. Therefore, 
when 

6<1 , w p u< (p-l)(7-<5) 2 (D.8) 

the dynamics of the low temperature system is not reliable and a resolution of the singularity 
is needed to access it. The relevant region is shown in figure ^ in the main text. 

When b > 1, the argument of the Bessel function is large near r = and the solutions 
behave as 

b-a ±- 



h± ~ r — e 16-ik"- 1 , (D.9) 

Therefore one is always normalizable and the other is always not normalizable. we conclude 
that when b > 1, the Sturm-Liouville problem for spin two is always well defined. 

There are a few marginal cases that need to be checked separately. In the a = 1 
case the two linearly independent solutions are Iq and Kq. For b < 1, both solutions are 
normalizable and therefore this case is not reliable. For b > 1 only one is normalizable. 
The other case is b = 1. In this case the two solutions behave as 

h± ~ r ^iH (D.10) 

Since h + is always normalizable the constraint that h- is non-normalizable reads, (1 — 
a) 2 + \w 2 > 4. Since w ~ ui in order for this to hold for all u>, me finally must have 

6 = 1 , a>3 (D.ll) 

This is still in agreement with ( p.8| ) with 6 = 1. 

D.2 The spectrum of current excitations 

The spectrum of current excitations at zero momentum k = can be computed in the parity 
odd sector, in which the Ai and gu fluctuations decouple from the rest of the fluctuation 
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spectrum, from eq. ( t4,5|) . In the domain wall coordinate system used in (8.8) eq. (4.5) 
reads in Euclidean signature 



= ^'+(| + y + (p-2M'M! 



u- 



+ 



f 2 e 2A fZe 2 (P~V A 



A; 



A; H A; 



m\ 77*2 



A; 



(D.12) 
(D.13) 



p(j — 5) 2 h — d) 2 . a ph — S) 2 ,„ 

a ^ 2- ^r — +- < 2 , b = ^- — ^ - v > , 6=- + ^rr: — tt = x > ( D - 14 ) 



mi 



2(p - 1) 
8(p-l)(-A^ 2 ) 



2(p - 1) 



2 4(p - 1) 



„2 2 (5-7)00 

2/;26 2 9 W P U e 



8(p-l)(-A) 



(D.15) 



For 6 > 1 the frequency term dominates and, as in the spin-2 case, the Sturm-Liouville 
problem is always acceptable. For b < 1, the charge term dominates as r — > and the 
solutions behave as 



i-a±-y/(i-a) 2 +4m| 
Af ~ r 2 



(D.16) 

Again the A^~ solution is always normalizable. For A~ to be non-normalizable, we must 
have 

^J(l-a) 2 + Am 2 >2 (D.17) 

As rri2 ~ g this is always satisfied for sufficiently large charge density 



2(p - l)e(T- <5 ) < ^ ^ A 



■ 2 >a 2 ^ ^~^: - l-^)(3 P(7 " 5) 



2(p-l 



2n / p(7 _5)a 



2(p - 1) 



(D.18) 



When a < — 1, this is satisfied for any charge density. In the opposite case 2 > a > — 1, 
the constraint ( p. 18 ) is necessary. 

For b = b' = 1, the solutions behave as 



Af ~ r" 



-a±y / (l-a) 2 +4(m2+ m 2) 



(D.19) 



Asking that this leads to an acceptable problem for all frequencies we end up with the 
constraint ( p. 17 ). 

To summarize, there are no constraints on 7, 5 coming from the spin-1 Sturm-Liouville 
problem. However for b < 1 there are constraints on the charge density summarized in 
( |DT7D . 



The various constraints are summarized in figure E3 for p = 3 and figure 44 for p = 4 



E. The gravitational system in phase variables: fixed points and stability. 

In this appendix we will formulate and study a different form of the system of equations 
presented in appendix [A]. They involve three phase variables satisfying a non-linear first- 
order system of equations from which all the background functions can be reconstructed 
by direct integration. The advantage of this formulation is that it is largely coordinate- 
independent and it is useful to study local stability of fixed points and the global profile of 
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the solution space. It is also the most efficient formulation for finding numerical solutions. 



This is a generalization of the approached developed in [17],[20|. A related approach was 
developed [J59]|, but it involves five first-order equations and is coordinate-dependent. This 
work investigated the global profile of the space of solutions and we find results qualitatively 
similar. 



We will define the phase variables first in the domain wall system (2.1, 2.2, iOD and we 
will then show that they can be defined in a coordinate independent fashion. 
In the domain-wall coordinate system, the equations of motion read 

0=(p-l)A" + ±<p a , (E.l) 

„ = ^ + /+ ^-_ 5 _ I ^__, (E.2) 

- (p - 1)AV + P A>) - e->V W - \r + 2g2(p _ 1 ,l, z ■ (E.3) 



along with 



F rt = A' t = — j^—- . (E.4) 



To proceed further we will need the second derivatives 

rh' 2 a 2 

All V II /2 / \l " tT? r\ 

A =-W^T) ' 9 = ~ 9 - P9A+ e 2 ( P -^Z^) > (K5) 

f = -(g' + pAW - e-°V> - 2e2{ fZ +9Z , ( E - 6 ) 

where Z' = %,V'=%. 

We define the three phase variables as 

1 6' d Q 2 -2(p-l)A 

xw = vm^rA m = ^ ' ww= ~^v— (EJ) 

In terms of these we can express A, g as 



1 r d<t> I p rd<t>Y _ . 



These determine A, g as functions of <p. 
We then use the following relations 



(E.9) 



to obtain 



# _ V{W - 1) f 

Tr ~ ±X V p{p-l){X 2 -Y-l) " (E - 10) 



113 



which upon integration gives (f>(r). Therefore knowledge of X ((f)) ,Y ((f)) , W(4>) permits the 
reconstruction of A(r),g(r), <f)(r) from (E.8) and ( E.lOj ) by standard integration. 

Although we defined the phase variables in the domain wall coordinate system, they 
can be defined in the generic coordinate system (|2.5|) as follows 



XI 



y/2p(p - 1) d\og(C) ' 



= 2 dlog(D/C) 
W) p dlog(C) 



W((f>) 



AV *"* 



2BDV 

(E.ll) 

The final ingredient in this approach is the first-order system satisfied by the phase 
variables that is summarized below. 



(4 



^2 



dX 



dW 




P 



20-1) 
dY X 2 



(X' 



Y-D 



p 



T^ + §W 



Y -\ 



y/2 P (p - 1)X 



pY + 



2p X(W-1) 

2(p-l)W 
W -1 



(E.12) 

(E.13) 
(E.14) 



This system has as input functions the potential V and the gauge coupling constant 
Z. For the system studied in this paper, V = Vo e j Z = e ^^ an d the equations above 
become 



dW 



7' 




2(p-l) 1 



P 



X 



w. 



dX 



P 



2(p 
dY 



-(X 2 -Y -1) 



X 2 - Y - 1 
^2p(p - l)X 



pY + 



p — 1 jW — 5 
2p X(W-1) 

2(p-l)W 
W -1 



(E.15) 



(E.16) 



(E.17) 



These equations have three integration constants. The determination of A, g in terms 
of (f> in equation flE.SD has another two. The extra integration constant to determine A(r) 
is a diffeomorphism artifact in this case although it can be traded for an initial value for the 
scalar (f)Q. The physical set of integration constants are three: <pQ,T,Q. The two last ones 
are decided in the system above. The third constant of the first-order system corresponds 
to the scalar vev and should be tuned so as to avoid a (bad) naked singularity. 

Eqs. flE.12j )-( jE.14 ) can be inverted and give 



dlog(W-l) 



v 



+ 



p — i 



X 



Y \ d\og(X 2 
2XJ ( 



Y-l) 



Integrating, we obtain 



W((f)) = 1 



(X 2 -Y-l) f^f* ^x-Z) 



V 



(E.18) 



(E.19) 
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The scalar invariants of the background can be expressed in terms of the phase vari- 



ables. The scalar curvature is given by 

p + 1 U 



R 



-2U + 



V 



P -1X 2 -Y-1 
while the scalar and gauge-field kinetic terms are 



2 " -x* 



p + 1 



2V(W-1)- 



X 2 



ZF Z 



Y -\ 



AWV . 



(E.20) 



(E.21) 



X 2 - Y - 1 

Singularities of the scalar invariants are signaled by the vanishing of X 2 — Y — 1 or the 
divergence of W, V —> oo. The horizon is signaled by a divergence of Y. In particular 
X 2 = Y + 1 is a bad generic singularity. 

The first-order equations above are equivalent to the original equations if (f) is not 
constant. It can be shown that only in the case 7 = 5 a solution <j> = <j>q constant exists, 
and therefore this case must be studied separately. In what follows we assume 7 7^ 5. 



There are special solutions to the system ( E.15 )-( E.17 ) that relate to the vanishing 
of the right-hand sides of equations. These are fixed point solutions to the equations and 
as we will see they coincide with the charged and uncharged extremal solutions as well as 
with completely singular (and unacceptable) solutions. The rest of the solutions interpolate 
between various fixed points solutions. 

E.l Uncharged solutions, W = 

W = is a fixed point of ( E.12J ) and describes uncharged black holes. 



The rest of the equations become 

dX I p~ 



(X 5 



Y-D 



2(p-l) 

dY X 2 - Y - 1 



p-l 6 



2p X 



d<j) ^/2p{p - l)X 
Note that the second equation can be solved as 



pY . 



(E.22) 
(E.23) 



Y{4>) 



I p r4> du 

~\J 2(p-l) J X(u) 



C + 



I P C^ ^ v 1 - x2 ( u ) r ~\l W^=T) ^ 4fJT\ 
Y 2(p-l) J aa X(u) e 



(E.24) 



For an extremal solution we must have Y = 0. 
The fixed point solution here has 



Y = 



X = a 



p — 1 
2p 



6. 



(E.25) 



This corresponds to the extremal solution in (2.47) for tq = 0. We can also deform Y by 

keeping X = a to obtain 

I -a 2 
y = \ a . . (E.26) 



p(l-a 2 )(0-0 o ) 

g V 2 p(p _1 ) a 



1 
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This has a singularity hidden by a horizon. This coincides with the solution in ( 2.47 ) 
In this case we can write the general solution in the implicit form as 



(X-a 



,2a 



{X -i) a + 1 {x + i) a - 1 



O k(4>-(j>o) 



\X\>1, 



and 



(X - a) 2a X + 1 
(l-X 2 ) a X -1 



3 fc(</>-</> ) 



\X\ < 1 



A: 



2(p-l) 



p-_^ 2 
2p 



1 



(E.27) 



(E.28) 



This solution captures all the flows among the three fixed points ±l,a and ±oo. In 
particular, the Gubser bound is satisfied when \a\ < 1. The fixed points at X = ±1 are 
bad singularities whose nature is independent of the potentials. 

The solutions for X > 1 are unphysical due to ( |E.9| ). Therefore all these solutions 
are unacceptable as they interpolate between a bad and a good singularity. They become 
acceptable if we introduce a UV cutoff (or change the UV potential) , and we keep the flow 
to the good IR singularity. 

E.2 Charged solutions, W / 

There is a single acceptable fixed point solution given by 

2(p-l)2 + 1 5-5 2 



X r 



'2(p-l) 1 



P 



(7-5) 



Yr 



P 



( 7 -sy 



W r 



2 + ,5(7 - 5) 
2 + 7(7-5) 



(E.29) 



This is the extremal charged solution described in ( |8.8[ )-( pT2 ) in section || with zero 
temperature. One can also solve for a non-trivial Y keeping X = X c and W = W c describing 
the BH solution: 

1 

- . (E.30) 



Y 



2(p-i) w c r 1+ t^ s) (<P-M , X 2 

V Wc-1 "+" ^c 



e i±fe-l(^ ) _ x 



Finally there is the completely singular generic fixed point X 2 =Y + 1. This exhausts 
all fixed points of the equations. 

E.3 Summary of fixed point solutions and their properties. 

The only acceptable fixed-point (extremal) solutions to equations ( |E.15 )-( [E.17 ) we found 
were the uncharged solution 



X,, 



P-1 

2p 



W u = , Y u = 0. 



(E.31) 



and the charged solution 

X c = 



'2(p-l) 1 



P 



(7-5) 



Y, 



W, 



2(p - 1) 2 + 7 5 - 5 2 
P (7 - 5) 2 

2 + 5(7 - 5) 
2 + 7(7-5) " 



2(p-l) 



P 



(7-5)2 



- >0. 
u 



(E.32) 
(E.33) 
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We can calculate that 



which implies that always 



XT" 1 



2 + 5(j-d) _v 
2 2 



X u 

~Xr 



< 1. 



(E.34) 



(E.35) 



In table |] the position of the fixed points is summarized as a function of the signs of 
two parameters. 



6 


7 — 5 




X u 


Y u 


w u 




x c 


Y c 


w c 


+ 


+ 




+ 










- 


+ 


+ 


- 


- 




- 










+ 


+ 


+ 


- 


+ 




- 










- 


+ 


+ 


+ 


- 




+ 










+ 


+ 


+ 



Table 1: A table showing the positivity properties of the fixed point solutions as a function of the values 
of the parameters. We have used the fact that the charged solutions exist when w p > 0, u > 0, 



v < in (8.8) 



Another important ingredient for the study of the solutions is the behavior of the 
dilaton as well as the potential in such solutions. This is tabulated in table 0. 
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Table 2: We indicate the behavior of the scalar (j> and the potential V for the uncharged and charged 
extremal solutions. An arrow /* indicates an increase towards the IR while \ denotes a decrease 
towards the IR. Note that when 8(^ — 5) > 0, V c decreases towards the IR. Since in the uncharged 
solution Vu always increases towards the IR, this suggests that in these cases </>' and X should 
vanish somewhere in between in order for an interpolation between the uncharged and the charged 
solution to exist, as indicated by the table above. 



E.4 Stability analysis of fixed-point solutions 

We will now proceed to the stability analysis around the fixed-point solutions with or 
without charge. 



E.4.1 Uncharged fixed-point solution 

We start around the uncharged solution, expand W 
and linearize the equations to obtain 

/•\ / \ / (p-l)5 2 -2p 



W u + w,X -+X u + x,Y^-Y u + y 



Mi, 



M„ 
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((p-l)6i-2p)(-y-5) \ 
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(E.36) 
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There are two equal eigenvalues Ai = A2 with respective left eigenvectors v\ given by 



{p - l)5 2 - 2p 
2(p-l)«5 



and a third eigenvalue 



ui 



(p - l)5 2 - 1p 
2(p-l)<J 



t'2 




(E.37) 



2 + 5( 7 - 8) 
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2(p-i) ((p-i)a 2 -2p) 



V 



p (p-l)<5(2 7 -<5)+2(p-2) 
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(E.38) 
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By combining the data of tables [l] and g we obtain that in all cases the uncharged 
solution is unstable to flow away from it, if S 2 < — ^-. To reach this conclusion we must 
take into account whether <f> increases or decreases towards the IR as this crucial for the sign 
of the appropriate eigenvalues. Generically, the endpoint of such a flow is the singularity 
X 2 = Y + 1. With appropriately-tuned initial conditions the charged fixed point solution 
is obtained. 

Finally we would like to point out that we can perturbatively solve the charged equa- 
tions of motion by perturbing in q 2 around the uncharged solution. If we denote collectively 
by $0 the metric components and scalar of the extremal uncharged solution and by $1 the 
first 0(q 2 ) correction, the domain- wall perturbations around the uncharged extremal back- 
ground 



An 



imply that 



(p-l)5 2 ) 



<£ = 3>n 



lnr . 



$1 



50 =0. 



<I>i 



2, 
— lnr . 




2 2v 
q r^ 1 



(E.39) 



(E.40) 



As v < for the existence of acceptable charged extremal solutions, this indicates that 
near the boundary (r — > 00) the charge induces a subleading perturbation to the uncharged 
solution. 



E.4.2 Charged fixed-point solution 

Next we expand and linearize around the charged solution to obtain 




/ 



, M c 



2(p-l)( 7 -A") 








2 % /2p(p-l)( 7 -5) 4 

"Opt' 3 



An obvious left eigenvector and eigenvalue is 



2(p-l)(7-<$) p(7-5) 5 











(E.41) 



Ai 



2(p-l)( 7 -tf) 



Vl 



(E.42) 
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The other two eigenvalues satisfy 



2 w n u , w„u v 



A -7u TT7 K A + 77 h H2=°> ( E - 43 ) 

2(p-l)(7-«J) 4(p-l)( 7 -5) 2 

with solution 

tUpit ± * w 2 u 2 — 4(p — l)w p u 2 v 

A ± = VtkT^) • (E ' 44) 

The charged extremal solutions exist for v < 0, w p > 0, u > 0. As w < 0, A+ has the sign 
of 7 — 5 and A_ the opposite one. 

Similarly, an analysis of the stability of the charged solution taking into account the 
results of table g shows that the solutions have two stable directions and one that is 
unstable. This is in accordance with what we expect from AdS/CFT: we must "tune" one 
initial condition in order to arrive at an acceptable solution in the IR. The tuning allows 
precisely the avoidance of the negative eigenvalue in the IR that drives the solution to be 
bad there. 

Therefore a "well-tuned" solution interpolates between the uncharged extremal solution 
in the UV and the charged extremal solution in the IR. In the next subsection we will express 
the 7<5 = 1 solutions of section as concrete such interpolations. Numerical interpolations 
were found in other concrete cases in |69|], where this issue has been studied in detail, and 
also in [£JJ where numerical interpolation between the charged extremal fixed point in the 
IR and AdS asymptotics in the UV is provided for the case of a flat potential, 5 = 0. 

E.5 The extremal jS = 1 solution in terras of phase variables 

The solution was given in (6.1) in section || Its extremal version reads 



r 2 4(1-<S 2 ) /r \2<5 2 -2 4<5 2 (<S 2 -1) „ 2(6 2 -!) 2 

ds 2 = -—V (3-* 2 )(i+* 2 )di 2 + (-) y(3~* 2 )(i+« 2 ) V + rVc-^X 1 ^) (dx 2 + dy 2 ) , 

(E.45) 
rr \ 3 ~ s2 o a2 v/1 + 5 2 afi- s2 

nr)-i-£) . *'= M p_ P ) ■ < E - 46 » 

r\ 2 <5 ^(^-i) 25 

-) y(3-* 2 ) ( i+^) A= , : Vdt . (E.47) 

The solution can be brought in the domain wall frame r-^uby 

2(g 2 -l) 2 1 AS 2 

e 2A = r 2 y(3-« 2 )d+^) , / = -1/oTP) (E.48) 

(<S 2 -l)(3+,5 2 ) ^ r 

dn = yd+^ 2 )(3-^) (E.49) 

A direct evaluation of the phase variables gives 

jr 2S(S 2 -1) r y> 
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where we used 



1 + (3-<5 2 )(l+<5 2 ) V K + (1+5 2 ) V 1 ^J 

W = ^^ 1 (l-Vf, (E.52) 

rV' o 1 - V 

-^ = (3-5 2 )-^-. (E.53) 

The solution interpolates between (X U ,Y U ,W U ) at V = 1 and (X c , Y^, Wc) at V = 0. The 
relation between V and can be calculated as 

/ r \ 25 _ 2,5 4<S(<S 2 -1) 

e^ 1 = (-^) 0--V) T^^ y(3-^)(i+^> . (E.54) 

F. UV asymptotics 

In this appendix we will investigate the UV asymptotics of the solutions at finite density. 
This information is useful among other things in order to understand the spectra of charged 
fluctuations in such backgrounds. We will remain in p = 3 as this is the only subtle case 
where such an analysis is needed. 



The action we consider is (1.1) with 



V=^-^ + - > Z = l- k -^ + .... (F.l) 

The field 4> corresponds to an operator with UV dimension A so that 

m 2 f = A(A - 3) . (F.2) 

We work in the general diagonal ansatz ( |2.5| ) 

The UV fixed point is given by <j> = 0, and the general finite density solution is the RN 
black-hole given by 

P 2 P 2 P 2 n 2 r A 

D = - 2 f(r) , B = 7W) , C = - 2 , f{r) = l-Mr* + q - w . (F.3) 

We now consider a relevant perturbation associated with the scalar (p. We turn on 
an infinitesimal source for the dual operator and we will work to second order in this 
parameter. 

To leading order in the perturbation the solution is 

D = D (l + x 2 d{r) + 0(x 3 )) , B = B {1 + x 2 b{r) + 0(x 3 )) , (F.4) 

C = C (l + x 2 c(r) + 0{x 3 )) , 4> = xr A g(r), 

where x is a small parameter, the UV coupling of the operator dual to the scalar. The 
functions have a regular expansion around the boundary r = 0: 

b(r) = b + b 3 r 3 + 6 4 r 4 + C(r 6 ) , c{r) = c + c 3 r 3 + c 4 r 4 + O (r 6 ) , (F.5) 
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d{r) = d + d 3 r 3 + d 4 r 4 + 0(r 6 ) , g(r) = 1 + 53 r 3 + 9i r A + C(r 6 



(F.6) 



One of the functions can be trivialized via a coordinate transformation. We will keep it 
explicit and write the general solution as a function of the b(r) coefficients. 
Solving the equations ( |A.1| )-( |A73| ), we obtain: 
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16(2A + l)^ 2 

3M6 63 



.</<> 



_ A(A + 3) 2 M 2 
" 36(2A + 3) ' 

(4A 2 + 6A - 9)M 
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((A 2 + 4A + 2)(A-1) 
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48(2A + 3) 
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64(A + 2)(2A + 1)£ 2 



(F.7) 

(F.8) 
(F.9) 



do 



d\ 



1 _ _6o_ , 

"8 2A ' 3 



(2A + 3)q 2 b 



3(A + l)Mb b 3 

2A(2A + 3) 2A + 3 



(2A 2 + 3A + 9)M 
24(2A + 3) 



+ 



((A 3 + 3A 2 + 10A + 6) - 2(A - 6)k)q 2 



(F.10) 



(F.ll) 



8A(A + 2)£ 2 2(A + 2) 64(A + 2)(2A + l)^ 2 

These formulae are sufficient to compute the UV asymptotics of effective Schrodiger 
potentials. 
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